Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,296,060 C→T 39.0% intergenic (+266/+376) gltP → / ← yjcO glutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,296,0600CT39.0% 137.0 / 214.1 201intergenic (+266/+376)gltP/yjcOglutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base C (56/62);  new base T (53/24);  total (115/86)
Fisher's exact test for biased strand distribution p-value = 4.90e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-02

Only 100 of 201 total aligned reads displayed.

CGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCT  >  NC_000913/4295914‑4296208
                                                                                                                                                  |                                                                                                                                                    
 ggATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCt                                                                                                                                                  >  2:1672828/1‑150 (MQ=17)
  taTACGGCGTTTATGCCGCTTAAGAGTTCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTg                                                                                                                                                 <  2:1376020/149‑1 (MQ=11)
  gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCTACTCTTTACGCTTCTTATCAGGCCTACGCCATACAGCTCAATAGACTGATTTAGCTTTATTTTGTATGTCGTATAAGGCGTTTACGCCGCATCCGACATCAATGCCTg                                                                                                                                                 <  2:652202/150‑1 (MQ=11)
  gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTg                                                                                                                                                 >  1:1956523/1‑150 (MQ=32)
  gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCGGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGGTTACGCCGCATCCGACATCAACGCCTg                                                                                                                                                 >  1:1289829/1‑150 (MQ=17)
   aTAAGGCTTTTATTCCGCATCCTACTTCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGa                                                                                                                                                <  2:1240246/150‑1 (MQ=255)
   aTAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGa                                                                                                                                                >  2:219087/1‑150 (MQ=32)
     aaGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATg                                                                                                                                              <  1:1554393/150‑1 (MQ=32)
       ggCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCg                                                                                                                                            >  1:3072/1‑150 (MQ=32)
        gCGTTTATGCCGCTTAAGTCATCAACGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGa                                                                                                                                           <  2:1455481/150‑1 (MQ=11)
              aTGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGAGACGCTTg                                                                                                                                     >  1:1085796/1‑150 (MQ=17)
               tGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGc                                                                                                                                    >  1:903746/1‑150 (MQ=18)
               tGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGc                                                                                                                                    >  1:1622501/1‑150 (MQ=18)
                gCCTCATACTACATCAAAGCCTGATGCGACGCTTTACGCGTCTTATCAGGCCTACGCCAGACAGATCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGcc                                                                                                                                   <  2:354746/150‑1 (MQ=11)
                gCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGcc                                                                                                                                   <  1:715707/150‑1 (MQ=21)
                gCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGcc                                                                                                                                   >  2:454150/1‑150 (MQ=21)
                gCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGcc                                                                                                                                   >  1:254525/1‑150 (MQ=21)
                gCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGcc                                                                                                                                   >  1:1702958/1‑150 (MQ=21)
                gCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGcc                                                                                                                                   >  1:1642378/1‑150 (MQ=21)
                gCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGcc                                                                                                                                   >  2:229098/1‑150 (MQ=21)
                gCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGc                                                                                                                                    <  2:1556769/149‑1 (MQ=33)
                 ccGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcg                                                                                                                                  >  1:1200174/1‑150 (MQ=25)
                  cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcgc                                                                                                                                 >  2:399102/1‑150 (MQ=25)
                   gCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCcgcg                                                                                                                                >  2:1834132/1‑150 (MQ=25)
                     aTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTc                                                                                                                              <  1:1407897/150‑1 (MQ=37)
                      tCCGACATCAACGCCTTATGCTACGCTTATCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCt                                                                                                                             <  1:880524/150‑1 (MQ=255)
                      tCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCt                                                                                                                             >  1:1448242/1‑150 (MQ=25)
                       ccGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCtt                                                                                                                            >  1:1168172/1‑150 (MQ=25)
                       ccGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCtt                                                                                                                            >  2:866216/1‑150 (MQ=25)
                        cGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                           >  1:1273172/1‑150 (MQ=25)
                        cGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTa                                                                                                                           >  2:1064211/1‑150 (MQ=25)
                         gACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                          >  1:1029136/1‑150 (MQ=25)
                         gACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAAGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                          >  1:825949/1‑150 (MQ=37)
                            aTCAATGCCTTATTCGACTCTTTACTCGTCTTATCAGGCCTACTCCTTACAGATCAATATCATGATTTAGCTTGATTTTGTATTTAGGATATGGCTTTTACTCCGCTTAAGACTTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                       <  2:1956724/150‑1 (MQ=2)
                            aTCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                       >  2:779721/1‑150 (MQ=25)
                             tCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAgg                                                                                                                      >  2:1028404/1‑150 (MQ=25)
                              cAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGGTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                     >  1:1724399/1‑150 (MQ=25)
                               aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                    >  1:19791/1‑150 (MQ=25)
                               aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                    >  1:47381/1‑150 (MQ=25)
                               aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                    >  1:54707/1‑150 (MQ=25)
                               aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                    <  2:1489885/150‑1 (MQ=37)
                               aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                    <  2:1574558/150‑1 (MQ=37)
                               aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                    >  1:746641/1‑150 (MQ=37)
                               aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                    <  1:1745037/150‑1 (MQ=37)
                                aCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCt                                                                                                                   >  1:926273‑M2/2‑150 (MQ=255)
                                 ccccTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATACCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTa                                                                                                                  <  1:895157/148‑1 (MQ=255)
                                   ccTGATGCGACGCTTGACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACg                                                                                                                                                     <  2:1725758/113‑1 (MQ=17)
                                   ccTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                                                                                                                <  1:1462663/150‑3 (MQ=18)
                                   ccTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                                                                                                                >  1:79559/1‑148 (MQ=18)
                                     tGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGcc                                                                                                              <  2:626265/150‑1 (MQ=37)
                                         gCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                         <  1:218300/135‑1 (MQ=25)
                                         gCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGAc                                                                                                          <  1:1668190/150‑1 (MQ=37)
                                          cGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACa                                                                                                         >  1:1561445/1‑150 (MQ=37)
                                             cGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgcg                                                                                                      <  1:1042115/150‑1 (MQ=37)
                                                  gccgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAAGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAgcgcgagg                                                                                                 >  1:1249912/3‑146 (MQ=12)
                                                  aaCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAAGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAAGGGAGACAGCGCATTa                                                                                                 >  2:1435837/1‑150 (MQ=11)
                                                  aaCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGAGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAAGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAAGAGAGCCAGCGCAATa                                                                                                 >  2:1428810/1‑150 (MQ=11)
                                                  aaCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCATGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAAGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAAGAAAGACAGCGCAATa                                                                                                 >  2:266872/1‑150 (MQ=11)
                                                     gcttattATGAGGCCTACGCGATTGTGAGCATTAGTCTCATTTAGGTTGATTTTGTTTTTCTGTTAAGGCTTTTACTCCGGATACGACTTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGGCTACGCCAGACAGCGCAATAGcc                                                                                              <  2:29952/145‑1 (MQ=2)
                                                     gcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGc                                                                                                                                                    >  1:1455720/1‑96 (MQ=33)
                                                       gTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTg                                                                                            >  2:725583/1‑150 (MQ=37)
                                                        tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCt                                                                                                                   <  2:673277/126‑1 (MQ=25)
                                                        tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCt                                                                                                                   >  1:673277/1‑126 (MQ=25)
                                                         cTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCt                                                                                                                             >  1:434798/1‑115 (MQ=17)
                                                                   cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                    >  1:1834535/1‑114 (MQ=37)
                                                                       gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCt                                                                                                                                       <  1:1939503/91‑1 (MQ=18)
                                                                       gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                                                                                                                <  1:527289/114‑3 (MQ=21)
                                                                       gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAAAGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                                                                                                                >  2:527289/1‑112 (MQ=255)
                                                                            aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTg                                                                                                                                                 >  1:1583303/1‑76 (MQ=21)
                                                                            aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTg                                                                                                                                     >  1:1003349/1‑88 (MQ=25)
                                                                            aCAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGCTTg                                                                                                                                     <  2:1003349/88‑1 (MQ=25)
                                                                               gcgcAATAGCCTGATTTAGCGTGATTTTGTAGGTCAGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                          >  1:621564/1‑96 (MQ=255)
                                                                                 gcAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATg                                                                                                                                              <  2:1886069/74‑1 (MQ=34)
                                                                                   aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGAc                                                                                                                                          <  2:1015908/76‑1 (MQ=34)
                                                                                   aaTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCcg                                                                                                                                  <  2:1219745/84‑1 (MQ=39)
                                                                                    aTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGc                                                                                                                     >  2:1627685/1‑96 (MQ=40)
                                                                                       gCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAATGCCTGATGCGACGc                                                                                                                                        <  2:1080065/74‑1 (MQ=21)
                                                                                             tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAg                                                                         >  2:1800742/1‑131 (MQ=37)
                                                                                                     gATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCg                                                                                    >  2:413389/1‑112 (MQ=39)
                                                                                                     gATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAAGCCGCATcc                                              >  2:1834636/1‑150 (MQ=37)
                                                                                                        tttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGAc                                           >  1:495051/1‑150 (MQ=37)
                                                                                                         ttGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACTACTGATTCGTCGCTTTCCGCGTCTTATCAGTCCTACGCCAGACAGCGCAATAGCCTGATTTAGCTTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACa                                          <  2:606272/150‑1 (MQ=11)
                                                                                                          tGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACAt                                         >  2:748837/1‑150 (MQ=21)
                                                                                                            tAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCaa                                      >  1:383149/1‑150 (MQ=37)
                                                                                                            gatgTCGGATGCGGCGTAAACGCCTTATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCa                                       >  2:1662261/4‑150 (MQ=11)
                                                                                                              ggTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGc                                                                                     <  1:768313/102‑1 (MQ=17)
                                                                                                              ggTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGc                                                                                                                                             >  1:1229003/1‑46 (MQ=38)
                                                                                                                    aTAAGGCGTTTACGCCGCATCCGACATCAACGCCTGTTTTTATTATTGACTCGTCTTATCAGGCCTACTACAGTCATCGAAATTGTTTGATTTATCGTGTTTTTTTAGTTCTGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGa                               <  2:1289829/150‑1 (MQ=2)
                                                                                                                       aGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGc                            >  1:652202/1‑150 (MQ=37)
                                                                                                                       aGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAGGGCGTTTATGCCGCATCCGACATCAACGCCTGATgg                            <  1:659749/150‑2 (MQ=17)
                                                                                                                       aGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAGGGCGTTTATGCCGCATCCGACATCAACGCCTGATGc                            >  2:659749/1‑150 (MQ=17)
                                                                                                                           gTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGATCGGATAAGGACTATAAGCAGCAACCGACAACAACGCCTGATGCGACg                        >  1:1021960/1‑150 (MQ=11)
                                                                                                                                      aTCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTc             <  1:1860866/150‑1 (MQ=34)
                                                                                                                                         cGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGGCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATAAGACGCGTTAAGCGTcgca          <  1:398283/150‑4 (MQ=11)
                                                                                                                                           aCATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATc        <  2:1321790/150‑1 (MQ=34)
                                                                                                                                              tCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAgg     >  1:1537021/1‑150 (MQ=33)
                                                                                                                                               cAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAAGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc    >  2:715707/1‑150 (MQ=17)
                                                                                                                                               cAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGACAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGc    >  2:1929771/1‑150 (MQ=17)
                                                                                                                                               cAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTAAGCAAGACAGCGCAAGAGCCTGATTTAGCGTGATATTTTAGGTAGGATAAGGAGTTTAAGCAGAATAAGACATAAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGc    >  2:895157/1‑150 (MQ=11)
                                                                                                                                                 aCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCt  <  1:611333/150‑1 (MQ=31)
                                                                                                                                                  |                                                                                                                                                    
CGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCT  >  NC_000913/4295914‑4296208

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