Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,318,613 A→G 5.4% Y153H (TAC→CAC) ‡ ilvC ← ketol‑acid reductoisomerase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,318,6130AG 5.4% 1002.8 / 27.4 311Y153H (TAC→CAC) ‡ilvCketol‑acid reductoisomerase
Reads supporting (aligned to +/- strand):  ref base A (141/153);  new base G (11/6);  total (152/159)
Fisher's exact test for biased strand distribution p-value = 2.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

Only 100 of 316 total aligned reads displayed.

GGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAGGTCGGCACGCGGAACAACCTGGTTGTAGTGGATCGAGAAGCCGTGAGAGAAGGCCA  >  NC_002947/5318464‑5318762
                                                                                                                                                      |                                                                                                                                                     
   aCCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGAt                                                                                                                                                    >  2:8246953/1‑150 (MQ=255)
     cGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAg                                                                                                                                                  >  2:7500161/1‑150 (MQ=255)
      gCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGc                                                                                                                                                 >  2:8067612/1‑150 (MQ=255)
         gcaaaGAACAGCCTGCTCACCGAACAGGTCAGTTTCGCACTCTTCCTTGAAGGTGGTTTCAATAACGCCCGAACGACCGCCACCAACGCCCGAAGCGTAAGACAGAGCGAGATCATGTGCGTTGCCCGAAGCATCCTGATGGATCGCGAt                                                                                                                                              >  2:2944029/5‑150 (MQ=255)
          cacaGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGAtc                                                                                                                                             >  1:1351904/1‑150 (MQ=255)
          cacaGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGAtc                                                                                                                                             >  2:1904532/1‑150 (MQ=255)
           acaGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGAtca                                                                                                                                            >  1:1853339/1‑150 (MQ=255)
             aGAACAGCCTGCTCACCGAACAGGTCAGTTTCGCACTCTTCCTTGAAGGTGGTTTCAATAACGCCCGAACGACCGCCACCAACGCCCGAAGCGTAAGACAGAGCGAGATCATGTGCGTTGCCCGAAGCATCCTGATGGATCGCGAtaagg                                                                                                                                          >  1:1124571/1‑150 (MQ=255)
                aCAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCTATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGAtcaggtca                                                                                                                                       <  2:5331212/150‑1 (MQ=255)
                 cAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGAtcaggtcag                                                                                                                                      <  2:4553486/150‑1 (MQ=255)
                 cAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGAtcaggtcag                                                                                                                                      <  1:8421947/150‑1 (MQ=255)
                 cAGCCTGCTCACCGAACAGGTCAGTTTCGCACTCTTCCTTGAAGGTGGTTTCAATAACGCCCTAACGACCGCCACCAACGCCCGAAGCGTAAGACAGAGCGAGATCATGTGCGTTGCCCGAAGCATCCTGATGGATCGCGAtaaggcaag                                                                                                                                      <  2:1551534/150‑5 (MQ=255)
                  aGCCTGCTCACCGAACAGGTCAGTTTCGCACTCTTCCTTGAAGGTGGTTTCAATAACGCCCGAACGACCGCCACCAACGCCCGAAGCGTAAGACAGAGCGAGATCATGTGCGTTGCCCGAAGCATCCTGATGGATCGCGAtaaggcaagg                                                                                                                                     >  2:2235839/1‑145 (MQ=255)
                      tGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATg                                                                                                                                 <  1:982107/150‑1 (MQ=255)
                        cTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATgcc                                                                                                                               <  2:7452911/150‑1 (MQ=255)
                          cACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATgccgc                                                                                                                             >  1:7259570/1‑150 (MQ=255)
                           aCCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGACAGGGATgccgcc                                                                                                                            >  2:1193534/1‑150 (MQ=255)
                            ccGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATgccgccg                                                                                                                           >  2:1789614/1‑150 (MQ=255)
                               aaCAGGTCGGGTGCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCt                                                                                                                        <  2:9105279/150‑1 (MQ=255)
                                aCAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCtt                                                                                                                       <  2:3460187/150‑1 (MQ=255)
                                   ggTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGa                                                                                                                    <  2:4284323/150‑1 (MQ=255)
                                   ggTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGa                                                                                                                    <  1:2047394/150‑1 (MQ=255)
                                    gTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGAc                                                                                                                   <  2:2458588/150‑1 (MQ=255)
                                    gTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGAc                                                                                                                   >  2:1615798/1‑150 (MQ=255)
                                    gTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGAc                                                                                                                   <  2:8584347/150‑1 (MQ=255)
                                          ttcggtttcgTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTc                                                                                                             >  2:2067627/1‑150 (MQ=255)
                                             ggtttcgTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGa                                                                                                          <  1:4370708/150‑1 (MQ=255)
                                             ggtttcgTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGa                                                                                                          <  1:643372/150‑1 (MQ=255)
                                                   gTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGa                                                                                                                                   >  2:2427050/1‑119 (MQ=255)
                                                    tCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGGGCGAACg                                                                                                   <  2:4674151/150‑1 (MQ=255)
                                                    tCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACg                                                                                                   >  1:4068857/1‑150 (MQ=255)
                                                      ctggAAGGTGGTTTCGATGTTACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGTAGCGTCCTGGTAGTTGGCGATCAGGTCAGGGATGCCGCCGCCTGTGACGAACTCGGAGCGAACGgt                                                                                                 >  2:3004768/4‑150 (MQ=255)
                                                      cTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGTCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGgt                                                                                                 <  1:3757459/150‑1 (MQ=255)
                                                      cTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGgt                                                                                                 >  1:5852149/1‑150 (MQ=255)
                                                      cTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGgt                                                                                                 <  2:3404231/150‑1 (MQ=255)
                                                       ttGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGgtg                                                                                                >  2:2688828/1‑150 (MQ=255)
                                                        tGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGgtgt                                                                                               <  1:4543161/150‑1 (MQ=255)
                                                        tGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGgtgt                                                                                               <  1:8819239/150‑1 (MQ=255)
                                                         gAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGtgtg                                                                                              <  2:446843/150‑1 (MQ=255)
                                                         gAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGtgtg                                                                                              <  1:639131/150‑1 (MQ=255)
                                                            ggtggtTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGcc                                                                                           >  2:6463989/1‑150 (MQ=255)
                                                            ggtggtTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGcc                                                                                           <  1:622947/150‑1 (MQ=255)
                                                              tggtTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCg                                                                                         <  2:4831854/150‑1 (MQ=255)
                                                              tggtTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCg                                                                                         <  2:5881945/150‑1 (MQ=255)
                                                               ggtTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCgg                                                                                        >  1:743982/1‑150 (MQ=255)
                                                                gtTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCggg                                                                                       <  1:1222721/150‑1 (MQ=255)
                                                                     gatgatACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACTGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCttt                                                                                  <  2:3822611/150‑1 (MQ=255)
                                                                     gatgatACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCttt                                                                                  >  2:7361320/1‑150 (MQ=255)
                                                                     gatgatACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCttt                                                                                  >  1:2506480/1‑150 (MQ=255)
                                                                      atgatACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTc                                                                                 <  2:5629473/150‑1 (MQ=255)
                                                                       tgatACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCg                                                                                <  1:1003549/150‑1 (MQ=255)
                                                                       tgatACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCg                                                                                <  1:300555/150‑1 (MQ=255)
                                                                       tgatACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCg                                                                                >  2:576568/1‑150 (MQ=255)
                                                                           aCCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGgcgc                                                                            <  2:8814713/150‑1 (MQ=255)
                                                                           aCCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGgcgc                                                                            >  1:7122408/1‑150 (MQ=255)
                                                                            ccGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGcgcg                                                                           <  1:341935/150‑1 (MQ=255)
                                                                            ccGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGcgcg                                                                           >  2:710402/1‑150 (MQ=255)
                                                                              ggTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGAt                                                                         >  1:2296338/1‑150 (MQ=255)
                                                                               gTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATc                                                                        <  1:6578671/150‑1 (MQ=255)
                                                                               gTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATc                                                                        >  2:4059559/1‑150 (MQ=255)
                                                                               gTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATc                                                                        >  2:7588210/1‑150 (MQ=255)
                                                                                 aCGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCAt                                                                      <  1:4226541/150‑1 (MQ=255)
                                                                                 aCGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCAt                                                                      <  2:6824613/150‑1 (MQ=255)
                                                                                  cGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATg                                                                     <  2:2897325/150‑1 (MQ=255)
                                                                                    gCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGAt                                                                   <  2:5471257/150‑1 (MQ=255)
                                                                                     ccaccaCCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATc                                                                  >  2:2878864/1‑150 (MQ=255)
                                                                                     ccaccaCCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATc                                                                  <  1:3192517/150‑1 (MQ=255)
                                                                                          accTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTc                                                             >  1:5102773/1‑150 (MQ=255)
                                                                                            cTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGa                                                           <  2:7885992/150‑1 (MQ=255)
                                                                                            cTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGa                                                           >  2:5957898/1‑150 (MQ=255)
                                                                                             tACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAg                                                          >  1:2177303/1‑150 (MQ=255)
                                                                                             tACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAg                                                          >  1:3947726/1‑150 (MQ=255)
                                                                                                  ccGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAGGTCgg                                                     >  1:6142356/1‑150 (MQ=255)
                                                                                                  ccGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAGGTCgg                                                     <  1:8246953/150‑1 (MQ=255)
                                                                                                   cGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAGGTCgac                                                    <  2:6432954/150‑3 (MQ=255)
                                                                                                         aTAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAGGTCGGCACGCgg                                              >  2:6216731/1‑150 (MQ=255)
                                                                                                          tAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAGGTCGGCACGCGGa                                             >  1:6604191/1‑150 (MQ=255)
                                                                                                          ggcatgaGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAGGTCGGCACGCGGa                                             <  1:7004714/144‑1 (MQ=255)
                                                                                                                  tGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAGGTCGGCACGCGGAACAACCTg                                     <  1:7705771/150‑1 (MQ=255)
                                                                                                                            ttGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAGGTCGGCACGCGGAACAACCAGGTTGttccgc                           <  1:2878864/150‑6 (MQ=255)
                                                                                                                                  gTTGCCGGAAGTGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAGGTCGGCACGCGGAACAACCTGGTTGTAGTGGATCga                      >  2:8484136/1‑150 (MQ=255)
                                                                                                                                     gCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAGGTCGGCACGCGGAACAACCTGGTTGTAGTGGATCGAGaa                   <  1:1407395/150‑1 (MQ=255)
                                                                                                                                      ccGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAGGTCGGCACGCGGAACAACCTGGTTGTAGTGGATCGAGAAg                  >  1:900403/1‑150 (MQ=255)
                                                                                                                                      ccGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAGGTCGGCACGCGGAACAACCTGGTTGTAGTGGATCGAGAAg                  >  1:5809604/1‑150 (MQ=255)
                                                                                                                                        ggAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAGGTCGGCACGCGGAACAACCTGGTTGTAGTGGATCGAGAAGcc                <  1:7511609/150‑1 (MQ=255)
                                                                                                                                          aaGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAGGTCGGCACGCGGAACAACCTGGTTGTAGTGGATCGAGAAGCCGt              <  1:2243868/150‑1 (MQ=255)
                                                                                                                                          aaGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAGGTCGGCACGCGGAACAACCTGGTTGTAGTGGATCGAGAAGCCGt              <  1:1615798/150‑1 (MQ=255)
                                                                                                                                          aaGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAGGTCGGCACGCGGAACAACCTGGTTGTAGTGGATCGAGAAGCCGt              <  2:5265672/150‑1 (MQ=255)
                                                                                                                                            gCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAGGTCGGCACGCGGAACAACCTGGTTGTAGTGGATCGAGAAGCCGTga            <  1:1564154/150‑1 (MQ=255)
                                                                                                                                            gCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAGGTCGGCACGCGGAACAACCTGGTTGTAGTGGATCGAGAAGCCGTga            <  2:3669048/150‑1 (MQ=255)
                                                                                                                                            gCGTCCTGGTAGAGAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAGGTCGGCACGCGGAACAACCTGGTTGTAGTGGATCGAGAAGCCGTga            <  2:6668333/150‑1 (MQ=255)
                                                                                                                                              gTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAGGTCGGCACGCGGAACAACCTGGTTGTAGTGGATCGAGAAGCCGTgaga          <  1:3437067/150‑1 (MQ=255)
                                                                                                                                              gTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAGGTCGGCACGCGGAACAACCTGGTTGTAGTGGATCGAGAAGCCGTgaga          >  1:8873214/1‑150 (MQ=255)
                                                                                                                                              gTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAGGTCGGCACGCGGAACAACCTGGTTGTAGTGGATCGAGAAGCCGTgaga          >  1:6051039/1‑150 (MQ=255)
                                                                                                                                                 cTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAGGTCGGCACGCGGAACAACCTGGTTGTAGTGGATCGAGAAGCCGTGAGAGaa       <  2:1351904/150‑1 (MQ=255)
                                                                                                                                                  tGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAGGTCGGCACGCGGAACAACCTGGTTGTAGTGGATCGAGAAGCCGTGAGAGAAg      <  2:5546775/150‑1 (MQ=255)
                                                                                                                                                  tGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAGGTCGGCACGCGGAACAACCTGGTTGTAGTGGATCGAGAAGCCGTGAGAGAAg      >  1:6019964/1‑150 (MQ=255)
                                                                                                                                                   ggTAGATAGCGATCAGGACAGGGATGCCGCCGCCTGTGACGAACTCGGAGCGAACGGTGGGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAGGTCGGCACGCGGAACAACCGGGTAGTAGGGGAGCGAGACTCCGTgagagacgc     >  2:7878108/1‑147 (MQ=255)
                                                                                                                                                    gTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAGGTCGGCACGCGGAACAACCTGGTTGTAGTGGATCGAGAAGCCGTGAGAGAAGGc    >  2:3790334/1‑150 (MQ=255)
                                                                                                                                                      aGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAGGTCGGCACGCGGAACAACCTGGTTGTAGTGGATCGAGAAGCCGTGAGAGAAGGCCa  >  1:5257851/1‑150 (MQ=255)
                                                                                                                                                      |                                                                                                                                                     
GGTACCGCCACACAGAACAGCCTGCTCACCGAACAGGTCGGTTTCGGTTTCGTCCTTGAAGGTGGTTTCGATGATACCGGTACGGCCACCACCTACGCCCGAAGCATAGGACAGTGCGACGTTCTTGGCGTTGCCGGAAGCGTCCTGGTAGATAGCGATCAGGTCAGGGATGCCGCCGCCTTTGACGAACTCGGAGCGAACGGTGTGGCCCGGGGCTTTCGGCGCGATCATGATCACGTCGAGGTCGGCACGCGGAACAACCTGGTTGTAGTGGATCGAGAAGCCGTGAGAGAAGGCCA  >  NC_002947/5318464‑5318762

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