Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,040,289 C→T A271V (GCG→GTG)  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,2890CT98.1% 173.5 / NA 53A271V (GCG→GTG) msrPprotein‑L‑methionine sulfoxide reductase catalytic subunit MsrP
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (24/28);  minor base G (1/0);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 4.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

CAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCA  >  NC_000913/2040043‑2040529
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                
caggccagCAGGACCTCTTTATGGGCGGGGGGCTAAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCCCACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTTGCAGTGCCt                                                                                                                                                                                                                                              <  1:473986/251‑1 (MQ=255)
                         gggcgggtaaaattctcgttATGTCGAAGGTTTGCGTCTCGACGAACAAAGGCACCCGCCCACCCTTGTGCCTGTGGGTTTTTATGGCAAGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACg                                                                                                                                                                                                                     <  2:162271/233‑1 (MQ=255)
                                  aaaaTATCCTTATGTCGAAGGAGTGCGTCTCGAGGAAGCAATGCTCCCGCTCACGCTTATGACCGGAGGTTTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTGTACGCCAACGTTAATCCGt                                                                                                                                                                                                            <  1:194856/250‑1 (MQ=255)
                                        tcttttggccaagggATTGCGTCAGGACGAAGAAATGCACCCGCCAACTCTGATGACGGTAGTTTGTTATGGAATGGCGTTTCCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAAGCGATAGTCAGTATGAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACTAATTCGGTTTTTACGCCAACGTGAATCCGTATGtt                                                                                                                                                                                                       <  2:575299/237‑1 (MQ=255)
                                        tCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTg                                                                                                                                                                                                      >  1:241216/1‑251 (MQ=255)
                                        cccTTATGTCGAAGGATTGCCTCCCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGTTTTTTATGGCAAGTCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGtt                                                                                                                                                                                                       <  1:285990/249‑1 (MQ=255)
                                            ttggTCGAAGAATTGCGTCTCGACGAAACAATGCATCCGCCCACATTGAGGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATc                                                                                                                                                                                                   <  2:60487/247‑1 (MQ=255)
                                                 agAAGTATTGCGTCCGGTCGAAGCAATGCACCCGCTCACTTTTGTGGCGGTAGTTGTTTATGGAAGGGTGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCcg                                                                                                                                                                                              <  1:431202/249‑1 (MQ=255)
                                                  gAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTAAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCcgcg                                                                                                                                                                                            >  2:47276/1‑251 (MQ=255)
                                                         tgggtcccgacaagccaatacccccGCTCACAGGGATGCCTGTGGGTGTTTAGGGAAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTc                                                                                                                                                                                      <  2:403400/227‑1 (MQ=255)
                                                          gcgtcaggacgaAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAACTCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCCCCCGCGCTGGTcac                                                                                                                                                                                    >  2:75215/8‑251 (MQ=255)
                                                           cGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTcac                                                                                                                                                                                    >  1:223973/1‑250 (MQ=255)
                                                            gTCTCGACGAAGCAATGCCTCCGCCCACACTGATGACGGTAGGTGTTTATGGAAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAg                                                                                                                                                                                  <  1:436767/251‑1 (MQ=255)
                                                                 gacgaAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTAc                                                                                                                                                                             >  1:401247/1‑251 (MQ=255)
                                                                  acgaAGAAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTAcc                                                                                                                                                                            <  2:142178/251‑1 (MQ=255)
                                                                    gaAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCg                                                                                                                                                                           >  1:335381/1‑250 (MQ=255)
                                                                                  gCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGCAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTTCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGtt                                                                                                                                                            >  2:173360/1‑251 (MQ=255)
                                                                                                    aggggtTTATGGCAAGGAGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCAGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACAACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGATAATCCGCATGTTGATCACCCGCCCTGGCCACAGGCTACCGAAAGATCTATTGGTTCAGGCGGCCTCCTCGATg                                                                                                                                          >  2:350760/5‑251 (MQ=255)
                                                                                                           tagggCAGGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAgc                                                                                                                                    <  1:475462/247‑1 (MQ=255)
                                                                                                                aaaGGCTTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCgccagc                                                                                                                              <  2:420034/250‑1 (MQ=255)
                                                                                                                   ggCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGAATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGCTCAGGCGGCATCCTCGATGTACAGCgccagcca                                                                                                                            >  1:366837/1‑250 (MQ=255)
                                                                                                                       ttACCGACACAAAATGGCGCGCCGGCGCGACTGATTGTGCAGTGGAAATATGGGTTTAAAGGGATTAAATCGATAGTCCGTATTAAGATGCCCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTAGGTTGATCACCAGCGCTGGTCACCGGCTACCAAACGATTAATTGGTTCGGTCGGCATCATCGATGTCCAgcgcacgcccacgct                                                                                                                       >  2:269680/1‑240 (MQ=255)
                                                                                                                              cacaAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGttt                                                                                                                <  1:420443/251‑1 (MQ=255)
                                                                                                                                        cgcgcCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTAcg                                                                                                      <  1:484241/251‑1 (MQ=255)
                                                                                                                                        cgcgcCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTAcg                                                                                                      <  1:62795/251‑1 (MQ=255)
                                                                                                                                             cGGTGCTACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCGGTCCCGCGAGCCTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTGTTTACGCCAACGTTAATGCGTATGTTGATCACCCGCGCTGGTCATAGGCTACCGAACGATTTAGTGGTGCAGGCGGCATCCTCGATGTACAGCGCCAGCAAACGCGACGGTCTAATGGTTAcgccgac                                                                                                 <  2:375633/251‑3 (MQ=255)
                                                                                                                                                 gCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAg                                                                                              <  1:79615/250‑1 (MQ=255)
                                                                                                                                                               cGTGGAAATATGGCTTTAAAGGGATGAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTTCGCTAACCGCCTGGAATCTGGCAGTGCCTGACGAATACGGGTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGACAGCCAACGCTACTGTTTAATGGTTACGCCGATCAGGTGGCATCGCTGTAt                                                                               <  1:433336/251‑1 (MQ=255)
                                                                                                                                                                gTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCTGCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTTGTCACAGGCTACCGAACGATTTATTGGTGCAGTCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGTCGACCAGGTGGCATCGCTGTAt                                                                               <  2:524104/250‑1 (MQ=255)
                                                                                                                                                                gTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTAt                                                                               >  1:524104/1‑250 (MQ=255)
                                                                                                                                                                gTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATc                                                                              >  2:21682/1‑251 (MQ=255)
                                                                                                                                                                       tatGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCACCAACCAGGTGGAATCTGGCAAGGCATGACGAATACGATTTTTACGCCAACGTTAATCCGTAAGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCt                                                                       >  2:283286/1‑251 (MQ=255)
                                                                                                                                                                        atGGCTTTAAAGGGATTAAATCGATAGTGAGTATTAAGCTGACCCGCGAGCGTCCGCCACACACCTGGAATATGGCAGTGCCTGACGACTACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTg                                                                      >  2:584581/1‑251 (MQ=255)
                                                                                                                                                                           gCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTggct                                                                   >  2:465930/1‑249 (MQ=255)
                                                                                                                                                                               tAAAGAGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGAGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGCTCCCCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTCCAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGc                                                               >  2:353691/1‑251 (MQ=255)
                                                                                                                                                                                    ggATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGCTGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCggg                                                            >  2:554400/1‑249 (MQ=255)
                                                                                                                                                                                      aTTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGATTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCATGGATTTGCGGgaga                                                         >  1:273092/1‑250 (MQ=255)
                                                                                                                                                                                         aaaTCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAAtt                                                      <  1:47276/250‑1 (MQ=255)
                                                                                                                                                                                            tCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCt                                                   >  1:424148/1‑250 (MQ=255)
                                                                                                                                                                                             cGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGAGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTa                                                  >  1:235666/1‑250 (MQ=255)
                                                                                                                                                                                              gATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACTCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTaaa                                                <  2:273092/251‑1 (MQ=255)
                                                                                                                                                                                                     aGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCt                                         <  1:173360/251‑1 (MQ=255)
                                                                                                                                                                                                            tgCTGACTCGCGAGCGTCCCCCAACAACCTGGAATCTGGCAGTGCCTGACGAATACGGCTTTTACGCCAACGTTCATCCGTATGTTGATCACCCGCGCTGGTCCCAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGc                                   >  2:92971/2‑250 (MQ=255)
                                                                                                                                                                                                             gCTGACCCGCGAGCGTCCGGCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCCCCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCaa                                 <  1:89403/251‑1 (MQ=255)
                                                                                                                                                                                                              cTGACCCGAGAGCGTCAGCCAACCACCTGGAATCTGGCCGTGCCTGCCGAATACGGTTTTTACGCCCACGTTCATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCATCGCTCCTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCaa                                 >  2:500648/1‑250 (MQ=255)
                                                                                                                                                                                                                   ccGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAgg                           <  2:36081/251‑1 (MQ=255)
                                                                                                                                                                                                                        aGCGTCCGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTTGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTTACAGCAAAACAGGTGACa                      <  2:510462/251‑1 (MQ=255)
                                                                                                                                                                                                                              cGCCAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTg                <  1:403199/251‑1 (MQ=255)
                                                                                                                                                                                                                                ccAACCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGa               >  1:2430/1‑250 (MQ=255)
                                                                                                                                                                                                                                  aaCCACCTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGAGACATGGCTGaa              >  1:573507/1‑249 (MQ=255)
                                                                                                                                                                                                                                   accaccTGGAATCTTGCAGTGTCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGt           <  2:511750/251‑1 (MQ=255)
                                                                                                                                                                                                                                    ccaccTGGAATCTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGtt          >  1:325627/1‑251 (MQ=255)
                                                                                                                                                                                                                                            aagcTGGCAGTGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCa  <  1:22223/248‑1 (MQ=255)
                                                                                                                                                                                                                                                      |                                                                                                                                                                                                                                                
CAGGCCAGCAGGACCGCTTTATCGGCGGTGGGCTGAAATATCCTTATGTCGAAGGATTGCGTCTCGACGAAGCAATGCATCCGCTCACACTGATGACCGTAGGTGTTTATGGCAAGGCGTTACCGCCACAAAATGGCGCGCCGGTGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCA  >  NC_000913/2040043‑2040529

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