Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,296,380 1 bp→GCG intergenic (+586/+56) 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,3791.G91.9% 232.2 / 33.4 123intergenic (+585/+57)gltP/yjcOglutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base . (9/1);  new base G (42/71);  total (51/72)
Fisher's exact test for biased strand distribution p-value = 1.55e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.16e-01
*NC_0009134,296,3792.C91.9% 232.4 / 34.0 123intergenic (+585/+57)gltP/yjcOglutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base . (9/1);  new base C (42/71);  total (51/72)
Fisher's exact test for biased strand distribution p-value = 1.55e-03
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.01e-01
*NC_0009134,296,3800AG99.2% 254.9 / NA 123intergenic (+586/+56)gltP/yjcOglutamate/aspartate : H(+) symporter GltP/Sel1 repeat‑containing protein YjcO
Reads supporting (aligned to +/- strand):  ref base A (0/0);  major base G (50/72);  minor base C (1/0);  total (51/72)
Fisher's exact test for biased strand distribution p-value = 4.15e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 149 total aligned reads displayed.

CCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACC  >  NC_000913/4296341‑4296482
                                         |||||                                                                                                    
 cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTaa                                             <  2:7326211/100‑1 (MQ=11)
 cTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAc                                            <  2:7326173/101‑1 (MQ=1)
  tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                                           <  2:7326206/101‑1 (MQ=11)
  tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                                           >  1:7326504/1‑101 (MQ=11)
  tGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAAcc                                           >  1:7326503/1‑101 (MQ=11)
   gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                                          <  1:7325624/101‑1 (MQ=11)
   gATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCg                                          >  2:7326506/1‑101 (MQ=11)
    aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGt                                         <  2:7325861/101‑1 (MQ=11)
    aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGt                                         <  1:7326057/101‑1 (MQ=11)
    aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGt                                         <  1:7326180/101‑1 (MQ=11)
    aTTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTTGCC‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGt                                         >  2:7326507/1‑101 (MQ=11)
     tttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGtt                                        <  2:7326170/101‑1 (MQ=11)
      ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                                       <  2:7326185/101‑1 (MQ=11)
      ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                                       <  2:7326165/101‑1 (MQ=11)
      ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                                       >  1:7326508/1‑101 (MQ=11)
      ttAGCGTGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTg                                       <  2:7325390/101‑1 (MQ=11)
       tAGCGTGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                                      <  1:7325821/101‑1 (MQ=11)
       tAGCGTGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTgg                                      >  1:7326510/1‑101 (MQ=11)
        aGCGTGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                                     <  1:7326144/101‑1 (MQ=11)
        aGCGTGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGt                                     >  1:7326512/1‑101 (MQ=11)
         gCGTGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGtt                                    >  1:7326513/1‑101 (MQ=14)
         gCGTGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGtt                                    <  1:7325866/101‑1 (MQ=14)
          cGTGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTa                                   >  2:7326514/1‑101 (MQ=14)
          cGTGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTa                                   <  2:7325854/101‑1 (MQ=14)
           gTGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTaa                                  >  1:7326517/1‑101 (MQ=14)
           gTGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTaa                                  <  2:7326186/101‑1 (MQ=14)
           gTGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTaa                                  >  1:7326518/1‑101 (MQ=14)
           gTGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTaa                                  <  1:7325848/101‑1 (MQ=14)
           gTGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTaa                                  <  2:7326124/101‑1 (MQ=14)
            tGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt                                 <  2:7326172/101‑1 (MQ=14)
            tGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt                                 <  1:7325550/101‑1 (MQ=14)
            tGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt                                 >  2:7326520/1‑101 (MQ=14)
            tGATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt                                 >  2:7326523/1‑101 (MQ=14)
            tGATTTTGTAGGTCGGATAAGGCGTTT‑‑GCGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAt                                 >  1:7326522/1‑101 (MQ=14)
             gATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt                                >  2:7326526/1‑101 (MQ=14)
             gATTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtt                                <  2:7326089/101‑1 (MQ=14)
              aTTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttt                               <  2:7325641/101‑1 (MQ=14)
              aTTTTGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttt                               <  2:7326013/101‑1 (MQ=14)
               ttttGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAtttt                              >  2:7326527/1‑101 (MQ=14)
                tttGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                             >  1:7326530/1‑101 (MQ=14)
                tttGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                             <  1:7326065/101‑1 (MQ=14)
                tttGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                             >  2:7326529/1‑101 (MQ=14)
                tttGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAAttttt                             >  2:7326531/1‑101 (MQ=14)
                 ttGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTc                            >  1:7326533/1‑101 (MQ=17)
                  tGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa                           <  2:7325838/101‑1 (MQ=17)
                  tGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa                           <  2:7325647/101‑1 (MQ=17)
                  tGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa                           >  2:7326535/1‑101 (MQ=17)
                  tGTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCa                           >  1:7326536/1‑101 (MQ=17)
                   gTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa                          <  2:7325885/101‑1 (MQ=17)
                   gTAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaa                          <  2:7326139/101‑1 (MQ=17)
                    tAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                         >  2:7326537/1‑101 (MQ=17)
                    tAGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCaaa                         <  2:7325828/101‑1 (MQ=17)
                     aGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc                        <  1:7325869/101‑1 (MQ=17)
                     aGGTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAc                        <  1:7325802/101‑1 (MQ=17)
                      ggTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAAct                       >  2:7326539/1‑101 (MQ=17)
                       gTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc                      <  1:7326153/101‑1 (MQ=17)
                       gTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc                      <  2:7326200/101‑1 (MQ=17)
                       gTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc                      >  2:7326540/1‑101 (MQ=18)
                       gTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc                      <  2:7326079/101‑1 (MQ=17)
                       gTCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActc                      <  2:7325897/101‑1 (MQ=17)
                        tCGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAActct                     <  2:7325474/101‑1 (MQ=18)
                         cGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                    >  2:7326544/1‑101 (MQ=18)
                         cGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                    >  2:7326543/1‑101 (MQ=18)
                         cGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                    <  2:7326001/101‑1 (MQ=18)
                         cGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                    <  1:7326111/101‑1 (MQ=18)
                         cGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                    <  1:7325818/101‑1 (MQ=18)
                         cGGATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCtt                    <  2:7326198/101‑1 (MQ=18)
                          ggATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc                   <  1:7326094/101‑1 (MQ=18)
                           gATAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                  <  2:7326513/101‑1 (MQ=18)
                           gATAAGGCGTTT‑‑GCGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTca                  >  1:7326546/1‑101 (MQ=18)
                            aTAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac                 <  1:7326181/101‑1 (MQ=18)
                            aTAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac                 <  1:7326112/101‑1 (MQ=18)
                            aTAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac                 <  2:7326207/101‑1 (MQ=18)
                            aTAAGGCGTTTGCG‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac                 <  2:7325591/101‑1 (MQ=18)
                            aTAAGGCGTTT‑‑GCGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac                 >  1:7326547/1‑101 (MQ=18)
                            aTAAGGCGTTT‑‑GCGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac                 >  2:7326549/1‑101 (MQ=18)
                            aTAAGGCGTTT‑‑GCGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcac                 >  2:7326550/1‑101 (MQ=18)
                             tAAGGCGTTT‑‑GCGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcaca                >  1:7326551/1‑101 (MQ=18)
                             tAAGGCGTTT‑‑GCGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcaca                >  1:7326553/1‑101 (MQ=18)
                             tAAGGCGTTT‑‑GCGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcaca                <  1:7325971/101‑1 (MQ=18)
                             tAAGGCGTTT‑‑GCGCCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTcaca                >  2:7326554/1‑101 (MQ=18)
                                       ttgcgccGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTAt        <  2:7325826/96‑1 (MQ=255)
                                       ttgcgccGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTAt        <  1:7325884/96‑1 (MQ=255)
                                       ttgcgccGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTAt        <  2:7326086/96‑1 (MQ=255)
                                        tgcgccGCATCCGACATCAATGCCTGATGCGACGCTTGTCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTc                   <  1:7326078/85‑1 (MQ=255)
                                        tgcgccGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATc       <  1:7325651/97‑1 (MQ=255)
                                        tgcgccGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATc       <  2:7326093/97‑1 (MQ=255)
                                         gcgccGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCa      <  1:7325849/98‑1 (MQ=255)
                                         gcgccGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCa      >  2:7326570/4‑101 (MQ=255)
                                           gccGCATCCGACATCAATGCCTGATGCGACGTTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCa      >  1:7326574/2‑99 (MQ=255)
                                           gccGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa    >  2:7326573/2‑101 (MQ=255)
                                           gccGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCaaa    >  2:7326572/2‑101 (MQ=255)
                                            ccGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc   >  2:7326578/1‑101 (MQ=255)
                                            ccGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc   >  2:7326577/1‑101 (MQ=255)
                                            ccGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc   <  1:7325429/101‑1 (MQ=255)
                                            ccGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAc   >  1:7326576/1‑101 (MQ=255)
                                             cGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATc       <  2:7325515/96‑1 (MQ=255)
                                             cGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc  <  2:7326231/101‑1 (MQ=255)
                                             cGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc  <  1:7326009/101‑1 (MQ=255)
                                             cGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAAcc  <  1:7325892/101‑1 (MQ=255)
                                         |||||                                                                                                    
CCTGATTTAGCGTGATTTTGTAGGTCGGATAAGGCGTTT‑‑A‑‑CCGCATCCGACATCAATGCCTGATGCGACGCTTGCCGCGTCTTATCAGGCCTATCTTAACCGTTGGTTAATTTTTCAAACTCTTCACACCCGGTATCAAACC  >  NC_000913/4296341‑4296482

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