Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,174,080 T→C 6.3% intergenic (‑116/+12) mfd ← / ← ycfT transcription‑repair coupling factor/inner membrane protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,174,0800TC 6.3% 1715.2 / 68.6 510intergenic (‑116/+12)mfd/ycfTtranscription‑repair coupling factor/inner membrane protein
Reads supporting (aligned to +/- strand):  ref base T (230/248);  new base C (12/20);  total (242/268)
Fisher's exact test for biased strand distribution p-value = 2.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.24e-01

Only 100 of 513 total aligned reads displayed.

GACGGGCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAG  >  NC_000913/1173932‑1174229
                                                                                                                                                    |                                                                                                                                                     
  cGGGCAGCGTAAAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGAAATGCTGTTCATTATCCTTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTggg                                                                                                                                                    <  1:1379495/150‑1 (MQ=255)
   gggCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTgggg                                                                                                                                                   >  1:3845676/1‑150 (MQ=255)
   gggCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTgggg                                                                                                                                                   >  2:3043355/1‑150 (MQ=255)
    ggCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGTATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTTTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGc                                                                                                                                                  <  1:888570/150‑1 (MQ=255)
      cAGCTTATAACGATTTTGTTCTGGCATACCGGTTGTCAGAGTCTGGTTAGTATTTGCCTCAACATTTGGGGGAATATCTCTTATAAGAAATGCTGTTAATTATCCGTAATCGACATGGTTTAGAAATTTGAATCGCCCCGGCTGGGGCGa                                                                                                                                                <  1:2228602/150‑1 (MQ=255)
      cAGCGTATAACGATATTGTTCCGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGa                                                                                                                                                <  2:6246767/150‑1 (MQ=255)
      cAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCAGCCGGGGCGa                                                                                                                                                <  2:7914256/150‑1 (MQ=255)
       aGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGAt                                                                                                                                               >  2:1541729/1‑150 (MQ=255)
       aGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCAGCCGGGGCGAt                                                                                                                                               <  1:6765865/150‑1 (MQ=255)
        gCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATg                                                                                                                                              <  2:3603206/150‑1 (MQ=255)
        gCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATg                                                                                                                                              <  2:310973/150‑1 (MQ=255)
        gCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCAGCCGGGGCGAtt                                                                                                                                              <  1:6713945/150‑2 (MQ=255)
           tataACGATATTTTTAAGTCATAACGGTTTTAATATTCTCTTTAGTATATTCCACAAAATATGGGGTAATATCTCTGATAAGAAATGCTGTTCATTATCCTTAATCGACATGTTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCt                                                                                                                                           <  1:5487661/150‑1 (MQ=255)
              aaCGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAg                                                                                                                                        >  2:5150758/1‑150 (MQ=255)
                 gATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCgg                                                                                                                                     >  2:3228202/1‑150 (MQ=255)
                  atatTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCggg                                                                                                                                    <  2:1312035/150‑1 (MQ=255)
                       gTTCAGGCATAACGGTTGTCAGTTTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCTAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTaa                                                                                                                               >  1:1864071/1‑150 (MQ=255)
                        ttCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAAc                                                                                                                              >  1:6127795/1‑150 (MQ=255)
                        ttCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAAc                                                                                                                              <  2:656483/150‑1 (MQ=255)
                        ttCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAAc                                                                                                                              >  2:5394383/1‑150 (MQ=255)
                          cAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTg                                                                                                                            >  1:3303539/1‑150 (MQ=255)
                           cggCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGa                                                                                                                           <  1:5013410/149‑1 (MQ=255)
                           aGGCATACCGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGa                                                                                                                           <  1:767940/150‑1 (MQ=255)
                           aGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCAGCCGGGGCGa                                                                                                                                                <  2:3487544/129‑1 (MQ=255)
                           aGGCATAACGGTTGTCAGATTCTCGTTAGGATAGGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTGAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGa                                                                                                                           <  2:6393450/150‑1 (MQ=255)
                            ggCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGAc                                                                                                                          <  2:1785739/150‑1 (MQ=255)
                                tAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCa                                                                                                                      <  1:3740395/150‑1 (MQ=255)
                                       tGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCa                                                                                                               <  2:4574634/150‑1 (MQ=255)
                                        gTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAg                                                                                                              <  1:4806268/150‑1 (MQ=255)
                                         tCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGt                                                                                                             >  2:5277669/1‑150 (MQ=255)
                                          cAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTg                                                                                                            >  1:3183653/1‑150 (MQ=255)
                                            gATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGaa                                                                                                          >  1:6252147/1‑150 (MQ=255)
                                            gATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGaa                                                                                                          >  1:2738650/1‑150 (MQ=255)
                                              ttCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAgg                                                                                                        >  1:248822/1‑150 (MQ=255)
                                                 tcGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCgg                                                                                                     >  2:7295980/1‑150 (MQ=255)
                                                 tcGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATAGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCgg                                                                                                     >  1:7889389/1‑150 (MQ=255)
                                                 tcGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATAGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCgg                                                                                                     >  1:3955936/1‑150 (MQ=255)
                                                  cGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCggg                                                                                                    <  1:5098093/150‑1 (MQ=255)
                                                  cGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCggg                                                                                                    <  1:8176241/150‑1 (MQ=255)
                                                   gTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCgggg                                                                                                   >  1:3060224/1‑150 (MQ=255)
                                                    ttAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGga                                                                                                  >  1:3543580/1‑150 (MQ=255)
                                                    ttAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGga                                                                                                  <  2:6074762/150‑1 (MQ=255)
                                                       ggATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGaa                                                                                               <  2:7521197/150‑1 (MQ=255)
                                                        gATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGaac                                                                                              >  1:8012536/1‑150 (MQ=255)
                                                        gATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGATTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGaac                                                                                              >  2:442972/1‑150 (MQ=255)
                                                               cTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGAt                                                                                       >  1:5196066/1‑150 (MQ=255)
                                                               cTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGAt                                                                                       >  1:1025491/1‑150 (MQ=255)
                                                               cTCAAAATTTGGGGGCATATCTCTGATAAGCAATGCTTTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGAt                                                                                       <  1:1241628/150‑1 (MQ=255)
                                                                 cAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCg                                                                                     >  1:5641435/1‑150 (MQ=255)
                                                                   aCATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCt                                                                                   <  1:6910316/150‑1 (MQ=255)
                                                                     atatGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTCACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGa                                                                                 >  2:4854214/1‑150 (MQ=255)
                                                                       atGGGGGCATATCTCTGATACGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAAt                                                                               <  2:6780624/150‑1 (MQ=255)
                                                                       atGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAAt                                                                               <  2:1119780/150‑1 (MQ=255)
                                                                        tGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATg                                                                              <  1:1585886/150‑1 (MQ=255)
                                                                         gggggCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGc                                                                             >  1:4166270/1‑150 (MQ=255)
                                                                            ggCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACg                                                                          <  2:2918770/150‑1 (MQ=255)
                                                                             gCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGc                                                                         >  1:1102220/1‑150 (MQ=255)
                                                                               atatCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCtg                                                                       <  1:5263959/150‑1 (MQ=255)
                                                                               atatCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCtg                                                                       >  2:3307812/1‑150 (MQ=255)
                                                                                tatCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCtgt                                                                      <  1:7128077/150‑1 (MQ=255)
                                                                                  tctcTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGa                                                                    >  2:8119466/1‑150 (MQ=255)
                                                                                  tctcTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGa                                                                    >  1:641053/1‑150 (MQ=255)
                                                                                       gATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGtt                                                               <  2:4691047/150‑1 (MQ=255)
                                                                                         tAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGtttt                                                             <  2:238189/150‑1 (MQ=255)
                                                                                              aaTGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTAc                                                        <  2:3954959/150‑1 (MQ=255)
                                                                                                tGCTGTTCATTATCCTTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCa                                                      <  1:2959254/150‑1 (MQ=255)
                                                                                                tGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCa                                                      <  1:2626158/150‑1 (MQ=255)
                                                                                                   tGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCa                                                   >  2:4526506/1‑150 (MQ=255)
                                                                                                     ttCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAg                                                 >  1:5549692/1‑150 (MQ=255)
                                                                                                        aTTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGccc                                              >  2:2748380/1‑150 (MQ=255)
                                                                                                         ttATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCg                                             <  1:1441240/150‑1 (MQ=255)
                                                                                                            tCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCaa                                          >  2:7420765/1‑150 (MQ=255)
                                                                                                             ccGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAAc                                         <  1:2662536/150‑1 (MQ=255)
                                                                                                              cGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACa                                        <  1:5802979/150‑1 (MQ=255)
                                                                                                                        aTGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGa                              >  1:1679730/1‑150 (MQ=255)
                                                                                                                         tGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGaa                             <  1:7462089/150‑1 (MQ=255)
                                                                                                                           gTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACa                           >  2:7712984/1‑150 (MQ=255)
                                                                                                                           gTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACa                           >  1:6574933/1‑150 (MQ=255)
                                                                                                                               aGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACt                       <  1:8229758/150‑1 (MQ=255)
                                                                                                                               aGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACt                       <  1:4358619/150‑1 (MQ=255)
                                                                                                                                 cAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGa                     >  2:4995631/1‑150 (MQ=255)
                                                                                                                                 cAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGa                     <  2:2193531/150‑1 (MQ=255)
                                                                                                                                  aaaTTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGAc                    >  2:5679152/1‑150 (MQ=255)
                                                                                                                                   aaTTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACa                   <  2:6475530/150‑1 (MQ=255)
                                                                                                                                   aaTTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACa                   >  2:2071778/1‑150 (MQ=255)
                                                                                                                                   aaTTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACa                   >  2:7189064/1‑150 (MQ=255)
                                                                                                                                   aaTTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACa                   <  1:4679913/150‑1 (MQ=255)
                                                                                                                                   aaTTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACa                   >  1:4587084/1‑150 (MQ=255)
                                                                                                                                        aaTCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAggg              >  2:3824508/1‑150 (MQ=255)
                                                                                                                                        aaTCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAggg              <  1:1763025/150‑1 (MQ=255)
                                                                                                                                        aaTCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAggg              <  1:5150758/150‑1 (MQ=255)
                                                                                                                                        aaTCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAggg              >  2:1906872/1‑150 (MQ=255)
                                                                                                                                         aTCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGa             >  1:3371172/1‑150 (MQ=255)
                                                                                                                                          tCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGAt            >  1:2399572/1‑150 (MQ=255)
                                                                                                                                               ccGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGAACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGAACAGATGAACAAACTGACAAAGGGATAAAcc       >  2:5782567/1‑150 (MQ=255)
                                                                                                                                                cGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACca      <  1:6009156/150‑1 (MQ=255)
                                                                                                                                                 ggCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACcag     <  2:4137523/150‑1 (MQ=255)
                                                                                                                                                 ggCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACcag     <  1:998919/150‑1 (MQ=255)
                                                                                                                                                    tGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACcagcag  >  1:212295/1‑150 (MQ=255)
                                                                                                                                                    tGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACcagcag  <  2:255674/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GACGGGCAGCGTATAACGATATTGTTCAGGCATAACGGTTGTCAGATTCTCGTTAGGATATGCCTCAACATATGGGGGCATATCTCTGATAAGCAATGCTGTTCATTATCCGTAATCGACATGGTTTAGCAAATTGAATCGCCCCGGCTGGGGCGATGGCTTAGCGGGAGTAACTGACGGCAGCGGGCAGTGAAGGCGGGGAGAACAACAGATCGCTGAATGCACGCTGTGAAAGTTTTCTTACCAGCAAGCCCGCAACAGTACAGATGAACAAACTGACAAAGGGATAAACCAGCAG  >  NC_000913/1173932‑1174229

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