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,4330CA94.4% 89.1 / 0.0 37A319D (GCC→GAC) msrPprotein‑L‑methionine sulfoxide reductase catalytic subunit MsrP
Reads supporting (aligned to +/- strand):  ref base C (1/1);  new base A (18/16);  total (20/17)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.71e-02

TGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACC  >  NC_000913/2040187‑2040668
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            
tGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGCTTAAATCGATAGTCAGTATTAAGCTGACCCGCGCGCGTCCGCCAACCACCTGGAATCTGGAAGCGCCTGACGAATCCGGGTATTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGACACCGAACGTTTTATTGGTTCAGGCGGCATCCTCGATGTATAGCGCCAGCCAACGCGACTGTTTATTGGTTACgcagccacg                                                                                                                                                                                                                                          >  2:266258/1‑248 (MQ=255)
                 tGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGCTCACCCGCGCTGGTAACAGGCTACCGAACGATTTATTGGTTCAGGCGACATCCTCGATGTACAGCGCCAGCCAACGATCCTGTTTAAAGGTTACGCAGACCAAGTGGCATccacgtatcg                                                                                                                                                                                                                         >  2:2271/1‑242 (MQ=255)
                     aaTATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGc                                                                                                                                                                                                                     >  1:327161/1‑251 (MQ=255)
                                    ggATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGgag                                                                                                                                                                                                      <  1:163861/251‑1 (MQ=255)
                                      aTTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGgaga                                                                                                                                                                                                     >  1:229364/1‑250 (MQ=255)
                                          aaTCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGCCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCAAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTc                                                                                                                                                                                                >  1:337274/1‑251 (MQ=255)
                                             gattagcagttatcaatatcaCCCGCGAGCTCTCGCCAACCCCCTGGCATTTGGCAGCTCCTGACGAAGACGGGTTTTTCTCCAACGTTACTCCGTATGTTGTTCACCTTCGCTGGTCACAGGCTACCGAACAATTTTTTGGTTCAGGTGCCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTGAATGGTTACGTCGACCAGGTGGCATCACTGTCTCGTGGTCTGGATTTGCGGGAGAATTTCTaa                                                                                                                                                                                             <  2:245451/231‑1 (MQ=255)
                                                  gTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCg                                                                                                                                                                                        >  1:345643/1‑251 (MQ=255)
                                                     aGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCt                                                                                                                                                                                     >  1:16702/1‑251 (MQ=255)
                                                          tAAGCTGACCCGCGGGCGTCCGCCAACCACCGGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACa                                                                                                                                                                                 <  2:96432/250‑1 (MQ=255)
                                                             gCTGACCCGCGAGCGTCCGCCAACCTCCTGGAATCTGGCAGCGCCTGACGAATACGGTTTGTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCaa                                                                                                                                                                             <  2:95302/251‑1 (MQ=255)
                                                                       gAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTAAGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCAGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCAGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAACTGCGTCTGACAGCAAAACAGGTGAc                                                                                                                                                                   >  1:315264/1‑251 (MQ=255)
                                                                                                               gAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGCTGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAACACAGGTGACATGGCTGAAAGTTTGCCTGCAt                                                                                                                                             >  2:61540/1‑234 (MQ=255)
                                                                                                               gAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCAt                                                                                                                                             <  1:61540/234‑1 (MQ=255)
                                                                                                               gAATACGGTTTTTAAGCCAACGTTATTCCGTATGTTGATCACCCGCGCTGGTCACCGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCATCCAACGCTACTGTTTAATGGTTACGCCGACCAGATGGCATCGCTGTATCGTGGCCTGGCTTTGCGGGATAATTTCTAAATGCGTCTGACAGCCACCAAGGTGACATGGCTGAACGTTTGCCGGCATCTTACCGGCTTGTTGcc                                                                                                                            >  2:88499/1‑251 (MQ=255)
                                                                                                                            aCGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTctggctg                                                                                                               <  1:198495/251‑1 (MQ=255)
                                                                                                                              cctctaccgttctcgtttGTTGATCACCCGCGCTTGTCACAGGCTACCTAACGATGTAGTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGTATGTCTAAATGCTTCTGACAGCAAATCAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTGTCTCTGGCTGGt                                                                                                             <  2:205302/234‑1 (MQ=255)
                                                                                                                                           cGTATGTTGTTCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCTACCAGGTGTCATCGCTGTATCGTGTCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGCCATGGCTGAAAGTTTGCCTGCCTCTTGCCGGATTGTTGCCGTTTCTCTGGCTGATCTGCGCGatcaat                                                                                                >  1:150219/1‑251 (MQ=255)
                                                                                                                                           cGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGTGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGatcaat                                                                                                >  1:202320/1‑251 (MQ=255)
                                                                                                                                           cGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGt                                                                                                             >  1:205302/1‑238 (MQ=255)
                                                                                                                                                gTTGATCTCCCGGGCTGGTCCCAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCCGGAGGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCATGTGGCATCGTTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGTGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTGTGCCTGCATCTTGCCTGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACgg                                                                                           <  2:202320/251‑1 (MQ=255)
                                                                                                                                                   gATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACCCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTgg                                                                                        <  2:66723/251‑1 (MQ=255)
                                                                                                                                                      cACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACt                                                                                    >  1:133856/1‑251 (MQ=255)
                                                                                                                                                                    acaGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCCGGTGGCATCGCTGTATCGTGGCCTGGATATGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAAAAGGTGACATGGCTGAACGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTTTCTGGCTGGTCTGGGCGTTCACTCACGTGGGCATGGGTGCCGATCCgg                                                                      >  2:337165/1‑251 (MQ=255)
                                                                                                                                                                                 cGATGTATTGGTTCTGTCGGCATCCTCGATGTACATTGGCAGCCAACTCTACTGTTTAATGGTTACTCCGATCAGGTTGCATCGCTGTATCGTGGCCGGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCGGCATCTTGCCGGATTGTGGCCGTTTCTCTGGCTTGACTGGGTGATCAATCTCAGTGGACTGGGGGCCGATCCGTTGAAAGATATTCa                                                         <  2:229364/251‑1 (MQ=255)
                                                                                                                                                                                      tATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAACTGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCGCATCTTGCCGGATTGTTGCCGTTTCTTTGGCTGGTCTGGGCGATCAATCCAGGTGGACTGGGTGCCGATCCGGGGAAAGATATTCAGCCttt                                                   >  2:133898/1‑251 (MQ=255)
                                                                                                                                                                                          ggTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTAc                                                <  1:88499/251‑1 (MQ=255)
                                                                                                                                                                                               aGGCTGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTc                                           <  1:136442/251‑1 (MQ=255)
                                                                                                                                                                                                  cggcATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTgg                                                                                        >  1:107086/1‑204 (MQ=255)
                                                                                                                                                                                                  cggcATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTGTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTgg                                                                                        <  2:107086/204‑1 (MQ=255)
                                                                                                                                                                                                                            caACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGAcc              <  1:253428/251‑1 (MQ=255)
                                                                                                                                                                                                                              aCGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCtt            <  1:206643/251‑1 (MQ=255)
                                                                                                                                                                                                                               cGCTACTGTTTAATGGTTACGCAGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGAGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGGTATTCAGAATTTTACTGGTCGACCTGCACTGAACTTTTTGCTGGCGAACttt           >  2:280573/1‑250 (MQ=255)
                                                                                                                                                                                                                                  tACTGTTTAATGGTTACGCCGGCCAGGTGGCATCGCTGAATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCAGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGtt         >  1:241008/1‑250 (MQ=255)
                                                                                                                                                                                                                                    cTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCCCTGCACTGAAATTTTTGCTGGCGACCTTGTTAAt      >  2:284168/1‑251 (MQ=255)
                                                                                                                                                                                                                                         tAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCAcc  <  1:188980/250‑1 (MQ=255)
                                                                                                                                                                                                                                            tGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGCTTTGCGGGAGCATTTCTAAATGCGTCTGACAGCAAAACAGGTGCCATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGATGCCGTTTCTCTGGCTTGTCTGGGCGATCAATCCCGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCCTTTTAATGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTa        >  2:82000/1‑241 (MQ=255)
                                                                                                                                                                                                                                            tGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTa        <  1:82000/241‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            
TGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACC  >  NC_000913/2040187‑2040668

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