Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,040,433 C→A A319D (GCC→GAC)  msrP → protein‑L‑methionine sulfoxide reductase catalytic subunit MsrP

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,040,4330CA95.1% 185.2 / ‑6.1 61A319D (GCC→GAC) msrPprotein‑L‑methionine sulfoxide reductase catalytic subunit MsrP
Reads supporting (aligned to +/- strand):  ref base C (3/0);  new base A (20/38);  total (23/38)
Fisher's exact test for biased strand distribution p-value = 4.92e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCCTCTGGCA  >  NC_000913/2040187‑2040677
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
tGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGACAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGACCGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGGTACGCCGccc                                                                                                                                                                                                                                                    >  2:392923/1‑249 (MQ=255)
  cGACTGATTGTGCCGTGGAAATCTGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAggg                                                                                                                                                                                                                                                >  1:426980/1‑250 (MQ=255)
  cGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGt                                                                                                                                                                                                                                                <  1:367972/251‑1 (MQ=255)
       gATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCGGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTGTAATGGTTACGCCTACCAGGTGGCa                                                                                                                                                                                                                                            <  2:483604/250‑1 (MQ=255)
             gCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGt                                                                                                                                                                                                                                     <  1:22550/251‑1 (MQ=255)
                      atatGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGcc                                                                                                                                                                                                                            <  2:434697/251‑1 (MQ=255)
                        atGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTg                                                                                                                                                                                                                          <  1:12308/251‑1 (MQ=255)
                        aggtcTTTAATGGGATTAATTCGATCGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCTCCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGCTAATCCGTATGTTGATCCCCCGCGCTGTTCACCGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCt                                                                                                                                                                                                                           <  2:131348/246‑1 (MQ=255)
                            ctttacggggATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGATTCTGGCAGCGCCTGACGAATACGGTTTTTACGCCCACGTTAATCCGTATGTTGATCCCCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCATCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATAGCTGTATCGTGGCCTGGAtt                                                                                                                                                                                                                      >  1:435544/8‑251 (MQ=255)
                              ttAAATGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGGCCGCCAACCACCTGGAAGCTGGCATCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGAGCCCCCGCGCTGGTTACAGGCTACCGAACGATTTATTGGTTCAGGCGGTAGCCTCTATGTACAGCGCCAGCCAACGCTACTGTTTAATTGTGACTCCGACCATTTGGCATCGCTGCATCGTGGCCTGGAttt                                                                                                                                                                                                                     <  2:35564/250‑1 (MQ=255)
                                     gATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGGAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGTATTTGCGGgaga                                                                                                                                                                                                             <  2:255255/251‑1 (MQ=255)
                                             cGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTa                                                                                                                                                                                                      >  1:393680/1‑250 (MQ=255)
                                               aTAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAAt                                                                                                                                                                                                   >  1:120548/1‑251 (MQ=255)
                                               aTAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAAt                                                                                                                                                                                                   <  2:393680/251‑1 (MQ=255)
                                                          tAAGCTGACCCTCGAGCGTCCGCCAACCACCTGGAATCTGGCCGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGTCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACa                                                                                                                                                                                         <  2:398590/250‑1 (MQ=255)
                                                            aGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGc                                                                                                                                                                                       >  1:357027/1‑250 (MQ=255)
                                                                gACCCGCGAGCGTCTGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTAGTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGTAGAATTTCTAAATGCGTCTGACAGCaaaa                                                                                                                                                                                   <  2:436064/250‑1 (MQ=255)
                                                                        aGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACTGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTATTGTTTAATGGTTACGTCGACCAGGTGGCACCGCTGTGTCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACa                                                                                                                                                                          <  1:66148/251‑1 (MQ=255)
                                                                                  aaCCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAg                                                                                                                                                                <  1:533158/251‑1 (MQ=255)
                                                                                      cccTGAAATCTGGCAGCGCCTGAAGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTg                                                                                                                                                            <  1:282160/250‑1 (MQ=255)
                                                                                                 tgccgcgcCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCGCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCt                                                                                                                                                   <  2:312827/245‑1 (MQ=255)
                                                                                                     gcgcCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTTTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAAAGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCgg                                                                                                                                             <  1:57486/251‑1 (MQ=255)
                                                                                                      cgcCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGCTCACCCGCGATGGTCACAGGCTAACGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGAACAGGTGGCATCGCTGTATCGTGGCCTGGTTTTGCGGGAGACTTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCGGCATCTTGCCggg                                                                                                                                            >  2:524976/1‑250 (MQ=255)
                                                                                                         cTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGAttg                                                                                                                                         >  1:117005/1‑251 (MQ=255)
                                                                                                            gcGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAAGGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGAttgttg                                                                                                                                      <  2:225882/250‑1 (MQ=255)
                                                                                                              gAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGcc                                                                                                                                    <  1:392923/251‑1 (MQ=255)
                                                                                                                aTACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCTCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGt                                                                                                                                  >  2:6054/1‑251 (MQ=255)
                                                                                                                  aCGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGttt                                                                                                                                >  1:172464/1‑251 (MQ=255)
                                                                                                                  aCGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGttt                                                                                                                                >  1:39420/1‑251 (MQ=255)
                                                                                                                     gTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGATGGCACCAGGCTACCGTCCGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTACTGGTTACTCCGCCCAGGCGGCATCGCTGAATCg                                                                                                                                                                                                                                 >  1:236370/1‑151 (MQ=255)
                                                                                                                             gCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTCCCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGGATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGGTTGATGCCGATTCTCTGGCTGGt                                                                                                                     >  1:61853/1‑251 (MQ=255)
                                                                                                                                       aTCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCtt                                                                                                                                                  <  2:35946/212‑1 (MQ=255)
                                                                                                                                       aTCCGTATGTTGATCACCCGCGCTGGTCACAGGCTAACGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCtt                                                                                                                                                  >  1:35946/1‑212 (MQ=255)
                                                                                                                                          cGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCGCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTTTATCGTGGCCGGGATTTGCGGGATAATTTCTAAAGGCGTCTGACAGCAAATCGGGTGACATGGCGGATAGTTTGGCGGGATCTTTCTGGATGTTTGCCGGTTTGCTGGCTGTTCTGGGCGatcaat                                                                                                        <  2:229515/251‑1 (MQ=255)
                                                                                                                                                 tGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGt                                                                                                  <  2:172464/250‑1 (MQ=255)
                                                                                                                                                        ccGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTGTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGGATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTggg                                                                                          <  2:186377/251‑1 (MQ=255)
                                                                                                                                                            gcTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTTGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGcc                                                                                      <  1:524976/251‑1 (MQ=255)
                                                                                                                                                               ggTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGAt                                                                                   <  2:54146/251‑1 (MQ=255)
                                                                                                                                                                gTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATc                                                                                  <  2:88658/251‑1 (MQ=255)
                                                                                                                                                                      aGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTg                                                                             >  2:416261/1‑251 (MQ=255)
                                                                                                                                                                            ccGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGCTCCGGTGaa                                                                           >  2:367662/1‑247 (MQ=255)
                                                                                                                                                                             cGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCATCGCTACTGTTTAATGGTTTCGCCGACCAGGTGGCGTCGCTGTATCGTGGCCTGGATTTGCTGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGAGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGata                                                                      <  2:148111/251‑1 (MQ=255)
                                                                                                                                                                              gAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGatat                                                                     >  1:243582/1‑251 (MQ=255)
                                                                                                                                                                                 cGATTTATTGGTTCCGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTTGTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCa                                                                  <  2:122337/251‑1 (MQ=255)
                                                                                                                                                                                   agttATTGGTTCAGGCGGCATCCACGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGGCTGACAGCAAAACAGGTGACATGGCTGAAAGTGTTCCTGCAGCTGGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGc                                                                <  2:357027/249‑1 (MQ=255)
                                                                                                                                                                                       aTTGGTTCAGGCGGCATCCTCGATGTACAGCGCCTGCCATCGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTTTATCGTGGCCTGGATTGGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGATGAAAGTTTGCCTGCATCTTGCCGGATTGTCGCCGTTTCTCTGGCTGGTCTGGGTGATCAATCACGGTGGACGGGGTGCCGATCCGGTTAAAGATATTCAGCAttt                                                            <  2:61853/251‑1 (MQ=255)
                                                                                                                                                                                       aTTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTGCGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACTGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCAtt                                                             <  1:225330/250‑1 (MQ=255)
                                                                                                                                                                                          ggttcaggcgtgacacccGATGTACAGTGCCAGCCAACGCTACTGTGTAATGGTTACGCCGATCAGGTGGCATCGCTGGATCGTGGCCTGGATTTGCGGGTGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGAGTGTGGCCGTTGCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTAc                                                         <  2:120548/234‑1 (MQ=255)
                                                                                                                                                                                                ggcggcATCCTCGATGTACATCGCCAGCCAACGCTACGGTTTAATGGGTACTCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGAGTTTCGGGCGAATTTCTAAATGCGTCTGACAGCAAAACAGTTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGTCGTTTCTCTCGCTGGTCTGGTCGATCAATCACGTTTGACTGGGTGCCGATCCCGTAAAAGATATTCATCATGTTACTGGTCg                                                   <  2:474246/251‑1 (MQ=255)
                                                                                                                                                                                                        ccacggtGTACAGCGCCAGCCAACGCTACTGTTTAATGGTGACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGTTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATACGGTGAAAGATATTCAGCATTTTACTGGTCgcactgcac                                           <  2:358869/245‑1 (MQ=255)
                                                                                                                                                                                                           cGATGAACAACGCCAGCCCACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCCTCGCTGTATCGCGGCCTGGATTTGCGGGAGCATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGTGCGATCAATCACGGTGGAATGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGcactgcactg                                         >  2:67722/1‑250 (MQ=255)
                                                                                                                                                                                                            gATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTAc                                                         <  1:422396/233‑1 (MQ=255)
                                                                                                                                                                                                            gATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTAc                                                         >  2:422396/1‑233 (MQ=255)
                                                                                                                                                                                                             aTGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGaa                                       >  2:366172/1‑250 (MQ=255)
                                                                                                                                                                                                              tGTACAGCGCCAACCACCGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGaact                                     >  2:433861/1‑249 (MQ=255)
                                                                                                                                                                                                                   aGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTg                                >  1:77571/1‑251 (MQ=255)
                                                                                                                                                                                                                   aGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTg                                >  1:391359/1‑251 (MQ=255)
                                                                                                                                                                                                                      gccagccaACGCGACTGTTTAATGGTTACGCCGACCAGGTGGCAGCGCTGTATCGTGGCCTGGAGTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATTGCTGAAAGTGTGCCTGCATCTTGCCGGATTGTGGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGGGGACTGGGTGCCGTTCCGGTGAAAGATATTCAGCAGTGTACTGGTCGCACTGAACTGAAATTTTTGCTg                             <  2:435544/251‑1 (MQ=255)
                                                                                                                                                                                                                             aaCGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTGGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTGGCTGTTCGCACTGCACTGAAATTTTTGCTGGCGACCt                      <  1:531949/251‑1 (MQ=255)
                                                                                                                                                                                                                              aCGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCtt                     <  1:6054/251‑1 (MQ=255)
                                                                                                                                                                                                                                        ttAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCAATGAAATTTTTGCTGGCGACCTTGTTAATCAcc           >  2:345720/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 aCGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGGCTGACAGCTAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCAGTTTACTGGTCGTACTGCACTGAAAGTTTTGCTGGCGACCTTGTTAATCCCCCGTCTGGCa  <  2:39420/251‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
TGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCCTCTGGCA  >  NC_000913/2040187‑2040677

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