Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,206,622 A→G 6.4% T310A (ACA→GCA)  yehM → uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,206,6220AG 6.4% 519.7 / 19.5 157T310A (ACA→GCA) yehMuncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (69/78);  new base G (5/5);  total (74/83)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.44e-01

Only 100 of 160 total aligned reads displayed.

GGTGCTGGCAGTGGGGCTTGCAGAAGGCCGGTGAACAACTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGT  >  NC_000913/2206473‑2206771
                                                                                                                                                     |                                                                                                                                                     
ggTGCTGGCAGTGGGGCTTGCAGAAGGCCGGTGAACAACTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATa                                                                                                                                                       <  2:1065928/150‑1 (MQ=255)
     tGGCAGTGGGGCTTGCAGAAGGCCGGTGAACAACTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTa                                                                                                                                                  <  1:2372738/150‑1 (MQ=255)
      ggCAGTGGGGCTTGCAGAAGGCCGGTGAACAACTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTAc                                                                                                                                                 >  2:962144/1‑150 (MQ=255)
      ggCAGTGGGGCTTGCAGAAGGCCGGTGAACAACTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTAc                                                                                                                                                 >  2:962089/1‑150 (MQ=255)
      ggCAGTGGGGCTTGCAGAAGGCCGGTGAACAACTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTAc                                                                                                                                                 >  2:830946/1‑150 (MQ=255)
      ggCAGTGGGGCTTGCAGAAGGCCGGTGAACAACTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTAc                                                                                                                                                 >  2:489446/1‑150 (MQ=255)
       gCAGTGGGGCTTGCAGAAGGCCGGTGAACAACTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTAcc                                                                                                                                                >  2:2008201/1‑150 (MQ=255)
         aGTGGGGCTTGCAGAAGGCCGGTGAACAACTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGt                                                                                                                                              <  1:949367/150‑1 (MQ=255)
           tGGGGCTTGCAGAAGGCCGGTGAACAACTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTa                                                                                                                                            <  1:1075441/150‑1 (MQ=255)
           tGGGGCTTGCAGAAGGCCGGTGAACAACTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTa                                                                                                                                            >  1:518880/1‑150 (MQ=255)
           tGGGGCTTGCAGAAGGCCGGTGAACAACTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTa                                                                                                                                            >  1:1472259/1‑150 (MQ=255)
            ggggCTTGCAGAAGGCCGGTGAACAACTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTAc                                                                                                                                           <  1:1731566/150‑1 (MQ=255)
            ggggCTTGCAGAAGGCCGGTGAACAACTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTAc                                                                                                                                           <  1:90683/150‑1 (MQ=255)
            ggggCTTGCAGAAGGCCGGTGAACAACTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTAc                                                                                                                                           <  1:2110146/150‑1 (MQ=255)
            ggggCTTGCAGAAGGCCGGTGAACAACTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTAc                                                                                                                                           <  2:1352480/150‑1 (MQ=255)
             gggCTTGCAGAAGGCCGGTGAACAACTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACg                                                                                                                                          >  2:2184764/1‑150 (MQ=255)
             gggCTTGCAGAAGGCCGGTGAACAACTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACg                                                                                                                                          >  2:651774/1‑150 (MQ=255)
               gCTTGCAGAAGGCCGGTGAACAACTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCa                                                                                                                                        <  1:1734121/150‑1 (MQ=255)
                 ttGCAGAAGGCCGGTGAACAACTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCact                                                                                                                                      <  1:812719/150‑1 (MQ=255)
                 ttGCAGAAGGCCGGTGAACAACTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCact                                                                                                                                      <  1:1374190/150‑1 (MQ=255)
                          gCCGGTGAACAACTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGActggctg                                                                                                                             >  1:1169469/1‑150 (MQ=255)
                             ggTGAACAACTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGAt                                                                                                                          <  2:2363949/150‑1 (MQ=255)
                               tGAACAACTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGc                                                                                                                        <  2:400199/150‑1 (MQ=255)
                               tGAACAACTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGc                                                                                                                        <  2:459767/150‑1 (MQ=255)
                               tGAACAACTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGc                                                                                                                        <  2:400210/150‑1 (MQ=255)
                                gAACCACTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGcc                                                                                                                       <  2:2092726/150‑1 (MQ=255)
                                  acaacTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATGCACTACCGTTACGCACTGACTGGCTGGATGCCAt                                                                                                                     >  2:1453619/1‑150 (MQ=255)
                                  acaacTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATGCACTACCGTTACGCACTGACTGGCTGGATGCCAt                                                                                                                     >  2:1453548/1‑150 (MQ=255)
                                  acaacTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATGCACTACCGTTACGCACTGACTGGCTGGATGCCAt                                                                                                                     >  2:1453636/1‑150 (MQ=255)
                                    aacTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAg                                                                                                                   >  1:301326/1‑150 (MQ=255)
                                      cTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCa                                                                                                                 <  1:1835091/150‑1 (MQ=255)
                                       tGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAg                                                                                                                >  2:1447803/1‑150 (MQ=255)
                                       tGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAg                                                                                                                >  2:692040/1‑150 (MQ=255)
                                         cTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGc                                                                                                              >  2:995339/1‑150 (MQ=255)
                                            aaaaCTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCg                                                                                                           >  1:1149068/1‑150 (MQ=255)
                                                     cTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTaaa                                                                                                  >  1:1716085/1‑150 (MQ=255)
                                                     cTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTaaa                                                                                                  <  2:984657/150‑1 (MQ=255)
                                                     cTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTaaa                                                                                                  >  1:954997/1‑150 (MQ=255)
                                                     cTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTaaa                                                                                                  >  1:754821/1‑150 (MQ=255)
                                                      ttACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAaag                                                                                                 >  2:794464/1‑150 (MQ=255)
                                                       tACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAaaga                                                                                                >  1:52679/1‑150 (MQ=255)
                                                       tACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAaaga                                                                                                >  1:602423/1‑150 (MQ=255)
                                                        aCCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAaagaa                                                                                               <  2:2238762/150‑1 (MQ=255)
                                                                    cAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAAcgcg                                                                                   >  2:379052/1‑150 (MQ=255)
                                                                    cAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAAcgcg                                                                                   >  1:1516568/1‑150 (MQ=255)
                                                                       cacaAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCg                                                                                >  2:1215011/1‑150 (MQ=255)
                                                                          aaaTTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTg                                                                             <  2:98045/150‑1 (MQ=255)
                                                                            aTTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTTGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGcc                                                                           <  1:680891/150‑1 (MQ=255)
                                                                            aTTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGcc                                                                           <  1:2266154/150‑1 (MQ=255)
                                                                               gCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTg                                                                        <  2:1169469/150‑1 (MQ=255)
                                                                                cccGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTgg                                                                       <  1:2111603/150‑1 (MQ=255)
                                                                                 ccGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGa                                                                      <  1:2311351/150‑1 (MQ=255)
                                                                                 ccGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGa                                                                      <  2:169489/150‑1 (MQ=255)
                                                                                 ccGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGa                                                                      <  2:275959/150‑1 (MQ=255)
                                                                                   gCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGc                                                                    >  2:2155408/1‑150 (MQ=255)
                                                                                     ttCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTa                                                                  >  1:1793406/1‑150 (MQ=255)
                                                                                         aCCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATcgc                                                              <  1:2237561/150‑1 (MQ=255)
                                                                                            gcgGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCgcggc                                                           <  1:352887/150‑1 (MQ=255)
                                                                                             cgGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCgcggcg                                                          <  1:766137/150‑1 (MQ=255)
                                                                                              ggATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGt                                                         <  2:633271/150‑1 (MQ=255)
                                                                                                  atGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTAtt                                                     >  1:1391579/1‑150 (MQ=255)
                                                                                                  atGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTAtt                                                     <  2:1149273/150‑1 (MQ=255)
                                                                                                  atGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCGGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTAtt                                                     >  1:335357/1‑150 (MQ=255)
                                                                                                    ggCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCa                                                   <  2:1779645/150‑1 (MQ=255)
                                                                                                        gCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATccc                                               <  2:311051/150‑1 (MQ=255)
                                                                                                         ccGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCg                                              <  1:1978345/150‑1 (MQ=255)
                                                                                                          cGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGCCTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGa                                             <  2:1052758/150‑1 (MQ=255)
                                                                                                          cGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTGATTCATCCCGa                                             >  2:138555/1‑150 (MQ=255)
                                                                                                            cTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATa                                           <  1:692040/150‑1 (MQ=255)
                                                                                                            cTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATa                                           <  1:1447807/150‑1 (MQ=255)
                                                                                                               aTCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATAccg                                        <  2:2227986/150‑1 (MQ=255)
                                                                                                                tCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATAccga                                       <  2:308306/150‑1 (MQ=255)
                                                                                                                tCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATAccga                                       <  2:1514034/150‑1 (MQ=255)
                                                                                                                tCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATAccga                                       <  2:921553/150‑1 (MQ=255)
                                                                                                                tCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATAccga                                       <  1:2026254/150‑1 (MQ=255)
                                                                                                                  tGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATAccgatc                                     <  1:92373/150‑1 (MQ=255)
                                                                                                                   gCATGCGATGGCTCTGGCACAGTTGCGCGGTCATACACGACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATAccgatcc                                    <  1:1069864/150‑1 (MQ=255)
                                                                                                                   gCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATAccgatcc                                    <  1:2308118/150‑1 (MQ=255)
                                                                                                                    cATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATAccgatccg                                   >  1:1625266/1‑150 (MQ=255)
                                                                                                                     aTGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATAccgatccga                                  >  2:2037493/1‑150 (MQ=255)
                                                                                                                     aTGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATAccgatccga                                  >  1:578297/1‑150 (MQ=255)
                                                                                                                     aTGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATAccgatccga                                  >  1:2354232/1‑150 (MQ=255)
                                                                                                                     aTGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATAccgatccga                                  >  2:1425111/1‑150 (MQ=255)
                                                                                                                     aTGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATAccgatccga                                  >  1:908589/1‑150 (MQ=255)
                                                                                                                     aTGCGATGGCGCTGGCACAGTTGCGCAGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATAccgatccga                                  >  1:1460528/1‑150 (MQ=255)
                                                                                                                       gCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGAtt                                <  1:1776230/150‑1 (MQ=255)
                                                                                                                       gCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGAtt                                <  1:462746/150‑1 (MQ=255)
                                                                                                                       gCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGAtt                                <  1:1344249/150‑1 (MQ=255)
                                                                                                                        cGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTc                               >  2:2308702/1‑150 (MQ=255)
                                                                                                                        cGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTc                               >  2:1543504/1‑150 (MQ=255)
                                                                                                                          aTGGCGCTGGCACAGTTGCGCGGTCATGCACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTg                             <  1:1453587/150‑1 (MQ=255)
                                                                                                                          aTGGCGCTGGCACAGTTGCGCGGTCATGCACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTg                             <  1:1453619/150‑1 (MQ=255)
                                                                                                                          aTGGCGCTGGCACAGTTGCGCGGTCATGCACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTg                             <  1:337398/150‑1 (MQ=255)
                                                                                                                          aTGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAACAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTg                             <  1:1110004/150‑1 (MQ=255)
                                                                                                                              cgcTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTaa                         >  1:2275284/1‑150 (MQ=255)
                                                                                                                               gcTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAAc                        <  1:1272750/150‑1 (MQ=255)
                                                                                                                                       aGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAg                >  1:393121/1‑150 (MQ=255)
                                                                                                                                                   aTACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCgg    >  1:1855614/1‑150 (MQ=255)
                                                                                                                                                   aTACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCgg    >  2:1931093/1‑150 (MQ=255)
                                                                                                                                                     acacTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGt  >  1:994509/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGTGCTGGCAGTGGGGCTTGCAGAAGGCCGGTGAACAACTGCTAAAAACTATCCTTACCCGTTTGCGCCAGCACAAATTGCCCGCTTCTACCGCGGATATGGCTGCCGCTCATCTGCATGCGATGGCGCTGGCACAGTTGCGCGGTCATACACTACCGTTACGCACTGACTGGCTGGATGCCATAGCAGGCTCGCTGATTAAAGAAGCCCTGAACGCGCCGTTGCCGTGGAGCTATCGCGGCGTTATTCATCCCGATACCGATCCGATTCTGCTAACGTTGATAGACACATTAGCGGGT  >  NC_000913/2206473‑2206771

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