Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,131,375 A→G 5.3% M275T (ATG→ACG)  fimD ← type I pili subunit FimD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,131,3750AG 5.3% 218.6 / 2.1 75M275T (ATG→ACG) fimDtype I pili subunit FimD
Reads supporting (aligned to +/- strand):  ref base A (33/37);  new base G (2/2);  total (35/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.06e-01

TCAGGTCGTCGATGGCGTAGGGGCCGGCGCTGACGTACGTCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTG  >  NC_002947/2131229‑2131523
                                                                                                                                                  |                                                                                                                                                    
tCAGGTCGTCGATGGCGTAGGGGCCGGCGCTGACGTACGTCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATg                                                                                                                                                    >  8:233981/1‑149 (MQ=255)
   ggTCGTCGATGGCGTAGGGGCCGGCGCTGACGTACGTCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTcc                                                                                                                                                 >  6:207889/1‑149 (MQ=255)
     tcgtcgATGGCGTAGGGGCCGGCGCTGACGTACGTCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGa                                                                                                                                               <  3:90479/149‑1 (MQ=255)
         cgctGGCGTAGGGGCCGGCGCTGCCGTACGTCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCCCGCCGCGAATGATCGGGGCATAACCTCGCTGGTTGTCGGCCCGCATGCACAGGTACGAGGcc                                                                                                                                           >  8:230216/4‑149 (MQ=255)
         cgATGGCGTAGGGGCCGGCGCTGACGTACGTCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGcc                                                                                                                                           <  7:360645/149‑1 (MQ=255)
                  aGGGGCGGGCGCTGACGTCCGGCGAGGAGCGCGGATAACCATGCTGCCAGCTCTCGCCTTTCGCCCGGCGCTGGGCCACGCCGCGAAGGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACg                                                                                                                                  <  7:230216/149‑1 (MQ=255)
                    gggCCGGCGCTGACGTACGTCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGcc                                                                                                                                <  3:68039/149‑1 (MQ=255)
                     ggCCGGCGCTGACGTACGTCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCt                                                                                                                               <  3:333319/149‑1 (MQ=255)
                      gCCGGCGCTGACGTACGTCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCt                                                                                                                               >  1:57139/1‑148 (MQ=255)
                      gCCGGCGCTGACGTACGTCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTg                                                                                                                              <  5:327412/149‑1 (MQ=255)
                        cGGCGCTGACGTACGTCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGa                                                                                                                             >  3:154297/1‑148 (MQ=255)
                         ggCGCTGACGTACGTCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGat                                                                                                                            <  6:87647/148‑1 (MQ=255)
                         ggCGCTGACGTACGTCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGat                                                                                                                            >  5:87647/1‑148 (MQ=255)
                           ctcTGACGTACGTCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGTGGCATAACCGCGCTGGTTGTCTGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATc                                                                                                                         <  7:187730/147‑1 (MQ=255)
                           cgcTGACGTACGTCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCATCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATc                                                                                                                         <  6:38951/149‑1 (MQ=255)
                            gcTGACGTACGTCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCg                                                                                                                        >  1:360407/1‑149 (MQ=255)
                                  gtacgtCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGATTGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCAcc                                                                                                                   >  4:237839/1‑148 (MQ=255)
                                  gtacgtCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCAcc                                                                                                                   >  3:37698/1‑148 (MQ=255)
                                    acgtCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGaa                                                                                                                >  8:302453/1‑149 (MQ=255)
                                     cgtCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGaa                                                                                                                <  5:207889/148‑1 (MQ=255)
                                       tCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACgg                                                                                                             <  2:110121/149‑1 (MQ=255)
                                       tCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACgg                                                                                                             <  2:104691/149‑1 (MQ=255)
                                           gTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGaagg                                                                                                         >  1:144476/1‑147 (MQ=255)
                                                tCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTcc                                                                                                    <  2:57139/149‑1 (MQ=255)
                                                tCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTcc                                                                                                    <  2:217983/149‑1 (MQ=255)
                                                  ggATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCt                                                                                                  >  1:395755/1‑149 (MQ=255)
                                                             ctgcCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAgg                                                                                       >  4:316154/1‑149 (MQ=255)
                                                              tgcCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGt                                                                                      >  1:29159/1‑149 (MQ=255)
                                                              tgcCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGt                                                                                      >  4:97161/1‑149 (MQ=255)
                                                               gcCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGtt                                                                                     >  8:57875/1‑149 (MQ=255)
                                                               gcCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGACACCGAACGGAACACCTCCCTGTCCGTGAAGGtt                                                                                     >  2:293888/1‑149 (MQ=255)
                                                                ccAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGttt                                                                                    >  5:111140/1‑149 (MQ=255)
                                                                  aGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCg                                                                                  <  7:233981/149‑1 (MQ=255)
                                                                   gATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGc                                                                                 <  7:302453/149‑1 (MQ=255)
                                                                      ctcGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCaa                                                                              >  7:204652/1‑149 (MQ=255)
                                                                        cGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCg                                                                                                                        <  7:113817/105‑1 (MQ=255)
                                                                        cGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCg                                                                                                                        >  8:113817/1‑105 (MQ=255)
                                                                        cGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTAGCCaa                                                                              >  1:344350/1‑147 (MQ=255)
                                                                         gACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATc                                                                            <  1:336444/148‑1 (MQ=255)
                                                                         gACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCg                                                                           >  2:272697/1‑149 (MQ=255)
                                                                            tttCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGGTATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGATTCGCCAATCGt                                                                          >  3:305028/1‑147 (MQ=255)
                                                                                gCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGAtt                                                                    <  4:37698/149‑1 (MQ=255)
                                                                                gCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGACACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGAt                                                                     <  1:293888/148‑1 (MQ=255)
                                                                                    ggCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCg                                                                <  7:57875/149‑1 (MQ=255)
                                                                                     gCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGa                                                               <  2:360407/149‑1 (MQ=255)
                                                                                      ctctGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAg                                                              <  8:204652/149‑1 (MQ=255)
                                                                                           ggCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCTGACGATTGGCGaaac                                                            <  2:29159/146‑4 (MQ=255)
                                                                                            gCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTAcc                                                        >  2:248548/1‑149 (MQ=255)
                                                                                                     cgAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATc                                               >  5:144806/1‑149 (MQ=255)
                                                                                                     cgAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATc                                               <  2:197520/149‑1 (MQ=255)
                                                                                                     cgAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATc                                               <  7:320328/149‑1 (MQ=255)
                                                                                                      gAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATcg                                              <  4:154297/149‑1 (MQ=255)
                                                                                                        aTGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATcgcg                                            <  2:144476/149‑1 (MQ=255)
                                                                                                        aTGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATcgcg                                            <  2:344350/149‑1 (MQ=255)
                                                                                                          gATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCt                                          >  8:99779/1‑149 (MQ=255)
                                                                                                            tCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCtgt                                        >  6:157367/1‑149 (MQ=255)
                                                                                                              ggggCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCTGACGATTGGCGaaac                                                            >  1:356762/1‑124 (MQ=255)
                                                                                                              ggggCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCTGACGATTGGCGaaac                                                            <  2:356762/127‑4 (MQ=255)
                                                                                                                  cATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAg                                  <  1:272697/149‑1 (MQ=255)
                                                                                                                  cATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAg                                  <  3:237839/149‑1 (MQ=255)
                                                                                                                   aTAACCGCGCTGGTTGTCGGCCAGCATGCCCGTGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCg                                     <  4:335968/145‑1 (MQ=255)
                                                                                                                   aTAACCGCGCTGGTTGTCGGCCAGCATGCCCGTGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCg                                     >  3:341779/1‑145 (MQ=255)
                                                                                                                   aTAACCGCGCTGGTTGTCGGCCAGCATGCCCGTGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCg                                     >  3:335968/1‑145 (MQ=255)
                                                                                                                   aTAACCGCGCTGGTTGTCGGCCAGCATGCCCGTGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCg                                     <  4:341779/145‑1 (MQ=255)
                                                                                                                      aCCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATc                                               >  7:238746/1‑132 (MQ=255)
                                                                                                                      aCCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATc                                               <  8:238746/132‑1 (MQ=255)
                                                                                                                             tGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCg                       >  2:300590/1‑149 (MQ=255)
                                                                                                                               gTTGTCGGCCAGCATGCAGCTGTCCGTGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCTTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCgg                      <  4:305028/148‑1 (MQ=255)
                                                                                                                                     ggCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCAtt               <  3:316154/149‑1 (MQ=255)
                                                                                                                                         aGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAg            >  4:151717/1‑148 (MQ=255)
                                                                                                                                         aGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGt           <  3:97161/149‑1 (MQ=255)
                                                                                                                                          gCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACCCCTCCCTGTCCGTGAAGGTTTCGCCAATCGGCAGATTGGCGAGGGTACCCGGGATATCGCGCTGTGCGTAGGTCTGTGCCCGGGGCCATTCAGTc          >  7:101131/1‑149 (MQ=255)
                                                                                                                                             tGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTtg       <  5:245191/149‑1 (MQ=255)
                                                                                                                                                ccATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCAtt               >  8:172080/1‑138 (MQ=255)
                                                                                                                                                ccATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCAtt               <  7:172080/138‑1 (MQ=255)
                                                                                                                                                  aTGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTtgccctg  <  2:395755/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TCAGGTCGTCGATGGCGTAGGGGCCGGCGCTGACGTACGTCGAGTAGATCGGATAACCATGCTGCCAGATCTCGACTTTCGCCCGGCTCTGGGCCACGCCGCGAATGATCGGGGCATAACCGCGCTGGTTGTCGGCCAGCATGCCCATGTCCGAGGCCACGCGCACGCCTGATATCGGCACCGAACGGAACACCTCCCTGTCCGTGAAGGTTTCGCCAATCGTCAGATTGGCGAGGGTACCCGGTATATCGCGCTGTGCGTAGGTCTGTGCCCGGGTCCATTCAGTCTTGCCCTG  >  NC_002947/2131229‑2131523

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