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

TGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCCTCTGGCA  >  NC_000913/2040187‑2040677
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
tGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGACAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCAC‑AGGCTACCGACCGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGGTACGCCGccc                                                                                                                                                                                                                                                    >  2:392923/1‑249 (MQ=255)
  cGACTGATTGTGCCGTGGAAATCTGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAggg                                                                                                                                                                                                                                                >  1:426980/1‑250 (MQ=255)
  cGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGt                                                                                                                                                                                                                                                <  1:367972/251‑1 (MQ=255)
       gATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCGGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTGTAATGGTTACGCCTACCAGGTGGCa                                                                                                                                                                                                                                            <  2:483604/250‑1 (MQ=255)
             gCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGt                                                                                                                                                                                                                                     <  1:22550/251‑1 (MQ=255)
                      atatGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGcc                                                                                                                                                                                                                            <  2:434697/251‑1 (MQ=255)
                        atGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTg                                                                                                                                                                                                                          <  1:12308/251‑1 (MQ=255)
                        aggtcTTTAATGGGATTAATTCGATCGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCTCCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGCTAATCCGTATGTTGATCCCCCGCGCTGTTCAC‑CGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCt                                                                                                                                                                                                                           <  2:131348/246‑1 (MQ=255)
                            ctttacggggATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGATTCTGGCAGCGCCTGACGAATACGGTTTTTACGCCCACGTTAATCCGTATGTTGATCCCCCGCGCTGGTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCATCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATAGCTGTATCGTGGCCTGGAtt                                                                                                                                                                                                                      >  1:435544/8‑251 (MQ=255)
                              ttAAATGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGGCCGCCAACCACCTGGAAGCTGGCATCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGAGCCCCCGCGCTGGTTAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGTAGCCTCTATGTACAGCGCCAGCCAACGCTACTGTTTAATTGTGACTCCGACCATTTGGCATCGCTGCATCGTGGCCTGGAttt                                                                                                                                                                                                                     <  2:35564/250‑1 (MQ=255)
                                     gATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGGAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGTATTTGCGGgaga                                                                                                                                                                                                             <  2:255255/251‑1 (MQ=255)
                                             cGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTa                                                                                                                                                                                                      >  1:393680/1‑250 (MQ=255)
                                               aTAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAAt                                                                                                                                                                                                   >  1:120548/1‑251 (MQ=255)
                                               aTAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAAt                                                                                                                                                                                                   <  2:393680/251‑1 (MQ=255)
                                                          tAAGCTGACCCTCGAGCGTCCGCCAACCACCTGGAATCTGGCCGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCAC‑AGTCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACa                                                                                                                                                                                         <  2:398590/250‑1 (MQ=255)
                                                            aGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGc                                                                                                                                                                                       >  1:357027/1‑250 (MQ=255)
                                                                gACCCGCGAGCGTCTGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCAC‑AGGCTACCGAACGATTTAGTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGTAGAATTTCTAAATGCGTCTGACAGCaaaa                                                                                                                                                                                   <  2:436064/250‑1 (MQ=255)
                                                                        aGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACTGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTATTGTTTAATGGTTACGTCGACCAGGTGGCACCGCTGTGTCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACa                                                                                                                                                                          <  1:66148/251‑1 (MQ=255)
                                                                                  aaCCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAg                                                                                                                                                                <  1:533158/251‑1 (MQ=255)
                                                                                      cccTGAAATCTGGCAGCGCCTGAAGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTg                                                                                                                                                            <  1:282160/250‑1 (MQ=255)
                                                                                                 tgccgcgcCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCGCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCt                                                                                                                                                   <  2:312827/245‑1 (MQ=255)
                                                                                                     gcgcCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTTTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAAAGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCgg                                                                                                                                             <  1:57486/251‑1 (MQ=255)
                                                                                                      cgcCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGCTCACCCGCGATGGTCAC‑AGGCTAACGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGAACAGGTGGCATCGCTGTATCGTGGCCTGGTTTTGCGGGAGACTTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCGGCATCTTGCCggg                                                                                                                                            >  2:524976/1‑250 (MQ=255)
                                                                                                         cTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGAttg                                                                                                                                         >  1:117005/1‑251 (MQ=255)
                                                                                                            gcGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAAGGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGAttgttg                                                                                                                                      <  2:225882/250‑1 (MQ=255)
                                                                                                              gAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGcc                                                                                                                                    <  1:392923/251‑1 (MQ=255)
                                                                                                                aTACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCTCTGGTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGt                                                                                                                                  >  2:6054/1‑251 (MQ=255)
                                                                                                                  aCGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGttt                                                                                                                                >  1:172464/1‑251 (MQ=255)
                                                                                                                  aCGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGttt                                                                                                                                >  1:39420/1‑251 (MQ=255)
                                                                                                                     gTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGATGG‑CACCAGGCTACCGTCCGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTACTGGTTACTCCGCCCAGGCGGCATCGCTGAATCg                                                                                                                                                                                                                                 >  1:236370/1‑151 (MQ=255)
                                                                                                                             gCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCAC‑AGGCTCCCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGGATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGGTTGATGCCGATTCTCTGGCTGGt                                                                                                                     >  1:61853/1‑251 (MQ=255)
                                                                                                                                       aTCCGTATGTTGATCACCCGCGCTGGTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCtt                                                                                                                                                  <  2:35946/212‑1 (MQ=255)
                                                                                                                                       aTCCGTATGTTGATCACCCGCGCTGGTCAC‑AGGCTAACGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCtt                                                                                                                                                  >  1:35946/1‑212 (MQ=255)
                                                                                                                                          cGTATGTTGATCACCCGCGCTGGTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCGCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTTTATCGTGGCCGGGATTTGCGGGATAATTTCTAAAGGCGTCTGACAGCAAATCGGGTGACATGGCGGATAGTTTGGCGGGATCTTTCTGGATGTTTGCCGGTTTGCTGGCTGTTCTGGGCGatcaat                                                                                                        <  2:229515/251‑1 (MQ=255)
                                                                                                                                                 tGATCACCCGCGCTGGTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGt                                                                                                  <  2:172464/250‑1 (MQ=255)
                                                                                                                                                        ccGCGCTGGTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTGTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGGATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTggg                                                                                          <  2:186377/251‑1 (MQ=255)
                                                                                                                                                            gcTGGTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTTGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGcc                                                                                      <  1:524976/251‑1 (MQ=255)
                                                                                                                                                               ggTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGAt                                                                                   <  2:54146/251‑1 (MQ=255)
                                                                                                                                                                gTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATc                                                                                  <  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)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                                    
TGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCAC‑AGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCCTCTGGCA  >  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: ‑