Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,716,033 C→A 9.3% intergenic (+23/+33) eamB → / ← grcA cysteine/O‑acetylserine exporter EamB/stress‑induced alternate pyruvate formate‑lyase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,716,0330CA 9.3% 276.5 / 8.1 85intergenic (+23/+33)eamB/grcAcysteine/O‑acetylserine exporter EamB/stress‑induced alternate pyruvate formate‑lyase subunit
Reads supporting (aligned to +/- strand):  ref base C (34/43);  new base A (6/2);  total (40/45)
Fisher's exact test for biased strand distribution p-value = 1.40e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.62e-01

GTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCAGCTGCGGATATTTTTCCGG  >  NC_000913/2715884‑2716182
                                                                                                                                                     |                                                                                                                                                     
gTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGc                                                                                                                                                       >  1:134781/1‑150 (MQ=255)
 tttGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGcc                                                                                                                                                      >  1:7767/1‑150 (MQ=255)
  ttGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGccc                                                                                                                                                     >  2:263900/1‑150 (MQ=255)
       aaTGTGTGCTGGGCGCTGGCGGGGGAGCTGTTTCAGCGTTTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTc                                                                                                                                                <  2:677866/150‑1 (MQ=255)
          gtgtgCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGGCGACCTCTTCCGc                                                                                                                                             <  1:523301/150‑1 (MQ=255)
           tgtgCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCt                                                                                                                                            >  2:1534316/1‑150 (MQ=255)
            gtgCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCtt                                                                                                                                           >  1:1502415/1‑150 (MQ=255)
              gCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCTTTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAg                                                                                                                                         <  1:560903/150‑1 (MQ=255)
              gCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAg                                                                                                                                         >  1:243521/1‑150 (MQ=255)
               cTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAACTATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGt                                                                                                                                        <  1:167428/150‑1 (MQ=255)
                tcggCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTa                                                                                                                                       >  1:1251657/3‑150 (MQ=255)
                   gcgcTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACt                                                                                                                                    <  1:1082451/150‑1 (MQ=255)
                         gCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTAc                                                                                                                              <  1:1426231/150‑1 (MQ=255)
                         gCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTAc                                                                                                                              <  2:1502415/150‑1 (MQ=255)
                             ggCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTaa                                                                                                                          >  1:938211/1‑150 (MQ=255)
                              gCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTatg                                                                                                                         <  1:925558/150‑3 (MQ=255)
                                   tGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTa                                                                                                                    <  2:1076851/150‑1 (MQ=255)
                                    gTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTAc                                                                                                                   <  2:134781/150‑1 (MQ=255)
                                    gTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTAc                                                                                                                   >  2:1631596/1‑150 (MQ=255)
                                        cAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGc                                                                                                               >  1:151385/1‑150 (MQ=255)
                                         aGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCt                                                                                                              >  1:316499/1‑150 (MQ=255)
                                          gCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCtt                                                                                                             <  2:938211/150‑1 (MQ=255)
                                           cGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCttt                                                                                                            >  1:701255/1‑150 (MQ=255)
                                            gATTGTTTCGCCAGTATGGTCGCCATTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTc                                                                                                           <  1:204854/150‑1 (MQ=255)
                                            gATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTc                                                                                                           <  1:778418/150‑1 (MQ=255)
                                            gATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTc                                                                                                           <  1:678548/150‑1 (MQ=255)
                                            gATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGGCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTc                                                                                                           <  1:1895445/150‑1 (MQ=255)
                                              ttGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAg                                                                                                         <  1:244963/150‑1 (MQ=255)
                                              ttGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAg                                                                                                         <  2:316499/150‑1 (MQ=255)
                                               tGTTTCGCCTGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGt                                                                                                        <  2:908201/150‑1 (MQ=255)
                                               tGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGt                                                                                                        <  1:381579/150‑1 (MQ=255)
                                                gTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTa                                                                                                       >  2:1932933/1‑150 (MQ=255)
                                                gTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTa                                                                                                       >  1:560200/1‑150 (MQ=255)
                                                   tCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCTGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAg                                                                                                    <  2:1634536/150‑1 (MQ=255)
                                                   tCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTGGCGCTCTTCGGCTTAGTAACTTGCTACTTATGCTTTACAGGCTTTCAGTAAAg                                                                                                    <  2:1627533/150‑1 (MQ=255)
                                                     gCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGc                                                                                                                        <  1:513343/128‑1 (MQ=255)
                                                     gCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGc                                                                                                                        >  2:513343/1‑128 (MQ=255)
                                                     gCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGt                                                                                                  >  1:295695/1‑150 (MQ=255)
                                                      ccAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTa                                                                                                 <  1:1534316/150‑1 (MQ=255)
                                                      ccAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTCTTGCGGGTCTATTGAGCGGTACCCATTTTCTATTAACGATAAAAAGCTGAAGAGCTCGCCCTCTTCAGCTTAGTAACTTGCTACTTAAGCCATACAGCCTTTCAGTAAAggaa                                                                                                 >  2:1020976/1‑148 (MQ=255)
                                                           aTGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGc                                                                                            >  2:1273615/1‑150 (MQ=255)
                                                               tCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGACGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGGACGAGCGATa                                                                                        <  2:1688427/150‑1 (MQ=255)
                                                                cGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGGAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATaa                                                                                       >  2:1157621/1‑150 (MQ=255)
                                                                  ccAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACg                                                                                     <  2:1560671/150‑1 (MQ=255)
                                                                     gTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTcg                                                                                  <  2:151385/150‑1 (MQ=255)
                                                                     gTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTcg                                                                                  >  1:87055/1‑150 (MQ=255)
                                                                           tatCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCgctgc                                                                            <  1:1032238/150‑1 (MQ=255)
                                                                           tatCGTGCTTGCCCTGTTGCTGGTCTATGGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACATGCTTTCAGTAAAGGTACGAGCGATAACGTCGCgctgc                                                                            <  2:1776613/150‑1 (MQ=255)
                                                                                 gCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTcc                                                                      <  2:1887641/150‑1 (MQ=255)
                                                                                 gCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTcc                                                                      >  2:1829634/1‑150 (MQ=255)
                                                                                 gCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTcc                                                                      >  1:1735546/1‑150 (MQ=255)
                                                                                      cccTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGt                                                                 >  2:20630/1‑150 (MQ=255)
                                                                                         tattGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCag                                                              <  1:958686/148‑1 (MQ=255)
                                                                                           ttGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCagag                                                            <  1:273274/150‑1 (MQ=255)
                                                                                            tGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCagaga                                                           <  1:1940378/150‑1 (MQ=255)
                                                                                              cTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGt                                                         >  2:563011/1‑150 (MQ=255)
                                                                                                  tCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTaaa                                                     <  2:1324151/150‑1 (MQ=255)
                                                                                                    tATTGCGGGGTTCGTATTTTCTATTAACGAAAAAAAGCGTAAGAGGTCGCCGTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGc                                                   <  2:1628582/150‑1 (MQ=255)
                                                                                                    tATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGc                                                   <  1:198127/150‑1 (MQ=255)
                                                                                                       tGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGaa                                                >  1:1000714/1‑150 (MQ=255)
                                                                                                        gcgcGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAAc                                               >  1:7162/1‑150 (MQ=255)
                                                                                                           cgGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGc                                            <  1:1273615/150‑1 (MQ=255)
                                                                                                                  cATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCg                                     <  2:7162/150‑1 (MQ=255)
                                                                                                                    ttttCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGa                                   <  1:1631596/150‑1 (MQ=255)
                                                                                                                     tttCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGAt                                  <  1:263900/150‑1 (MQ=255)
                                                                                                                     tttCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGAt                                  <  2:1125237/150‑1 (MQ=255)
                                                                                                                           ttAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACgg                            <  1:676106/150‑1 (MQ=255)
                                                                                                                             aaCGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGAt                          >  1:1254563/1‑150 (MQ=255)
                                                                                                                             aaCGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGAt                          >  1:1776621/1‑150 (MQ=255)
                                                                                                                              aCGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATg                         >  2:1115814/1‑150 (MQ=255)
                                                                                                                               cGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATgg                        >  1:835988/1‑150 (MQ=255)
                                                                                                                               cGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATgg                        >  1:982775/1‑150 (MQ=255)
                                                                                                                                gAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGt                       >  2:1057961/1‑150 (MQ=255)
                                                                                                                                  tactaaGCGGAAGAGGGCGACCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCa                     >  2:1682016/5‑150 (MQ=255)
                                                                                                                                  aaaaaaGCGGAAGAGGGCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTcc                     >  1:662447/1‑149 (MQ=255)
                                                                                                                                     taaGCGGAAGAGGGCGACCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCAGCt                  >  2:289472/2‑150 (MQ=255)
                                                                                                                                       aGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCAGCTGc                >  1:364844/1‑150 (MQ=255)
                                                                                                                                       aGCGGAAGAGGGCGACCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCAGCTGc                >  1:793571/1‑150 (MQ=255)
                                                                                                                                        gCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCAGCTGCg               >  1:745636/1‑150 (MQ=255)
                                                                                                                                        gCGGAAGAGGGCGACCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCAGCTGCg               >  2:666532/1‑150 (MQ=255)
                                                                                                                                        gCGGAAGAGGGCGACCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCAGCTGCg               >  2:747540/1‑150 (MQ=255)
                                                                                                                                         cGGAAGAGGGCGACCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCAGCTGCgg              <  2:358559/150‑1 (MQ=255)
                                                                                                                                          ggAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCAGCTGCggt             >  2:1749580/1‑149 (MQ=255)
                                                                                                                                             agagGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGGCAGCTGCGGatat          <  1:20630/150‑1 (MQ=255)
                                                                                                                                              gagGGCGACCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCAGCTGCGGATAtt         >  1:1617926/1‑150 (MQ=255)
                                                                                                                                                ggTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCAGCTGCGGATAtttt       >  2:835136/1‑150 (MQ=255)
                                                                                                                                                   cGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCAGCTGCGGATATTTTTcc    <  1:1115814/150‑1 (MQ=255)
                                                                                                                                                    gCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCAGCTGCGGATATTTTTCCg   <  1:563011/150‑1 (MQ=255)
                                                                                                                                                     cccTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCAGCTGCGGATATTTTTCCgg  >  1:1559134/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GTTTGGCAATGTGTGCTGGGCGCTGGCGGGGCATCTGTTTCAGCGATTGTTTCGCCAGTATGGTCGCCAGTTAAATATCGTGCTTGCCCTGTTGCTGGTCTATTGCGCGGTACGCATTTTCTATTAACGAAAAAAAGCGGAAGAGGTCGCCCTCTTCCGCTTAGTAACTTGCTACTTAAGCCTTACAGGCTTTCAGTAAAGGTACGAGCGATAACGTCGCGCTGCTGTTCCGGAGTCAGAGAGTTAAAGCGAACTGCATAACCGGATACACGGATGGTCAGCTGCGGATATTTTTCCGG  >  NC_000913/2715884‑2716182

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