Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,131,199 G→C 8.2% S61T (AGT→ACT) ‡ flgB → flagellar basal‑body rod protein FlgB

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,131,1990GC 8.2% 498.2 / 26.0 153S61T (AGT→ACT) ‡flgBflagellar basal‑body rod protein FlgB
Reads supporting (aligned to +/- strand):  ref base G (65/73);  new base C (6/6);  total (72/81)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.29e-01

Only 100 of 158 total aligned reads displayed.

CAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCG  >  NC_000913/1131051‑1131348
                                                                                                                                                    |                                                                                                                                                     
 aacaaGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAgtg                                                                                                                                                     >  1:2335596/1‑150 (MQ=255)
   caaGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGGTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGtgtg                                                                                                                                                   <  1:1058974/150‑1 (MQ=255)
   caaGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGtgtg                                                                                                                                                   >  2:1500418/1‑150 (MQ=255)
        ggCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGc                                                                                                                                              >  2:2472511/1‑150 (MQ=255)
          cgcTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCAc                                                                                                                                            >  1:2345007/1‑150 (MQ=255)
          cgcTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCAc                                                                                                                                            >  2:486755/1‑150 (MQ=255)
          cgcTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCAc                                                                                                                                            <  2:1935092/150‑1 (MQ=255)
          cgcTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCACACTGGTTGCat                                                                                                                                            <  2:1599045/150‑2 (MQ=255)
           gcTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACt                                                                                                                                           <  1:1946514/150‑1 (MQ=255)
            cTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCACACTGGTTGCatcc                                                                                                                                          <  1:1247639/150‑4 (MQ=255)
            cTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCACACTGGTTGCatcc                                                                                                                                          <  1:994733/150‑4 (MQ=255)
            cTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCACACTGGTTGCatcc                                                                                                                                          <  2:1038108/150‑4 (MQ=255)
            cTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCACACTGGTTGCatcc                                                                                                                                          <  2:1628894/150‑4 (MQ=255)
              cAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGAc                                                                                                                                        <  1:1763829/150‑1 (MQ=255)
              cAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGAc                                                                                                                                        >  2:1258632/1‑150 (MQ=255)
               aaTCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACg                                                                                                                                       <  1:1084841/150‑1 (MQ=255)
               aaTCTGCGCGCCCAGCGGCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTGATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACg                                                                                                                                       <  2:2097984/150‑1 (MQ=255)
                aTCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGa                                                                                                                                      >  2:33511/1‑150 (MQ=255)
                  cTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATg                                                                                                                                    <  2:638265/150‑1 (MQ=255)
                   tGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGa                                                                                                                                   >  2:1494331/1‑150 (MQ=255)
                    gcgcgcCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCTCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGAc                                                                                                                                  <  2:1647137/150‑1 (MQ=255)
                    gcgcgcCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGAc                                                                                                                                  >  1:1635077/1‑150 (MQ=255)
                     cgcgcCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGAcc                                                                                                                                 <  2:2157275/150‑1 (MQ=255)
                     cgcgcCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGAcc                                                                                                                                 <  1:1042310/150‑1 (MQ=255)
                     cgcgcCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGAcc                                                                                                                                 >  2:1529646/1‑150 (MQ=255)
                         cccAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCaa                                                                                                                             >  1:1264599/1‑150 (MQ=255)
                           cAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACg                                                                                                                           <  1:486755/150‑1 (MQ=255)
                            aGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGc                                                                                                                          <  1:1251005/150‑1 (MQ=255)
                                tCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAaca                                                                                                                      <  1:33511/150‑1 (MQ=255)
                                     aaGTGCTGGCAGCAAACATCGCCAATGCCGATAGCCATGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTc                                                                                                                 <  1:2207839/150‑1 (MQ=255)
                                      aGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTcc                                                                                                                <  1:365117/150‑1 (MQ=255)
                                         gCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGgc                                                                                                             <  1:799549/150‑1 (MQ=255)
                                         gCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGgc                                                                                                             >  2:113514/1‑150 (MQ=255)
                                                  aaaCATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCt                                                                                                    >  2:1406194/1‑150 (MQ=255)
                                                    aCATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGa                                                                                                  >  2:468553/1‑150 (MQ=255)
                                                     cATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTACACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGAc                                                                                                 <  2:1716448/150‑1 (MQ=255)
                                                               gCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTAcc                                                                                       <  2:114757/150‑1 (MQ=255)
                                                                  gATACCCCTGGTTATCATGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCATTGTCGTTGCACGGACGATGACCTCAACGCAACACATTCAGGCGCAGGCGCTGACGCCTCCTACCGCa                                                                                    <  2:877109/150‑1 (MQ=255)
                                                                    tACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGa                                                                                  >  2:291089/1‑150 (MQ=255)
                                                                      ccccTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAAc                                                                                <  1:2472511/150‑1 (MQ=255)
                                                                      ccccTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAAc                                                                                <  2:1549447/150‑1 (MQ=255)
                                                                      ccccTGGTTATAAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAAc                                                                                >  2:2376732/1‑150 (MQ=255)
                                                                        ccTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTg                                                                              <  1:1490594/150‑1 (MQ=255)
                                                                          tGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCa                                                                            <  2:916602/150‑1 (MQ=255)
                                                                            gTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAAt                                                                          <  2:1275621/150‑1 (MQ=255)
                                                                              tATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATAc                                                                        >  2:313739/1‑150 (MQ=255)
                                                                                  aGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTa                                                                    >  1:185547/1‑150 (MQ=255)
                                                                                   ggCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTAt                                                                   >  1:2228326/1‑150 (MQ=255)
                                                                                       cgcgATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCg                                                               >  1:2283552/1‑150 (MQ=255)
                                                                                           atatCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGAcc                                                           <  1:931687/150‑1 (MQ=255)
                                                                                           atatCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGAcc                                                           >  2:1704657/1‑150 (MQ=255)
                                                                                           atatCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGAcc                                                           <  1:1641007/150‑1 (MQ=255)
                                                                                                gATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATAGCGTATTCCGGACCAGCCt                                                      >  1:1187946/1‑150 (MQ=255)
                                                                                                    ttGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGc                                                  <  2:81086/150‑1 (MQ=255)
                                                                                                     tGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCt                                                 >  1:2256526/1‑150 (MQ=255)
                                                                                                     tGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCt                                                 >  2:715214/1‑150 (MQ=255)
                                                                                                       ccAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTg                                               >  1:2061494/1‑150 (MQ=255)
                                                                                                       ccAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTg                                               <  2:1342903/150‑1 (MQ=255)
                                                                                                        cAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGa                                              >  2:1226199/1‑150 (MQ=255)
                                                                                                         aGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGAc                                             <  1:1733465/150‑1 (MQ=255)
                                                                                                         aGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGAc                                             <  1:203528/150‑1 (MQ=255)
                                                                                                           tGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACgg                                           <  1:1150740/150‑1 (MQ=255)
                                                                                                           tGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACgg                                           <  1:1143321/150‑1 (MQ=255)
                                                                                                               catAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTaaa                                       >  1:772339/3‑149 (MQ=255)
                                                                                                               cTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAAt                                       >  2:468220/1‑150 (MQ=255)
                                                                                                                ttAAAAAAGTCATGCAACGTGGACGGGATTCAAGCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATa                                      <  2:1194764/150‑1 (MQ=255)
                                                                                                                 tAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATAc                                     >  2:701169/1‑150 (MQ=255)
                                                                                                                     aaaGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCATTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTc                                 >  2:2302510/1‑150 (MQ=255)
                                                                                                                      aaGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCg                                <  1:672413/150‑1 (MQ=255)
                                                                                                                       aGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGa                               <  2:1834349/150‑1 (MQ=255)
                                                                                                                        gTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGat                              <  1:1105000/150‑1 (MQ=255)
                                                                                                                        gTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGat                              >  1:2172936/1‑150 (MQ=255)
                                                                                                                        gTCATGCAACGTCGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTATTACCGTcaag                              <  1:698696/150‑4 (MQ=255)
                                                                                                                         tCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTATTACCGTcaagc                             <  1:1863370/150‑5 (MQ=255)
                                                                                                                         gcATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTATTACCGTcaagc                             <  2:1890024/149‑5 (MQ=255)
                                                                                                                          cATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGatat                            <  1:1918234/150‑1 (MQ=255)
                                                                                                                           aTGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATg                           >  2:970194/1‑150 (MQ=255)
                                                                                                                           aTGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATg                           <  2:892901/150‑1 (MQ=255)
                                                                                                                              cAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGGCGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGAt                        >  2:2350107/1‑150 (MQ=255)
                                                                                                                                aCGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATcg                      <  2:688593/150‑1 (MQ=255)
                                                                                                                                    ggACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGaa                  >  1:2369061/1‑150 (MQ=255)
                                                                                                                                    ggACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGaa                  >  2:1108348/1‑150 (MQ=255)
                                                                                                                                      cgtcagtGCAACCACACTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACg                >  2:407468/7‑150 (MQ=255)
                                                                                                                                       gtcagtGCAACCACACTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGc               >  1:1886742/6‑150 (MQ=255)
                                                                                                                                         ggATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCAc             >  2:203919/1‑150 (MQ=255)
                                                                                                                                         ggATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCAc             >  2:612247/1‑150 (MQ=255)
                                                                                                                                         cagtGCAACCACACTGGTTGCGCTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCAc             >  2:686445/4‑150 (MQ=255)
                                                                                                                                         cagtGCAACCACACTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCAc             >  1:803958/4‑150 (MQ=255)
                                                                                                                                          agtGCAACCACACTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCAcc            >  1:2441473/3‑150 (MQ=255)
                                                                                                                                          agtGCAACCACACTGGTTGCACTGACGATGACCTCAACGCAACACAATCCGGCGCAGGCTCTGACGCCTCCTACCGCAGAACTGCAATACCGTAATCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCAcc            >  2:2043696/3‑150 (MQ=255)
                                                                                                                                              cAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTACCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGt        >  2:2247406/1‑150 (MQ=255)
                                                                                                                                               aaCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATTGATCGCGAATGCACCCAGtt       <  2:1967217/150‑1 (MQ=255)
                                                                                                                                               aaCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGtt       >  2:1353836/1‑150 (MQ=255)
                                                                                                                                                aCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGttt      <  1:662691/150‑1 (MQ=255)
                                                                                                                                                aCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGttt      <  1:1500418/150‑1 (MQ=255)
                                                                                                                                                aCCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGttt      >  1:1158605/1‑150 (MQ=255)
                                                                                                                                                 ccAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTATg     >  2:2134033/1‑150 (MQ=255)
                                                                                                                                                   aGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGcc   >  2:874194/1‑150 (MQ=255)
                                                                                                                                                   aGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGcc   <  2:185547/150‑1 (MQ=255)
                                                                                                                                                    gtgtGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCg  >  1:1205764/1‑150 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
CAACAAGAGGCGCTCAATCTGCGCGCCCAGCGTCAGGAAGTGCTGGCAGCAAACATCGCCAATGCCGATACCCCTGGTTATCAGGCGCGCGATATCGATTTTGCCAGTGAACTTAAAAAAGTCATGCAACGTGGACGGGATGCAACCAGTGTGGTTGCACTGACGATGACCTCAACGCAACACATTCCGGCGCAGGCGCTGACGCCTCCTACCGCAGAACTGCAATACCGTATTCCGGACCAGCCTTCGCTTGACGGTAATACCGTCGATATGGATCGCGAACGCACCCAGTTTGCCG  >  NC_000913/1131051‑1131348

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: