Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,794,088 T→A 5.8% intergenic (+73/‑165) gabT → / → gabP 4‑aminobutyrate aminotransferase GabT/4‑aminobutyrate:H(+) symporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,794,0880TA 5.8% 389.6 / 7.8 119intergenic (+73/‑165)gabT/gabP4‑aminobutyrate aminotransferase GabT/4‑aminobutyrate:H(+) symporter
Reads supporting (aligned to +/- strand):  ref base T (54/58);  new base A (3/4);  total (57/62)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

Only 100 of 121 total aligned reads displayed.

CGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGTGACACACTT  >  NC_000913/2793939‑2794237
                                                                                                                                                     |                                                                                                                                                     
cGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGttt                                                                                                                                                       >  1:634831/1‑150 (MQ=255)
 gCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTa                                                                                                                                                      >  2:76312/1‑150 (MQ=255)
 gCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTa                                                                                                                                                      >  2:1260400/1‑150 (MQ=255)
  cTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAc                                                                                                                                                     >  1:687438/1‑150 (MQ=255)
    cACCATTGAAGCCGCTCAGATCCGTCAGGGTCTGGAGATAACCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcgc                                                                                                                                                   >  1:161059/1‑150 (MQ=255)
     cccATTTAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGTAGAAGCAGTATCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCATGTCGGATAAGGCGTTTAcgcc                                                                                                                                                  <  2:1096692/149‑1 (MQ=255)
       cATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTAcgccgc                                                                                                                                                <  1:1156057/150‑1 (MQ=255)
           gAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                            <  1:1028181/150‑1 (MQ=255)
           gAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATcc                                                                                                                                            <  1:752779/150‑1 (MQ=255)
           gAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATcc                                                                                                                                            >  1:356608/1‑150 (MQ=255)
           gAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATcc                                                                                                                                            >  2:1476281/1‑150 (MQ=255)
             aGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGa                                                                                                                                          >  2:571182/1‑150 (MQ=255)
               aCGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATGCGGCGTAAACGCCTTATCCGAcc                                                                                                                                        <  1:1490088/150‑2 (MQ=255)
                cGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACaa                                                                                                                                       >  2:90773/1‑150 (MQ=255)
                  cTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGAc                                                                                                                                         >  2:474578/1‑146 (MQ=255)
                  cTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGAc                                                                                                                                         <  1:474578/146‑1 (MQ=255)
                  cTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATc                                                                                                                                     <  1:1351932/150‑1 (MQ=255)
                  cTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATc                                                                                                                                     >  2:532125/1‑150 (MQ=255)
                  cTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATc                                                                                                                                     >  2:1144073/1‑150 (MQ=255)
                    cAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTg                                                                                                                                   >  2:359944/1‑150 (MQ=255)
                     aGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGt                                                                                                                                  <  1:1031898/150‑1 (MQ=255)
                      gATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTa                                                                                                                                 >  1:87386/1‑150 (MQ=255)
                        tCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACg                                                                                                                               <  2:767835/150‑1 (MQ=255)
                           gTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGa                                                                                                                            <  1:1084788/150‑1 (MQ=255)
                           gTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGa                                                                                                                            <  2:692095/150‑1 (MQ=255)
                              aGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCTGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACa                                                                                                                         <  1:693787/150‑1 (MQ=255)
                                ggTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAgg                                                                                                                       >  1:1253398/1‑150 (MQ=255)
                                 gTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGga                                                                                                                      >  1:1294614/1‑150 (MQ=255)
                                 gTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGga                                                                                                                      >  1:1230736/1‑150 (MQ=255)
                                   cTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGgaag                                                                                                                    <  1:783395/150‑1 (MQ=255)
                                    tGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGgaaga                                                                                                                   >  1:230973/1‑150 (MQ=255)
                                    tGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGgaaga                                                                                                                   >  2:285980/1‑150 (MQ=255)
                                    tGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGgaaga                                                                                                                   <  1:441140/150‑1 (MQ=255)
                                     ggAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGgaagaa                                                                                                                  >  2:227449/1‑150 (MQ=255)
                                        cggcatcaGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATc                                                                                                               <  2:197454/147‑1 (MQ=255)
                                             tcaGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGt                                                                                                          <  1:1125668/150‑1 (MQ=255)
                                              caGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGtt                                                                                                         >  2:473962/1‑150 (MQ=255)
                                              caGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGtt                                                                                                         >  1:231553/1‑150 (MQ=255)
                                               aGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTg                                                                                                        <  2:276366/150‑1 (MQ=255)
                                                     tgtTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCggg                                                                                                  <  2:675625/150‑1 (MQ=255)
                                                      gtTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGt                                                                                                 >  1:1276844/1‑150 (MQ=255)
                                                         ttGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAg                                                                                              <  2:784658/150‑1 (MQ=255)
                                                         ttGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCAGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAACCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAg                                                                                              >  1:99822/1‑150 (MQ=255)
                                                          tgatgaGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCTGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAgg                                                                                             <  2:325849/150‑1 (MQ=255)
                                                           gatgaGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGc                                                                                            <  1:803296/150‑1 (MQ=255)
                                                           gatgaGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGc                                                                                            >  2:1399541/1‑150 (MQ=255)
                                                           gatgaGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGc                                                                                            >  2:572796/1‑150 (MQ=255)
                                                            atgaGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCg                                                                                           <  2:1128619/150‑1 (MQ=255)
                                                             tgaGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCgg                                                                                          >  1:1359774/1‑150 (MQ=255)
                                                              gaGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGt                                                                                                 >  2:1488244/1‑142 (MQ=255)
                                                              gaGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGt                                                                                                 <  1:1488244/142‑1 (MQ=255)
                                                                ggCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGc                                                                                       >  1:1148613/1‑150 (MQ=255)
                                                                ggCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGc                                                                                       >  1:1488895/1‑150 (MQ=255)
                                                                 gCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGcc                                                                                      <  2:687438/150‑1 (MQ=255)
                                                                 gCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGcc                                                                                      <  2:1148613/150‑1 (MQ=255)
                                                                  cGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCg                                                                                     >  1:902510/1‑150 (MQ=255)
                                                                   gAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGc                                                                                    <  2:666907/150‑1 (MQ=255)
                                                                   gAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGGCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGc                                                                                    <  1:1487579/150‑1 (MQ=255)
                                                                     aGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGcta                                                                                  >  2:488152/1‑149 (MQ=255)
                                                                        aGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTcc                                                                               >  2:1423098/1‑150 (MQ=255)
                                                                          tAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCgg                                                                             >  2:520067/1‑150 (MQ=255)
                                                                            gcgcCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCa                                                                           <  1:76312/150‑1 (MQ=255)
                                                                              gccgcTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCaaa                                                                         <  2:1253398/150‑1 (MQ=255)
                                                                               ccgcTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCaaaa                                                                        >  1:1091804/1‑150 (MQ=255)
                                                                                 gcTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAg                                                                      <  2:1276844/150‑1 (MQ=255)
                                                                                   tCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGaa                                                                    <  1:473962/150‑1 (MQ=255)
                                                                                    ccTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAAt                                                                   <  1:1423098/150‑1 (MQ=255)
                                                                                    ccTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAAt                                                                   <  1:1260400/150‑1 (MQ=255)
                                                                                      tATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCa                                                                 >  2:172763/1‑150 (MQ=255)
                                                                                       aTGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCaa                                                                <  1:1600197/150‑1 (MQ=255)
                                                                                            ggAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACa                                                           >  2:240479/1‑150 (MQ=255)
                                                                                              aGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAAt                                                         <  2:87386/150‑1 (MQ=255)
                                                                                               ggCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAAtt                                                        >  2:629198/1‑150 (MQ=255)
                                                                                                 cGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATtat                                                      >  1:726214/1‑150 (MQ=255)
                                                                                                      cTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATa                                                              >  2:318854/1‑137 (MQ=255)
                                                                                                      cTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATAcgcg                                                 <  1:594443/150‑1 (MQ=255)
                                                                                                          gcgTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGAc                                             <  1:625438/150‑1 (MQ=255)
                                                                                                              cTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCgg                                         >  2:55206/1‑150 (MQ=255)
                                                                                                                 gTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGc                                                                                    >  1:1624970/1‑104 (MQ=255)
                                                                                                                 gTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGc                                                                                    <  2:1624970/104‑1 (MQ=255)
                                                                                                                 gTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGgcgc                                      >  1:17829/1‑150 (MQ=255)
                                                                                                                        cTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAAt                               <  1:1144073/150‑1 (MQ=255)
                                                                                                                          aCGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCgg                                    >  1:774218/1‑143 (MQ=255)
                                                                                                                          aCGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCgg                                    <  2:774218/143‑1 (MQ=255)
                                                                                                                          aCGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGt                             <  1:240479/150‑1 (MQ=255)
                                                                                                                           cGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTc                            <  1:1109608/150‑1 (MQ=255)
                                                                                                                           cGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTc                            >  1:1359428/1‑150 (MQ=255)
                                                                                                                              ggATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCggg                         >  1:519079/1‑150 (MQ=255)
                                                                                                                               gATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCgggg                        <  2:948589/150‑1 (MQ=255)
                                                                                                                                    ttgTCGGATGCGGCGTAAACGCCTTATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGctc                   >  2:182761/3‑150 (MQ=255)
                                                                                                                                      gTCGGATGCGGCGTAAACGCCTTATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCtctc                 >  1:1614816/1‑150 (MQ=255)
                                                                                                                                        cGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTccc               <  2:726214/150‑1 (MQ=255)
                                                                                                                                         ggATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCa              >  1:1261033/1‑150 (MQ=255)
                                                                                                                                         ggATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCa              >  1:303349/1‑150 (MQ=255)
                                                                                                                                           aTAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAg            <  1:227449/150‑1 (MQ=255)
                                                                                                                                           aTAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAg            <  2:1230736/150‑1 (MQ=255)
                                                                                                                                              cggCGTAAACGCCTTATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCgg                                         <  2:929253/117‑1 (MQ=255)
                                                                                                                                               ggCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGTGac        <  1:1476281/150‑1 (MQ=255)
                                                                                                                                                gCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGTGaca       <  2:1359428/150‑1 (MQ=255)
                                                                                                                                                     tACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCAGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATCTCGGGGCGCTCTCCCAAGTGACACACtt  >  1:703053/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
CGCTCACCATTGAAGACGCTCAGATCCGTCAGGGTCTGGAGATCATCAGCCAGTGTTTTGATGAGGCGAAGCAGTAGCGCCGCTCCTATGCCGGAGGCGACGCTGCGCGTCTTGTCCGGCCTACGGGGATCCAGGTCGGATAAGGCGTTTACGCCGCATCCGACAATCTGTACGTGAACAGGAAGAAATCTATGTTGGCCGGGTAAGGCGGAGCCGCTCTCCGGCAAAAAGAATCAATAACAATTATACGCGTGACCCGGCGCGGGAAATGTCGGGGCGCTCTCCCAAGTGACACACTT  >  NC_000913/2793939‑2794237

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