Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,543,613 A→G 5.1% T167A (ACC→GCC)  fimA → type 1 fimbriae major subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,543,6130AG 5.1% 362.2 / 6.9 120T167A (ACC→GCC) fimAtype 1 fimbriae major subunit
Reads supporting (aligned to +/- strand):  ref base A (58/55);  new base G (3/3);  total (61/59)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.96e-01

Only 100 of 120 total aligned reads displayed.

TAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCG  >  NC_000913/4543465‑4543757
                                                                                                                                                    |                                                                                                                                                
tAGCGCAACAAATGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTTCGACATTTAGTTCAGAAACAACCCTGAATAACTGAACAAATTGCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGTGCCGCAAc                                                                                                                                                 <  2:575288/150‑1 (MQ=255)
 aGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcc                                                                                                                                                >  2:819750/1‑150 (MQ=255)
  gcgcAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAccc                                                                                                                                               >  1:1270913/1‑150 (MQ=255)
   cgcAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGATGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccc                                                                                                                                              >  2:1233008/1‑150 (MQ=255)
   cgcAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccc                                                                                                                                              >  2:1517508/1‑150 (MQ=255)
   cgcAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccc                                                                                                                                              >  1:180562/1‑150 (MQ=255)
   cgcAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAAcccc                                                                                                                                              >  1:43734/1‑150 (MQ=255)
    gcAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCg                                                                                                                                             >  1:757091/1‑150 (MQ=255)
    gcAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCg                                                                                                                                             >  1:751552/1‑150 (MQ=255)
    gcAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCg                                                                                                                                             >  2:1397275/1‑150 (MQ=255)
    gcAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCg                                                                                                                                             >  1:1446836/1‑150 (MQ=255)
     caacaaACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCgg                                                                                                                                            >  1:537995/1‑150 (MQ=255)
      aacaaACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCggg                                                                                                                                           <  2:147023/150‑1 (MQ=255)
       acaaACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGt                                                                                                                                          <  1:819750/150‑1 (MQ=255)
       acaaACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGt                                                                                                                                          <  1:1290220/150‑1 (MQ=255)
         aaaCGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAATCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCACCCGGGGTTGCGGCCCCGGTtgc                                                                                                                                        <  2:124724/150‑1 (MQ=255)
         aaaCGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCACCCGGGGTTGCGGCCCCGGTtgc                                                                                                                                        <  1:330782/150‑1 (MQ=255)
         aaaCGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGtgc                                                                                                                                        >  2:435749/1‑150 (MQ=255)
           aCGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGtgctg                                                                                                                                      <  2:1407415/150‑1 (MQ=255)
            cGTTGGTGTGCAGAGCCTGGACAGCACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTCGTTCAGAAACAACCCTGACTAACGGCACCACTACCATTCCGTTCCCGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGtgctgc                                                                                                                                     <  1:182653/150‑1 (MQ=255)
             gTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTgctgct                                                                                                                                    <  2:179020/150‑1 (MQ=255)
              ttGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTa                                                                                                                                   <  2:1473637/150‑1 (MQ=255)
                ggTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTGAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAAt                                                                                                                                 <  2:862430/150‑1 (MQ=255)
                 gtgtGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCTTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATg                                                                                                                                <  2:164195/150‑1 (MQ=255)
                  tgtgCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGc                                                                                                                               >  1:79704/1‑150 (MQ=255)
                   gtgCAGATCCTGAACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCg                                                                                                                              <  1:1129082/150‑1 (MQ=255)
                    tgCAGATCCTGGACAGAACGGGTGCTGCACTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCgg                                                                                                                             >  2:1445216/1‑150 (MQ=255)
                     gCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCACCCGGGGTTGCGGCCCCGGTtgcaa                                                                                                                                      >  2:1114746/1‑138 (MQ=255)
                     gCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCACCCGGGGTTGCGGCCCCGGTtgcaa                                                                                                                                      <  1:1114746/140‑3 (MQ=255)
                      cAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGAt                                                                                                                           >  1:417518/1‑150 (MQ=255)
                       aGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATg                                                                                                                          <  2:324755/150‑1 (MQ=255)
                       aGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATg                                                                                                                          >  1:336203/1‑150 (MQ=255)
                       aGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATg                                                                                                                          <  1:1310236/150‑1 (MQ=255)
                           ccTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAc                                                                                                                      <  1:813251/150‑1 (MQ=255)
                           ccTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAc                                                                                                                      >  1:686494/1‑150 (MQ=255)
                           ccTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAc                                                                                                                      >  2:1363205/1‑150 (MQ=255)
                           ccTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAc                                                                                                                      <  2:985310/150‑1 (MQ=255)
                           ccTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGAc                                                                                                                      <  1:191100/150‑1 (MQ=255)
                               gACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTc                                                                                                                  >  1:192101/1‑150 (MQ=255)
                                      cGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTc                                                                                                           <  2:180562/150‑1 (MQ=255)
                                         gTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGt                                                                                                        <  1:495402/150‑1 (MQ=255)
                                            ctgcGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGTGACCTTCAAGGTTCAGTATc                                                                                                     >  1:48916/1‑150 (MQ=255)
                                            ctgcGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATc                                                                                                     <  1:24915/150‑1 (MQ=255)
                                            ctgcGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATc                                                                                                     <  2:1226100/150‑1 (MQ=255)
                                               cgcTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAAt                                                                                                  <  2:746094/150‑1 (MQ=255)
                                                 cTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAgg                                                                                                              >  2:883800/1‑136 (MQ=255)
                                                  tGATGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAAc                                                                                               <  1:1233008/150‑1 (MQ=255)
                                                        tGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTAccc                                                                                         <  2:537995/150‑1 (MQ=255)
                                                            tGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGt                                                                                     >  1:432213/1‑150 (MQ=255)
                                                                gCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCGCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAg                                                                                 >  2:1073950/1‑150 (MQ=255)
                                                                gCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCATGTTCAg                                                                                 >  1:1356974/1‑150 (MQ=255)
                                                                gCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAg                                                                                 >  1:10131/1‑150 (MQ=255)
                                                                gCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAg                                                                                 >  1:103634/1‑150 (MQ=255)
                                                                 cGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAgg                                                                                <  1:109635/150‑1 (MQ=255)
                                                                  gACATTTAGTTCAGAAACAACCCTGAATAACGGAACCATTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCACAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAggg                                                                               <  2:1449251/150‑1 (MQ=255)
                                                                  gACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAggg                                                                               <  2:1356974/150‑1 (MQ=255)
                                                                  gACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAggg                                                                               <  2:130150/150‑1 (MQ=255)
                                                                        tAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCa                                                                         >  2:1174058/1‑150 (MQ=255)
                                                                        tAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCa                                                                         >  1:1185828/1‑150 (MQ=255)
                                                                         atttCATAAACAACCCTGAATAACGGAACCAATACCATTCCTTTCCAGGCGCGTTATTTTGCAACAGGGGCCGCAACCCCGGGTGCTGCTAATTCGGATTCGACCTTCAATGTTCATTATAAATAACCTACCCAGGTTCAGGTACGTCAt                                                                        <  2:1270913/148‑1 (MQ=255)
                                                                             cAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACg                                                                    >  2:984238/1‑150 (MQ=255)
                                                                              aGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCa                                                                                                          <  2:925609/111‑1 (MQ=255)
                                                                              aGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCa                                                                                                          >  1:925609/1‑111 (MQ=255)
                                                                              aGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACgg                                                                   >  1:1442616/1‑150 (MQ=255)
                                                                                   caacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAgg                                                              <  2:757091/150‑1 (MQ=255)
                                                                                    aacCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAgtg                                                             <  2:1442616/150‑3 (MQ=255)
                                                                                      cccTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGAt                                                           >  2:828177/1‑150 (MQ=255)
                                                                                        cTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAgg                                                                                                              >  2:875595/1‑97 (MQ=255)
                                                                                        cTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAgg                                                                                                              <  1:875595/97‑1 (MQ=255)
                                                                                        cTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGc                                                         <  1:1445216/150‑1 (MQ=255)
                                                                                         tGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGcc                                                        <  2:1109266/150‑1 (MQ=255)
                                                                                          gAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGccc                                                       <  2:1352666/150‑1 (MQ=255)
                                                                                               aCGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCt                                                  <  2:614644/150‑1 (MQ=255)
                                                                                                   aaCCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTtgtg                                              >  1:417061/1‑150 (MQ=255)
                                                                                                     ccAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCg                                            >  2:838254/1‑150 (MQ=255)
                                                                                                      cAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGa                                           >  2:71782/1‑150 (MQ=255)
                                                                                                       aaTACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTCTCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGAt                                          >  2:1001730/1‑150 (MQ=255)
                                                                                                          aCCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATaaa                                       <  1:1174058/150‑1 (MQ=255)
                                                                                                           ccATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATaaaa                                      >  2:1006097/1‑150 (MQ=255)
                                                                                                           ccATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATaaaa                                      >  1:532972/1‑150 (MQ=255)
                                                                                                           ccATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCACACCCTTGTGCGATaaaa                                      >  1:296940/1‑150 (MQ=255)
                                                                                                                ccGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAAc                                 >  1:1029467/1‑150 (MQ=255)
                                                                                                                     ccAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGa                            <  1:912159/150‑1 (MQ=255)
                                                                                                                      cAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGaa                           >  1:415377/1‑150 (MQ=255)
                                                                                                                        ggCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGaaaa                         >  1:564286/1‑150 (MQ=255)
                                                                                                                             gTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAag                    >  2:630513/1‑150 (MQ=255)
                                                                                                                               tATTTTGCAACCGGGGCCGCAACCGCGGGTGCTGCTAATGCGGATGCGACCTTCTAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAagag                  <  2:1419754/150‑1 (MQ=255)
                                                                                                                                  gcagCACCCGGGGTTGCGGCCCCGGTTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGatt               >  1:756686/4‑150 (MQ=255)
                                                                                                                                   ttGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGatta              <  2:686494/150‑1 (MQ=255)
                                                                                                                                   cagCACCCGGGGTTGCGGCCCCGGTTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGattt              >  1:170928/3‑149 (MQ=255)
                                                                                                                                    tGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCCTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGattat             >  2:629262/1‑150 (MQ=255)
                                                                                                                                      cAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTAttt           >  1:1193106/1‑150 (MQ=255)
                                                                                                                                       aaCCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTc          >  1:877070/1‑150 (MQ=255)
                                                                                                                                       aaCCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTc          <  2:1208666/150‑1 (MQ=255)
                                                                                                                                         ccGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGAGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTa        <  1:828177/150‑1 (MQ=255)
                                                                                                                                           ggggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTAtt      >  1:227551/1‑150 (MQ=255)
                                                                                                                                             ggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAg    <  1:146499/150‑1 (MQ=255)
                                                                                                                                             ggCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGGCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAg    <  2:159740/150‑1 (MQ=255)
                                                                                                                                              gCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAg  >  1:1433202/1‑149 (MQ=255)
                                                                                                                                               ccGCATCCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGTCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGcg  <  1:1174017/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                
TAGCGCAACAAACGTTGGTGTGCAGATCCTGGACAGAACGGGTGCTGCGCTGACGCTGGATGGTGCGACATTTAGTTCAGAAACAACCCTGAATAACGGAACCAATACCATTCCGTTCCAGGCGCGTTATTTTGCAACCGGGGCCGCAACCCCGGGTGCTGCTAATGCGGATGCGACCTTCAAGGTTCAGTATCAATAACCTACCCAGGTTCAGGGACGTCATTACGGGCAGGGATGCCCACCCTTGTGCGATAAAAATAACGATGAAAAGGAAGAGATTATTTCTATTAGCG  >  NC_000913/4543465‑4543757

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