Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 987,561 G→C 9.5% intergenic (‑579/+24) ompF ← / ← asnS outer membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913987,5610GC 9.5% 306.4 / 19.9 95intergenic (‑579/+24)ompF/asnSouter membrane porin 1a (Ia;b;F)/asparaginyl tRNA synthetase
Reads supporting (aligned to +/- strand):  ref base G (50/36);  new base C (5/4);  total (55/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.46e-01

AATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACC  >  NC_000913/987417‑987707
                                                                                                                                                 |                                                                                                                                                  
aaTAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGgcgc                                                                                                                                                 <  2:37903/148‑1 (MQ=255)
    aGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGct                                                                                                                                            <  3:706773/149‑1 (MQ=255)
     gAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGCCCGCTGGct                                                                                                                                            <  1:563361/148‑1 (MQ=255)
        tGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCtt                                                                                                                                         >  1:459416/1‑148 (MQ=255)
         gCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTg                                                                                                                                        <  2:451367/148‑1 (MQ=255)
           ggATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCATCGGGCGCTGGCTCTTTAAg                                                                                                                                     >  4:650998/1‑149 (MQ=255)
            gATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTGGACAGAACAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGa                                                                                                                                    <  4:632143/149‑1 (MQ=255)
                             aaaTGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTg                                                                                                                   >  3:574650/1‑149 (MQ=255)
                                  ttGCAATGTAATATTCCGTAACAGGATGATCGTTATGGTGGGGTTGCAAATAGATTGCTGGACTTAACAAACAGATACCTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGGCGCTAAATTACAAGCTGGCGt                                                                                                               <  2:15898/148‑1 (MQ=255)
                                    gCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTAc                                                                                                            <  4:822400/149‑1 (MQ=255)
                                      aaTGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTAcgc                                                                                                          <  4:486758/149‑1 (MQ=255)
                                       aTGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTAcgcg                                                                                                         >  3:361716/1‑149 (MQ=255)
                                        tGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGAAAAAAGATTGCTTGTCCTAACAAAGAGATAAATTGACAGAAAACGTAAAAAAAGAGACAGCGGGCGCCGTCCCCTGGCGACGATAAAATAGAAGCTAGCGTTAAGCgg                                                                                                        >  4:477850/1‑149 (MQ=255)
                                         gTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGa                                                                                                       <  4:2252/149‑1 (MQ=255)
                                               ttCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACg                                                                                                  >  1:423442/1‑148 (MQ=255)
                                                tCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTTATAAAGAGATAACTTGACGAAAAAGATAAAAAAAAGGGCAGCGGGCGCATGCTGTTGAAGAAGCGAAATTTGAAACTGACTTTTCGAGGAGTAGGt                                                                                                 >  1:462555/1‑148 (MQ=255)
                                                tCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTg                                                                                                <  4:361716/149‑1 (MQ=255)
                                                tCCGTAACAGGATGATCGTCATGGCGGGGTTGCAAATAGATTGCTGGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGt                                                                                                 <  2:342008/148‑1 (MQ=255)
                                                 ccGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTgg                                                                                               <  3:415827/149‑1 (MQ=255)
                                                 ccGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTg                                                                                                <  2:357909/148‑1 (MQ=255)
                                                   gTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGa                                                                                             >  3:917304/1‑149 (MQ=255)
                                                    tAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGaa                                                                                            <  4:874461/149‑1 (MQ=255)
                                                     aaCAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAAc                                                                                           <  4:917304/149‑1 (MQ=255)
                                                      aCAAGAAGATCGGTAATGCGGGGTTTCAAATTAGATGCTTGACCTAACAAAAACAATACTTGACAGAAAAGATAAAAAAAGAGCCAGCCGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAAc                                                                                           <  2:811895/148‑1 (MQ=255)
                                                         ggATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGa                                                                                                                                    <  1:108497/104‑1 (MQ=255)
                                                         ggATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACCTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGa                                                                                                                                    >  2:108497/1‑104 (MQ=255)
                                                          gatgatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGTAGTACGTGGGAACGGAat                                                                                      <  3:756925/149‑1 (MQ=255)
                                                          gatgatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAat                                                                                      >  4:638141/1‑149 (MQ=255)
                                                          gatgatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAat                                                                                      >  3:248855/1‑149 (MQ=255)
                                                           atgatCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatc                                                                                     <  3:122872/149‑1 (MQ=255)
                                                                cGTTATGGCGGGGTTGCAAATAGATTGCTTGGCTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacat                                                                                 >  2:527199/1‑148 (MQ=255)
                                                                cGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTg                                                                                                                                               >  2:385788/1‑86 (MQ=255)
                                                                cGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTg                                                                                                                                               <  1:385788/86‑1 (MQ=255)
                                                                 gTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacatca                                                                               <  3:581922/149‑1 (MQ=255)
                                                                   tATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAatcacatcac                                                                              >  1:540723/1‑148 (MQ=255)
                                                                       gCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgtac                                                                          >  1:788452/1‑148 (MQ=255)
                                                                        cGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCacgtacgt                                                                        <  3:358835/149‑1 (MQ=255)
                                                                          gggTTGCAAATAGATTGCTTGACTTAACAAACAGATACCTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGtt                                                                       <  1:527199/148‑1 (MQ=255)
                                                                             ttGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGc                                                                   >  3:678905/1‑149 (MQ=255)
                                                                              tGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGc                                                                   >  1:673102/1‑148 (MQ=255)
                                                                                cAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACg                                                                >  3:638426/1‑149 (MQ=255)
                                                                                 aaaTAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGc                                                               >  3:377048/1‑149 (MQ=255)
                                                                                   aTAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCa                                                             <  3:477850/149‑1 (MQ=255)
                                                                                        aTTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGtt                                                          >  1:413145/1‑148 (MQ=255)
                                                                                          tGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACg                                                       >  3:443195/1‑149 (MQ=255)
                                                                                             ttGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTa                                                     <  1:904575/148‑1 (MQ=255)
                                                                                              tGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTaa                                                    >  1:239277/1‑148 (MQ=255)
                                                                                              tGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTaa                                                    >  1:614670/1‑148 (MQ=255)
                                                                                                  ttAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCaa                                                >  1:141943/1‑148 (MQ=255)
                                                                                                  ttAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCaa                                                >  2:433441/1‑148 (MQ=255)
                                                                                                  ttAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAAt                                               >  3:262224/1‑149 (MQ=255)
                                                                                                  ttAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAAt                                               >  4:707738/1‑149 (MQ=255)
                                                                                                   tagcaaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAAt                                               >  1:833977/4‑148 (MQ=255)
                                                                                                   tAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAAt                                               >  2:183901/1‑148 (MQ=255)
                                                                                                      caaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAg                                            <  2:417476/148‑1 (MQ=255)
                                                                                                       aaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGAc                                          >  4:190345/1‑149 (MQ=255)
                                                                                                       aaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGAc                                          >  3:732330/1‑149 (MQ=255)
                                                                                                       aaacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGATCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGAc                                          >  4:888487/1‑149 (MQ=255)
                                                                                                        aacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGAc                                          >  2:215232/1‑148 (MQ=255)
                                                                                                        aacaGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGAc                                          >  2:638698/1‑148 (MQ=255)
                                                                                                          caGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGtt                                       <  4:248855/149‑1 (MQ=255)
                                                                                                            gATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCtt                                                                                                                                         <  3:492495/49‑1 (MQ=255)
                                                                                                            gATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCtt                                                                                                                                         >  4:492495/1‑49 (MQ=255)
                                                                                                            gATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTAc                                                                                                   >  3:710826/1‑87 (MQ=255)
                                                                                                            gATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTAc                                                                                                   <  4:710826/87‑1 (MQ=255)
                                                                                                            gATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTc                                      >  2:662470/1‑148 (MQ=255)
                                                                                                              tAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTa                                                                                                    <  2:379144/84‑1 (MQ=255)
                                                                                                              tAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTa                                                                                                    >  1:379144/1‑84 (MQ=255)
                                                                                                                aCTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCaaaa                                  <  2:413145/148‑1 (MQ=255)
                                                                                                                  ttGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCa                               >  3:567550/1‑149 (MQ=255)
                                                                                                                   tGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCaa                              <  4:567550/149‑1 (MQ=255)
                                                                                                                    gACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAg                             <  3:743119/149‑1 (MQ=255)
                                                                                                                    gACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAg                             >  4:141739/1‑149 (MQ=255)
                                                                                                                        gAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGAcc                          >  1:938107/1‑148 (MQ=255)
                                                                                                                        gAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCg                         >  3:47634/1‑149 (MQ=255)
                                                                                                                         aaaaGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGa                        >  3:687263/1‑149 (MQ=255)
                                                                                                                           aaGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGaaa                      <  4:638426/149‑1 (MQ=255)
                                                                                                                           aaGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGaa                       <  2:703978/148‑1 (MQ=255)
                                                                                                                               tAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGa                            >  3:227196/1‑139 (MQ=255)
                                                                                                                               tAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGa                            <  4:227196/139‑1 (MQ=255)
                                                                                                                                tcttcaaGAGCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTg                  >  2:728777/6‑148 (MQ=255)
                                                                                                                                aaaaaaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTg                  >  1:484455/1‑148 (MQ=255)
                                                                                                                                 aaaaaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGaa                >  3:89543/1‑149 (MQ=255)
                                                                                                                                 aaaaaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGa                 <  2:239277/148‑1 (MQ=255)
                                                                                                                                  aaaaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGaa                >  1:824160/1‑148 (MQ=255)
                                                                                                                                   aaaaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATg              >  3:506170/1‑149 (MQ=255)
                                                                                                                                     aaGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGc             <  1:433441/148‑1 (MQ=255)
                                                                                                                                     aaGAGCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGATCGAAACCTGAATGCg            >  4:842096/1‑149 (MQ=255)
                                                                                                                                      agagCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAAc                                                                                           <  4:126499/69‑1 (MQ=38)
                                                                                                                                      agagCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAAc                                                                                           >  3:126499/1‑69 (MQ=38)
                                                                                                                                      agagCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCaaaa                                  <  3:909678/126‑1 (MQ=255)
                                                                                                                                      agagCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCaaaa                                  >  4:909678/1‑126 (MQ=255)
                                                                                                                                        agCCAGCGCCCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGa          >  1:806809/1‑148 (MQ=255)
                                                                                                                                          ccAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAAc        >  2:688505/1‑148 (MQ=255)
                                                                                                                                          ccAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACa       <  3:378578/149‑1 (MQ=255)
                                                                                                                                           cAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAg      >  4:947860/1‑149 (MQ=255)
                                                                                                                                               gggCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTAcc  <  4:678905/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                  
AATAAGAATGCGGATTAGTTTTTGCGCGTAAATGTTGCAATGTAATATTCCGTAACAGGATGATCGTTATGGCGGGGTTGCAAATAGATTGCTTGACTTAACAAACAGATAACTTGACAGAAAAGATAAAAAAAGAGCCAGCGGGCGCTGGCTCTTGAAGACGCTAAATTAGAAGCTGGCGTTACGCGGAGTACGTGGGAACGGAATCACATCACGTACGTTTTGCACGCCAGTTACGTAAGCAATCAGACGTTCAAAACCAAGACCGAAACCTGAATGCGGAACAGTACC  >  NC_000913/987417‑987707

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 30 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: