Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,296,060 C→T 27.7% intergenic (+266/+376) gltP → / ← yjcO glutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,296,0600CT27.7% 503.1 / 240.7 317intergenic (+266/+376)gltP/yjcOglutamate/aspartate:proton symporter/Sel1 family TPR‑like repeat protein
Reads supporting (aligned to +/- strand):  ref base C (120/104);  new base T (54/32);  total (179/138)
Fisher's exact test for biased strand distribution p-value = 1.61e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.66e-01

Only 100 of 319 total aligned reads displayed.

GGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCT  >  NC_000913/4295911‑4296208
                                                                                                                                                      |                                                                                                                                                    
ggTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAAc                                                                                                                                                      >  1:2849103/1‑150 (MQ=32)
ggTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAAc                                                                                                                                                      <  2:3945073/150‑1 (MQ=32)
 gTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACg                                                                                                                                                     >  1:1357007/1‑150 (MQ=32)
  tCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAATGc                                                                                                                                                    >  1:3083713/1‑150 (MQ=17)
  tCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGc                                                                                                                                                    <  2:1222716/150‑1 (MQ=32)
    ggATATGGCGTTTATGCCGCTTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCt                                                                                                                                                  <  1:2388600/150‑1 (MQ=17)
     gATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAATGCCTg                                                                                                                                                 >  1:464738/1‑150 (MQ=17)
      atttggCTTTTATGCCGCATCCGACTTCAACGCCTGTTGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTTATTTTGTTGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGa                                                                                                                                                <  1:2345152/146‑1 (MQ=11)
      aTAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAATGcc                                                                                                                                                   >  2:81932/1‑147 (MQ=17)
        aaGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAATGCCTGATg                                                                                                                                              >  1:2200401/1‑150 (MQ=17)
        aaGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATg                                                                                                                                              >  2:3849205/1‑150 (MQ=32)
         aGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTCGCGTGATTTTGTCGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGc                                                                                                                                             <  1:3804027/150‑1 (MQ=12)
          ggCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTCACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCg                                                                                                                                            <  2:3115931/150‑1 (MQ=17)
           gCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAATGCCTGATGCGa                                                                                                                                           >  2:3216124/1‑150 (MQ=17)
           gCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGa                                                                                                                                           <  2:2738282/150‑1 (MQ=32)
           gCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCATGATTTATCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAATGCCTGATGCGa                                                                                                                                           >  2:1336426/1‑150 (MQ=12)
             ttttAGGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTCGGTCGGATAAGGCTTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACg                                                                                                                                         <  1:1718142/149‑1 (MQ=12)
              tttAGGCCGCATACGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGc                                                                                                                                        <  1:3356227/150‑1 (MQ=17)
                tATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAATGCCTGATGCGACGCtt                                                                                                                                      <  2:391104/150‑1 (MQ=17)
                tATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCtt                                                                                                                                      <  2:3240236/150‑1 (MQ=32)
                 aTGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAATGCCTGATGCGACGCTTg                                                                                                                                     >  1:3490807/1‑150 (MQ=17)
                  tGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAATGCCTGATGCGACGCTTGc                                                                                                                                    >  2:1932383/1‑150 (MQ=18)
                  tGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAATGCCTGATGCGACGCTTGc                                                                                                                                    >  2:2818073/1‑150 (MQ=18)
                   gCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGCCGGATAAGGCGTTTACGCCGCATTCGGCAACCAACGCCTGATGCGACGCTTGc                                                                                                                                    <  1:1967700/150‑1 (MQ=11)
                     cGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCcgc                                                                                                                                 >  1:824463/1‑150 (MQ=35)
                        aTCCGACATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAATGCCTGATGCGACGCTTGCCGCGTc                                                                                                                              <  2:1088896/150‑1 (MQ=12)
                        aTCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTc                                                                                                                              <  2:2403333/150‑1 (MQ=37)
                         tCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAATGCCTGATGCGACGCTTGCCGCGTCt                                                                                                                             >  2:1786550/1‑150 (MQ=25)
                         tCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCt                                                                                                                             >  2:2236107/1‑150 (MQ=37)
                          ccGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCtt                                                                                                                            >  2:752528/1‑150 (MQ=37)
                           cGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACCGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGGCGGATACGGCGTTTACGCCGCATCCGAC‑AGCAACGCCTGATGCGACGCTTGCCGCGTCtaa                                                                                                                           >  2:2500990/1‑148 (MQ=255)
                            gACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTAt                                                                                                                          >  1:1617137/1‑150 (MQ=37)
                             aCATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAATGCCTGATGCGACGCTTGCCGCGTCTTATc                                                                                                                         <  1:1536780/150‑1 (MQ=25)
                               aTCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                       <  1:3850348/150‑1 (MQ=25)
                               aTCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAg                                                                                                                       <  1:752528/150‑1 (MQ=25)
                                  aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGcc                                                                                                                    >  1:1148714/1‑150 (MQ=25)
                                  aaCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTACCGCGTCTTATCAGGcc                                                                                                                    >  2:1262807/1‑150 (MQ=17)
                                   aCGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCt                                                                                                                   >  2:73778/1‑150 (MQ=25)
                                    cGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTa                                                                                                                  >  1:2138/1‑150 (MQ=25)
                                    cGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTa                                                                                                                  <  2:1751110/150‑1 (MQ=37)
                                     gCCTGATAAGACGCGTTAAGCGTCGCATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGGAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTat                                                                                                                 >  2:1644999/1‑149 (MQ=11)
                                      ccTGATGCGACGCTTTACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACg                                                                                                                <  2:544234/150‑1 (MQ=37)
                                      ccTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACg                                                                                                                >  1:2211694/1‑150 (MQ=37)
                                       cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGc                                                                                                               >  1:1085348/1‑150 (MQ=37)
                                       cTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATGTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatct                                                                                                               >  2:2374153/1‑147 (MQ=12)
                                         gATGCGACGCTTTACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCa                                                                                                             <  1:2737073/150‑1 (MQ=37)
                                         gATGCGACGCTTGACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCa                                                                                                             <  2:3878910/150‑1 (MQ=37)
                                         gATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCa                                                                                                             <  2:3473471/150‑1 (MQ=37)
                                           tGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGa                                                                                                           <  2:2210065/150‑1 (MQ=37)
                                              gACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAg                                                                                                        >  1:2659312/1‑150 (MQ=37)
                                                 gCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACCCAAGACAgcgc                                                                                                     >  2:2436202/1‑150 (MQ=255)
                                                    tAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAAt                                                                                                  <  1:2557381/150‑1 (MQ=37)
                                                      ccgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAATGCCTGATGCGACGCtt                                                                                                                                      <  1:3801028/111‑1 (MQ=14)
                                                      ccgcgTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCtt                                                                                                                                      >  2:3801028/2‑112 (MQ=14)
                                                           tCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGAGGCGACGCTTGCCGCGTCTTATCAGGCCTACCCCCGACCGCGCAAAAGCCTGa                                                                                           >  2:3720656/1‑150 (MQ=11)
                                                               aTCAGGGCGACGCCAGTCAGCGAAATCCCCTTACTTAGGCTGATTTTGTACCTCGGAGATGGCGTTTACGCCGCATCCGTC‑ATGAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTa                                                                                       <  1:561567/150‑1 (MQ=11)
                                                                tCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAATGCCTGATg                                                                                                                                              <  1:3293401/94‑1 (MQ=18)
                                                                tCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAg                                                                                      <  2:3784941/150‑1 (MQ=37)
                                                                tCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAg                                                                                      <  2:802009/150‑1 (MQ=37)
                                                                 cAGGCCTACGCCAGACAGCGCATTAGCCTGATTTTGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGc                                                                                     <  2:1571851/150‑1 (MQ=255)
                                                                     ccTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTatc                                                                                                                >  2:2973869/1‑117 (MQ=21)
                                                                     ccTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAAAAGACTGATTAAGCGTGa                                                                                 >  1:3416455/1‑150 (MQ=255)
                                                                      cTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGAt                                                                                >  2:2124905/1‑150 (MQ=17)
                                                                          gCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTg                                                                            >  1:1833607/1‑150 (MQ=37)
                                                                            cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCt                                                                                                                                                  >  1:3177847/1‑78 (MQ=34)
                                                                            cAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTa                                                                          <  2:709235/150‑1 (MQ=37)
                                                                                       aTAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAg                                                               <  2:3013937/150‑1 (MQ=37)
                                                                                        tAGCCTGATTTTTCGTGAATTTTTATCTAGTTTAAGGCTTTTACGCCGCATCATTA‑TTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCGATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAgg                                                              <  2:2721559/150‑1 (MQ=11)
                                                                                        tAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAgg                                                              >  2:2325001/1‑150 (MQ=37)
                                                                                         aGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGc                                                             <  2:1833605/150‑1 (MQ=37)
                                                                                          gCCTGATATAGCGTGATTTGGTAGGTCGGCTATGGGTTTTACGCCGCCTACTCC‑TTGAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCg                                                            <  2:3581619/150‑1 (MQ=11)
                                                                                           ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAATGCCTGAt                                                                                                                                               <  1:3929865/66‑1 (MQ=25)
                                                                                           ccTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGt                                                           >  1:1017775/1‑150 (MQ=37)
                                                                                            cTGATTTAGCTTGATTTTGTTGTTATGTTACGGCTTTTACGCCGGTTACGCT‑TTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGtt                                                          <  2:2853624/150‑1 (MQ=11)
                                                                                               aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTAt                                                       <  1:2256443/150‑1 (MQ=37)
                                                                                                 ttAGCTTGATTTTGTAGGTCGGATATGGCTTTTACTCCGCATACGTC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGc                                                     <  1:3369855/150‑1 (MQ=11)
                                                                                                       tGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGGGAGACAGCGCTATAGGCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATc                                               >  2:3577330/1‑150 (MQ=255)
                                                                                                       tGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAATTCGTTTAAGCCGCATc                                               >  2:2672022/1‑150 (MQ=255)
                                                                                                         aTTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCg                                             >  1:3430667/1‑150 (MQ=37)
                                                                                                                 ggTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAAc                                     <  1:1261649/150‑1 (MQ=37)
                                                                                                                  gTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGt                                                                       <  2:2608996/115‑1 (MQ=39)
                                                                                                                   tCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTc                                                                                                                              >  1:269002/1‑59 (MQ=255)
                                                                                                                   tCGGATAAGGCGTTTAAGCCGCATCCGAC‑ATCAAAGCCTGATGCGACGCTTGCCGCGt                                                                                                                               >  1:3006535/1‑58 (MQ=14)
                                                                                                                     ggATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGACATCAACGCCt                                 >  1:3152153/1‑150 (MQ=21)
                                                                                                                      gATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTg                                >  1:1791092/1‑150 (MQ=37)
                                                                                                                       aTAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGa                               <  2:2839987/150‑1 (MQ=37)
                                                                                                                       aTAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGa                               >  2:3173878/1‑150 (MQ=37)
                                                                                                                            gCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGt                                                           >  2:3850021/1‑117 (MQ=39)
                                                                                                                                   cgccgcCTCCTGC‑GTCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTaa                   <  2:833298/150‑1 (MQ=255)
                                                                                                                                   cgccgcATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTaa                   >  2:2473448/1‑150 (MQ=35)
                                                                                                                                     ccgcATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAcg                 <  2:2659312/150‑1 (MQ=34)
                                                                                                                                      cgcATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAAcgc                >  2:3419016/1‑150 (MQ=34)
                                                                                                                                             gAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCtt                     <  1:3658811/138‑1 (MQ=34)
                                                                                                                                             gAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCtt                     >  2:3658811/1‑138 (MQ=34)
                                                                                                                                             gAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTAt         >  2:423375/1‑150 (MQ=34)
                                                                                                                                             gAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTAt         >  1:3790384/1‑150 (MQ=34)
                                                                                                                                             gAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTAt         >  1:292162/1‑150 (MQ=34)
                                                                                                                                                  tCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAgg     >  2:32848/1‑150 (MQ=33)
                                                                                                                                                    aaCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGcc   >  2:1536780/1‑150 (MQ=32)
                                                                                                                                                    aaCGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGcc   >  2:1454730/1‑150 (MQ=32)
                                                                                                                                                      |                                                                                                                                                    
GGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTACGCCGCATCCGAC‑ATCAACGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTACGCCAGACAGCGCAATAGCCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTATGCCGCATCCGACATCAACGCCTGATGCGACGCTTAACGCGTCTTATCAGGCCT  >  NC_000913/4295911‑4296208

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