Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 207,483 A→C 7.4% V786V (GTA→GTC dnaE → DNA polymerase III subunit alpha

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913207,4830AC 7.4% 514.1 / 29.5 175V786V (GTA→GTCdnaEDNA polymerase III subunit alpha
Reads supporting (aligned to +/- strand):  ref base A (82/80);  new base C (5/8);  total (87/88)
Fisher's exact test for biased strand distribution p-value = 5.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.64e-01

Only 100 of 177 total aligned reads displayed.

CGCTGAACTGGCGATGAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGG  >  NC_000913/207336‑207632
                                                                                                                                                   |                                                                                                                                                     
cGCTGAACTGGCGATGAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAAt                                                                                                                                                     >  1:2081217/1‑150 (MQ=255)
    gAACTGGCGATGAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGAcc                                                                                                                                                 >  1:332465/1‑150 (MQ=255)
         ggCGATGAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGa                                                                                                                                            >  2:1844983/1‑150 (MQ=255)
         ggCGATGAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGa                                                                                                                                            <  2:633221/150‑1 (MQ=255)
         ggCGATGAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGa                                                                                                                                            >  2:1597698/1‑150 (MQ=255)
          gCGATGAAAATCTTCGACCTGGTGGAGAAATTCGCTGTTTACGGATTTAACAAATCGCACTCTGAGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGat                                                                                                                                           <  1:1035440/150‑1 (MQ=255)
          gCGATGAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgcc                                                                                                                                           <  1:2225430/150‑3 (MQ=255)
          gCGATGAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGat                                                                                                                                           <  2:1397044/150‑1 (MQ=255)
          gCGATGAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATAGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGat                                                                                                                                           >  2:1777036/1‑150 (MQ=255)
           cGATGAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgcca                                                                                                                                          <  1:1492916/150‑4 (MQ=255)
           cGATGAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGata                                                                                                                                          >  2:2062002/1‑150 (MQ=255)
            gATGAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgccat                                                                                                                                         <  1:2036442/150‑5 (MQ=255)
            gATGAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTCATTAccgccgccat                                                                                                                                         <  1:1423099/150‑5 (MQ=255)
              tGAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATgg                                                                                                                                       <  1:1074426/150‑1 (MQ=255)
               gAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGa                                                                                                                                      <  1:2168409/150‑1 (MQ=255)
                aaaaTCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGac                                                                                                                                     <  1:209082/150‑1 (MQ=255)
                 aaaTCTTCGACCTGGTGGAGAAATTTGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGaca                                                                                                                                    >  2:2351269/1‑150 (MQ=255)
                 aaaTCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGaca                                                                                                                                    >  2:1182659/1‑150 (MQ=255)
                 aaaTCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGaca                                                                                                                                    <  1:1337406/150‑1 (MQ=255)
                 aaaTCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGaca                                                                                                                                    >  2:1220809/1‑150 (MQ=255)
                   aTCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGacaac                                                                                                                                  <  1:128778/150‑1 (MQ=255)
                      ttCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACAcc                                                                                                                               >  2:106486/1‑150 (MQ=255)
                         gcccTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGGGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCgag                                                                                                                            <  2:1054488/148‑1 (MQ=255)
                         gACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCgag                                                                                                                            >  2:1544346/1‑150 (MQ=255)
                          aCCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCgaga                                                                                                                           >  2:2154018/1‑150 (MQ=255)
                               gtggAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAggtg                                                                                                                      <  1:1844983/150‑1 (MQ=255)
                                  gagaAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGgtggtg                                                                                                                   <  2:510477/150‑1 (MQ=255)
                                  gagaAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGgtggtg                                                                                                                   <  1:58586/150‑1 (MQ=255)
                                           gCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtggtg                                                                                                          <  2:1089038/150‑1 (MQ=255)
                                            cTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCtggtgg                                                                                                         <  2:1464823/150‑1 (MQ=255)
                                             tGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGa                                                                                                        <  1:1658103/150‑1 (MQ=255)
                                                ttACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGa                                                                                                     <  1:1597698/150‑1 (MQ=255)
                                                 tACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAg                                                                                                    <  1:1524729/150‑1 (MQ=255)
                                                   cGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAg                                                                                                    <  2:631845/148‑1 (MQ=255)
                                                   cGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAg                                                                                                    >  1:631845/1‑148 (MQ=255)
                                                    ggATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGc                                                                                                 <  1:1220809/150‑1 (MQ=255)
                                                       tttAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTgg                                                                                              <  1:1809226/150‑1 (MQ=255)
                                                         tAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCg                                                                                            >  1:135537/1‑150 (MQ=255)
                                                          aaCAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCgg                                                                                           >  1:2450787/1‑150 (MQ=255)
                                                           aCAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGa                                                                                          >  2:2224060/1‑150 (MQ=255)
                                                            cAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGAt                                                                                         >  2:1732272/1‑150 (MQ=255)
                                                               aTCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATggg                                                                                      >  2:1790804/1‑150 (MQ=255)
                                                                tCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATTGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGTCGGTAATGCCCGCCGATATGGACAACACCGAGAATGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATgggg                                                                                     >  1:2386436/1‑150 (MQ=255)
                                                                tCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATgggg                                                                                     <  1:775126/150‑1 (MQ=255)
                                                                 cGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGc                                                                                    <  1:417989/150‑1 (MQ=255)
                                                                   cACTCTGCGGCCTATGCTTTGGTGTCATATAAAACGTTATGGCTTAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGGGGATGAGTGCGGGCGGATGGGGCTg                                                                                  >  2:826903/1‑150 (MQ=255)
                                                                     ctctGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGaa                                                                                <  1:1468720/150‑1 (MQ=255)
                                                                         gCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGATCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATc                                                                            >  1:58676/1‑150 (MQ=255)
                                                                            gCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCGGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTg                                                                         <  1:869560/150‑1 (MQ=255)
                                                                            gCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTg                                                                         >  1:1796170/1‑150 (MQ=255)
                                                                            gCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTg                                                                         <  1:1777146/150‑1 (MQ=255)
                                                                            gCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTg                                                                         <  1:1029675/150‑1 (MQ=255)
                                                                            gCCTAGGCTTTGGTGGCATATCAAACGGTATGGCTGAAAGCGCACTATCCTGCGGAGGTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGAGGGGGCTGAAAATCCTg                                                                         <  2:1753719/150‑1 (MQ=255)
                                                                             ccTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGc                                                                        >  1:24291/1‑150 (MQ=255)
                                                                             ccTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATACTGc                                                                        >  2:913097/1‑150 (MQ=255)
                                                                               tATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACCACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAAt                                                                             >  1:664770/1‑143 (MQ=255)
                                                                               tATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACCACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAAt                                                                             <  2:664770/143‑1 (MQ=255)
                                                                               tATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGcca                                                                      >  2:1692295/1‑150 (MQ=255)
                                                                               tATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGcca                                                                      >  1:1271821/1‑150 (MQ=255)
                                                                                aTGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGccac                                                                     <  1:1326896/150‑1 (MQ=255)
                                                                                  gCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGccacca                                                                   <  1:504523/150‑1 (MQ=255)
                                                                                        gtgtCATATCAAACGTTCTGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACGGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGAtata                                                             <  2:860773/150‑1 (MQ=255)
                                                                                        gtgtCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGCGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGAtata                                                             >  2:2332817/1‑150 (MQ=255)
                                                                                        gtgtCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGAtata                                                             <  1:378687/150‑1 (MQ=255)
                                                                                          gtCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATaaa                                                           >  2:1907626/1‑150 (MQ=255)
                                                                                          gtCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATaaa                                                           >  2:74984/1‑150 (MQ=255)
                                                                                           tCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAAc                                                          <  1:1761260/150‑1 (MQ=255)
                                                                                            cATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACt                                                         <  1:2060551/150‑1 (MQ=255)
                                                                                            cATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACt                                                         >  1:2340346/1‑150 (MQ=255)
                                                                                            cATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACt                                                         >  2:649461/1‑150 (MQ=255)
                                                                                               atCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCg                                                      >  2:774379/1‑150 (MQ=255)
                                                                                                tCAAACGTTATGTCTGAAAGCGCACTATCCAGCTGAGTTTATGGCGGCCGTAATGACCGCCGATATGGACAACACCGACAAGGTGGTGGGTCTTGTGGATGAGTGCTGGAGGATGGGGCTGAAAATCCTGACACCAGATTTAAACTCCgg                                                     >  2:2457160/1‑150 (MQ=255)
                                                                                                tCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCgg                                                     <  1:1072180/150‑1 (MQ=255)
                                                                                                tCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCgg                                                     >  2:1045657/1‑150 (MQ=255)
                                                                                                  aaaCGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTc                                                   >  1:1463218/1‑150 (MQ=255)
                                                                                                  aaaCGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTc                                                   >  1:949601/1‑150 (MQ=255)
                                                                                                  aaaCGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTc                                                   >  1:784052/1‑150 (MQ=255)
                                                                                                     cGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGTATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCttt                                                >  2:1367294/1‑150 (MQ=255)
                                                                                                      gTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTa                                               <  2:1669070/150‑1 (MQ=255)
                                                                                                           ggCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCAtt                                          >  1:1976745/1‑150 (MQ=255)
                                                                                                           ggCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCAtt                                          >  2:19944/1‑150 (MQ=255)
                                                                                                           ggCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCAtt                                          <  1:1064743/150‑1 (MQ=255)
                                                                                                           ggCAGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGTTCTTTACAAtt                                          >  1:1717133/1‑150 (MQ=255)
                                                                                                             cTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTc                                        >  2:1609151/1‑150 (MQ=255)
                                                                                                              tGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTcc                                       <  1:731937/150‑1 (MQ=255)
                                                                                                                aaaGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCAc                                     <  1:2062002/150‑1 (MQ=255)
                                                                                                                    cgcACTATCCTGCGGAGTTTTTGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCa                                 <  2:2471790/150‑1 (MQ=255)
                                                                                                                    cgcACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCa                                 <  1:1182659/150‑1 (MQ=255)
                                                                                                                     gcACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCaa                                >  2:1160428/1‑150 (MQ=255)
                                                                                                                        cTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAacga                             >  2:1197161/1‑150 (MQ=255)
                                                                                                                        cTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAacga                             <  2:602767/150‑1 (MQ=255)
                                                                                                                              tGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACgg                       <  1:100055/150‑1 (MQ=255)
                                                                                                                                 ggAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGa                    <  1:1197161/150‑1 (MQ=255)
                                                                                                                                      ttATGGCGGCGGTAATGACCGCCGATATGACCACCCCGAGCAGGTGGGGGGTCGGGTGGATGAGTTCGGGCGGGTGGGGGGGCAAATCCGGCCACCCGGTATAAACGCCGGTGGTTACCCTTTCCACGGACACCACGACGGCGCAATccg              >  2:481085/1‑148 (MQ=255)
                                                                                                                                       tATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCgt              <  2:4128/150‑1 (MQ=255)
                                                                                                                                        aTGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGGGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCgtg             <  1:2291707/150‑1 (MQ=255)
                                                                                                                                           gcggcgGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTAt          <  1:1266851/150‑1 (MQ=255)
                                                                                                                                             ggcggTCATTACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATgg        >  2:356448/1‑150 (MQ=255)
                                                                                                                                               cggTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTa      >  1:1987861/1‑150 (MQ=255)
                                                                                                                                                   aaTGACCGCCGATATGGACAACACCGAGAATGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTGCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTgg  <  1:1343253/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
CGCTGAACTGGCGATGAAAATCTTCGACCTGGTGGAGAAATTCGCTGGTTACGGATTTAACAAATCGCACTCTGCGGCCTATGCTTTGGTGTCATATCAAACGTTATGGCTGAAAGCGCACTATCCTGCGGAGTTTATGGCGGCGGTAATGACCGCCGATATGGACAACACCGAGAAGGTGGTGGGTCTGGTGGATGAGTGCTGGCGGATGGGGCTGAAAATCCTGCCACCAGATATAAACTCCGGTCTTTACCATTTCCACGTCAACGACGACGGCGAAATCGTGTATGGTATTGG  >  NC_000913/207336‑207632

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: