Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 922,314:1 +C 6.3% intergenic (+21/+52) macB → / ← cspD macrolide ABC transporter peremase/ATPase/inhibitor of DNA replication, cold shock protein homolog

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913922,3141.C 6.3% 392.9 / 23.4 96intergenic (+21/+52)macB/cspDmacrolide ABC transporter peremase/ATPase/inhibitor of DNA replication, cold shock protein homolog
Reads supporting (aligned to +/- strand):  ref base . (46/44);  new base C (3/3);  total (49/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.20e-01

GTTTTTCACCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTT  >  NC_000913/922169‑922463
                                                                                                                                                  |                                                                                                                                                     
gTTTTTCACCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGTACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGa                                                                                                                                                   >  2:30869/1‑150 (MQ=255)
gTTTTTCACCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGa                                                                                                                                                   >  1:611121/1‑150 (MQ=255)
      cACCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCCG‑CCGATCGGGc                                                                                                                                             >  2:451046/1‑150 (MQ=255)
      cACCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGc                                                                                                                                             >  1:1690458/1‑150 (MQ=255)
      cACCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGc                                                                                                                                             >  1:635403/1‑150 (MQ=255)
       aCCGTTGGCGGTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGAGGCGCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCt                                                                                                                                            <  1:1379269/150‑1 (MQ=255)
         cGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATTCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTgg                                                                                                                                          <  1:1234853/150‑1 (MQ=255)
            tGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCAGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCAt                                                                                                                                       <  2:714516/150‑1 (MQ=255)
             ggCGCTGCTGCGGGCGTTTCTCTGCTCGACGGTCACCGTGAGTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATTCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCAtt                                                                                                                                      <  1:114711/150‑1 (MQ=255)
              gcgcTGCTGCTGGCGTTTCTCTGCTCGACGGTCACGGGTATTTTATTTGGCTGGTTACCCGCACTAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAGCCCGATC‑GGCTGGCAttt                                                                                                                                     <  2:2009904/150‑1 (MQ=255)
               ggctgctgctGGCGTTTCTCTGCTCGACGGTCAGCGTTATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCTGCCCGATC‑GGCTGGCAtttt                                                                                                                                    <  1:674383/149‑1 (MQ=255)
                 ctgctgctgGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTCGCTGGTTACCCTCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGc                                                                                                                                  >  1:1255652/1‑150 (MQ=255)
                      gctgGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATTCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTa                                                                                                                             >  2:1103671/1‑150 (MQ=255)
                       ctgGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGGTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAg                                                                                                                            <  1:406835/150‑1 (MQ=255)
                         ggCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGa                                                                                                                          <  1:2055401/150‑1 (MQ=255)
                         ggCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGAGCGGGCTGGCATTTTGCCTTTAGGa                                                                                                                          >  2:499035/1‑150 (MQ=255)
                          gCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGAt                                                                                                                         >  1:1682186/1‑150 (MQ=255)
                           cGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATg                                                                                                                        <  1:1801442/150‑1 (MQ=255)
                            gTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGt                                                                                                                       <  2:969895/150‑1 (MQ=255)
                              ttCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTac                                                                                                                     <  1:942107/150‑1 (MQ=255)
                              ttCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCAAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTac                                                                                                                     >  1:1619203/1‑150 (MQ=255)
                                 tctGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTAcaca                                                                                                                  <  2:1735302/150‑1 (MQ=255)
                                 tctGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGATGGCATTTTGCCTTTAGGATGTAcaca                                                                                                                  >  2:1069158/1‑150 (MQ=255)
                                   tGCGCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAAt                                                                                                                <  1:464541/150‑1 (MQ=255)
                                     cTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATga                                                                                                              <  2:1706692/150‑1 (MQ=255)
                                                    ttttttATTTGGCGTGTTACCCGCTCGAAATGCGGCACGACTGGAGCCAGTCGATGCTCTGGCAAGAGAGTAATTTTTGAGTTAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATg                                                                                               <  1:210495/148‑1 (MQ=255)
                                                               gCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCt                                                                                    >  1:1057749/1‑150 (MQ=255)
                                                               gCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCt                                                                                    >  2:176953/1‑150 (MQ=255)
                                                                 gggTTACCCTCACGAAATGCGGCACGCCTGGCTCCAGTAGATGCGCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTc                                                                                  <  2:1756088/149‑1 (MQ=255)
                                                                  ggTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGTTCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTACt                                                                                 >  2:1282656/1‑150 (MQ=255)
                                                                  ggTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGCGCTATGCGACTGCCGCTTCt                                                                                 >  1:2025226/1‑150 (MQ=255)
                                                                  ggTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCt                                                                                 >  2:625679/1‑150 (MQ=255)
                                                                    ttACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTAc                                                                               <  1:653224/150‑1 (MQ=255)
                                                                     tACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACt                                                                              <  1:1275652/150‑1 (MQ=255)
                                                                       cccGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCCATCGGGCTGGATTTTTCCATTAAGGATGTACAAAATCAGAAATAATAGATATGATACTGCCCCttgaactac                                                                            >  2:1284404/1‑143 (MQ=255)
                                                                         cGCACGAAATGCTGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑ACTATCGGGCTGGCAGATTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGa                                                                          >  1:304833/1‑150 (MQ=255)
                                                                         cGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGa                                                                          >  1:1674910/1‑150 (MQ=255)
                                                                         cGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGa                                                                          >  2:1788816/1‑150 (MQ=255)
                                                                            aCGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACgg                                                                       >  1:2083273/1‑150 (MQ=255)
                                                                             cGAAAGGCGGCACGAGTGGATCCAGTAGATGCTCTGGGACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACTATGAGACAGAAGAGCTAGGCGACTGCCGCTTCTACTTCGACggg                                                                      <  1:1381163/150‑1 (MQ=255)
                                                                               aaaTGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGca                                                                    <  2:644799/150‑1 (MQ=255)
                                                                                aaTGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGAACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGcac                                                                   <  1:1375248/150‑1 (MQ=255)
                                                                                 aTGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGcaca                                                                  <  1:262431/150‑1 (MQ=255)
                                                                                     ggCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGAATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATaa                                                              >  1:741381/1‑150 (MQ=255)
                                                                                      gCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAac                                                             >  1:802431/1‑150 (MQ=255)
                                                                                      gCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAac                                                             >  1:868580/1‑150 (MQ=255)
                                                                                       cACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGTGCACAATAaca                                                            >  1:1896713/1‑150 (MQ=255)
                                                                                       cACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAaca                                                            >  1:1905425/1‑150 (MQ=255)
                                                                                           aCTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTgg                                                        <  1:872608/150‑1 (MQ=255)
                                                                                           aCTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTgg                                                        >  1:1206965/1‑150 (MQ=255)
                                                                                             tGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCg                                                      <  1:599791/150‑1 (MQ=255)
                                                                                             tGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCg                                                      <  2:628029/150‑1 (MQ=255)
                                                                                             tGGATCCAGTAGATGCTCAGGCACGAGAGTAATTTTTGATATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCg                                                      <  2:1083403/150‑1 (MQ=255)
                                                                                              ggATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGt                                                     <  1:857142/150‑1 (MQ=255)
                                                                                               gATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTg                                                    >  1:652051/1‑150 (MQ=255)
                                                                                               gATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTg                                                    <  1:1841704/150‑1 (MQ=255)
                                                                                               gATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTg                                                    >  2:69264/1‑150 (MQ=255)
                                                                                                aTCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCAGTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGa                                                   >  1:248516/1‑150 (MQ=255)
                                                                                                 tCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGAt                                                  <  1:1765942/150‑1 (MQ=255)
                                                                                                  ccAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAATAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGAtt                                                 <  2:1255652/150‑1 (MQ=255)
                                                                                                   cAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGACGCTTCTAATACGACGGCCAAAATAACACTGGCGTGATTg                                                >  1:136280/1‑150 (MQ=255)
                                                                                                    aGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGc                                               <  1:1060930/150‑1 (MQ=255)
                                                                                                     gTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGcc                                              <  1:176953/150‑1 (MQ=255)
                                                                                                             tctGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGcc                                      >  1:759070/1‑150 (MQ=255)
                                                                                                               tGGCACGAGAGTAATTTTTGAGATAAAAATGCAAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGcccc                                    <  2:1619203/150‑1 (MQ=255)
                                                                                                               tGGCACGAGAGTAAGTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGcccc                                    >  1:897094/1‑150 (MQ=255)
                                                                                                                 gCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCtg                                  <  2:494575/150‑1 (MQ=255)
                                                                                                                 gCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCtg                                  >  1:1928559/1‑150 (MQ=255)
                                                                                                                  cACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCtgg                                 <  2:1858864/150‑1 (MQ=255)
                                                                                                                     gagaGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGGGATTGCCTTTTGGCCCCtggtgg                              <  1:1522673/150‑1 (MQ=255)
                                                                                                                      agagTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGa                             <  1:561603/150‑1 (MQ=255)
                                                                                                                       tagTAATTTTTGAGTTAAAAATGCCCG‑CGGATCTGGCTGGCATTTTGCCTTTAGGAGGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGAc                            <  2:1647422/149‑1 (MQ=255)
                                                                                                                       gagTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGAc                            >  2:1927384/1‑150 (MQ=255)
                                                                                                                       gagTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGAc                            >  2:1063888/1‑150 (MQ=255)
                                                                                                                        agTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACa                           <  2:1617016/150‑1 (MQ=255)
                                                                                                                         gTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACAt                          >  2:665184/1‑150 (MQ=255)
                                                                                                                            agttttGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCGTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGCCTGACGCTACTACTTCGACGGGCACAATAACCCTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCaa                       >  1:205583/3‑150 (MQ=255)
                                                                                                                            aTTTTTGATATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGACTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCaa                       <  2:1308993/150‑1 (MQ=255)
                                                                                                                             tttttGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCaaa                      >  1:1917255/1‑150 (MQ=255)
                                                                                                                              ttttGAGATAATAATGCCAG‑CCGATCGGGCTGGCATTTTGACTTTAGGATGTACACAATGAGACATAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGTTTGCCTTTTGGCCCCTGGTGGACATCAAAc                     <  2:1069860/150‑1 (MQ=255)
                                                                                                                                  gagaTAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGaa                 >  2:1870951/1‑150 (MQ=255)
                                                                                                                                  gagaTAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGaa                 <  1:1103671/150‑1 (MQ=255)
                                                                                                                                   agaTAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGGCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACATCGACGGGCACAATAACACTGGCGTGATTGCCTTTTCGCCGCTGGTGGACATCCCACTGAAc                >  1:620442/1‑150 (MQ=255)
                                                                                                                                      gcaaaaTGCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGa             >  2:1619006/3‑150 (MQ=255)
                                                                                                                                        aaaaTGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCAGAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGAtt           <  2:1248338/150‑1 (MQ=255)
                                                                                                                                        aaaaTGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGAtt           >  1:762968/1‑150 (MQ=255)
                                                                                                                                         aaaTGCCAGCCCGATC‑GGCTGGCATATTGCCGTTAGGAGGGTCACAATGAGACAGAAGATATATTAGACTGCCGCTTCTACTTCGACCTGCAAATTAAAACTCTAGTGATTGGCTTTTGGCGCCTGGTGGACATCAACCGGAACGGATTg          >  2:398792/1‑150 (MQ=255)
                                                                                                                                         aaaTGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTCTGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTg          <  1:744051/150‑1 (MQ=255)
                                                                                                                                         aaaTGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTg          <  2:1682186/150‑1 (MQ=255)
                                                                                                                                         aaaTGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTg          <  2:840437/150‑1 (MQ=255)
                                                                                                                                            tGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTcc       >  1:1669765/1‑150 (MQ=255)
                                                                                                                                            tGCCAG‑CCGATCGGGCTGGCATTTTGCCATTACGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGAGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTcc       <  1:16453/150‑1 (MQ=255)
                                                                                                                                             gCCAGCCCGATC‑GGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCa      >  1:1878085/1‑150 (MQ=255)
                                                                                                                                              ccAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTAGGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGACTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCGAg     >  1:367137/1‑150 (MQ=255)
                                                                                                                                                aG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAATAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCATTTTGGCCCCTGGTGTACATCAAACTGAACGGATTGTCCAGCt   >  1:1288792/1‑150 (MQ=255)
                                                                                                                                                aG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCt   >  1:159652/1‑150 (MQ=255)
                                                                                                                                                 g‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCGACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACTTCAAACTGAACGGATTGTCCAGCtt  >  2:1029875/1‑150 (MQ=255)
                                                                                                                                                  |                                                                                                                                                     
GTTTTTCACCGTTGGCGCTGCTGCTGGCGTTTCTCTGCTCGACGGTCACCGGGATTTTATTTGGCTGGTTACCCGCACGAAATGCGGCACGACTGGATCCAGTAGATGCTCTGGCACGAGAGTAATTTTTGAGATAAAAATGCCAG‑CCGATCGGGCTGGCATTTTGCCTTTAGGATGTACACAATGAGACAGAAGAGCTATGCGACTGCCGCTTCTACTTCGACGGGCACAATAACACTGGCGTGATTGCCTTTTGGCCCCTGGTGGACATCAAACTGAACGGATTGTCCAGCTT  >  NC_000913/922169‑922463

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑