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

TGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTC‑CGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGAC‑TGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACC  >  NC_000913/2040187‑2040668
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            
tGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGCTTAAATCGATAGTCAGTATTAAGCTGACCCGCGCGCGTC‑CGCCAACCACCTGGAATCTGGAAGCGCCTGACGAATCCGGGTATTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGACACCGAACGTTTTATTGGTTCAGGCGGCATCCTCGATGTATAGCGCCAGCCAACGCGACTGTTTATTGGTTACgcagccacg                                                                                                                                                                                                                                          >  2:266258/1‑248 (MQ=255)
                 tGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTC‑CGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGCTCACCCGCGCTGGTAACAGGCTACCGAACGATTTATTGGTTCAGGCGACATCCTCGATGTACAGCGCCAGCCAACGATCCTGTTTAAAGGTTACGCAGACCAAGTGGCATccacgtatcg                                                                                                                                                                                                                         >  2:2271/1‑242 (MQ=255)
                     aaTATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTC‑CGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGc                                                                                                                                                                                                                     >  1:327161/1‑251 (MQ=255)
                                    ggATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTC‑CGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGgag                                                                                                                                                                                                      <  1:163861/251‑1 (MQ=255)
                                      aTTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTC‑CGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGgaga                                                                                                                                                                                                     >  1:229364/1‑250 (MQ=255)
                                          aaTCGATAGTCAGTATTAAGCTGACCCGCGAGCGTC‑CGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGCCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCAAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTc                                                                                                                                                                                                >  1:337274/1‑251 (MQ=255)
                                             gattagcagttatcaatatcaCCCGCGAGC‑TCTCGCCAACCCCCTGGCATTTGGCAGCTCCTGACGAAGACGGGTTTTTCTCCAACGTTACTCCGTATGTTGTTCACCTTCGCTGGTCACAGGCTACCGAACAATTTTTTGGTTCAGGTGCCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTGAATGGTTACGTCGACCAGGTGGCATCACTGTCTCGTGGTCTGGATTTGCGGGAGAATTTCTaa                                                                                                                                                                                             <  2:245451/231‑1 (MQ=255)
                                                  gTCAGTATTAAGCTGACCCGCGAGCGTC‑CGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCg                                                                                                                                                                                        >  1:345643/1‑251 (MQ=255)
                                                     aGTATTAAGCTGACCCGCGAGCGTC‑CGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCt                                                                                                                                                                                     >  1:16702/1‑251 (MQ=255)
                                                          tAAGCTGACCCGCGGGCGTC‑CGCCAACCACCGGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACa                                                                                                                                                                                 <  2:96432/250‑1 (MQ=255)
                                                             gCTGACCCGCGAGCGTC‑CGCCAACCTCCTGGAATCTGGCAGCGCCTGACGAATACGGTTTGTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCaa                                                                                                                                                                             <  2:95302/251‑1 (MQ=255)
                                                                       gAGCGTC‑CGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTAAGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCAGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCAGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAACTGCGTCTGACAGCAAAACAGGTGAc                                                                                                                                                                   >  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)
                                                                                                                                                      cACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGAC‑t                                                                                    >  1:133856/1‑251 (MQ=255)
                                                                                                                                                                    acaGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCCGGTGGCATCGCTGTATCGTGGCCTGGATATGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAAAAGGTGACATGGCTGAACGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTTTCTGGCTGGTCTGGGCGTTCACTCACG‑TGGGCATGGGTGCCGATCCgg                                                                      >  2:337165/1‑251 (MQ=255)
                                                                                                                                                                                 cGATGTATTGGTTCTGTCGGCATCCTCGATGTACATTGGCAGCCAACTCTACTGTTTAATGGTTACTCCGATCAGGTTGCATCGCTGTATCGTGGCCGGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCGGCATCTTGCCGGATTGTGGCCGTTTCTCTGGCTTGACTGGGTGATCAATCTCAGTGGAC‑TGGGGGCCGATCCGTTGAAAGATATTCa                                                         <  2:229364/251‑1 (MQ=255)
                                                                                                                                                                                      tATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAACTGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCC‑GCATCTTGCCGGATTGTTGCCGTTTCTTTGGCTGGTCTGGGCGATCAATCCAGGTGGAC‑TGGGTGCCGATCCGGGGAAAGATATTCAGCCttt                                                   >  2:133898/1‑251 (MQ=255)
                                                                                                                                                                                          ggTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGAC‑TGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTAc                                                <  1:88499/251‑1 (MQ=255)
                                                                                                                                                                                               aGGCTGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGAC‑TGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTc                                           <  1:136442/251‑1 (MQ=255)
                                                                                                                                                                                                  cggcATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTgg                                                                                        >  1:107086/1‑204 (MQ=255)
                                                                                                                                                                                                  cggcATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTGTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTgg                                                                                        <  2:107086/204‑1 (MQ=255)
                                                                                                                                                                                                                            caACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGAC‑TGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGAcc              <  1:253428/251‑1 (MQ=255)
                                                                                                                                                                                                                              aCGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGAC‑TGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCtt            <  1:206643/251‑1 (MQ=255)
                                                                                                                                                                                                                               cGCTACTGTTTAATGGTTACGCAGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGAGCGATCAATCACGGTGGAC‑TGGGTGCCGATCCGGTGAAAGGTATTCAGAATTTTACTGGTCGACCTGCACTGAACTTTTTGCTGGCGAACttt           >  2:280573/1‑250 (MQ=255)
                                                                                                                                                                                                                                  tACTGTTTAATGGTTACGCCGGCCAGGTGGCATCGCTGAATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCAGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGAC‑TGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGtt         >  1:241008/1‑250 (MQ=255)
                                                                                                                                                                                                                                    cTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGAC‑TGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCCCTGCACTGAAATTTTTGCTGGCGACCTTGTTAAt      >  2:284168/1‑251 (MQ=255)
                                                                                                                                                                                                                                         tAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGAC‑TGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCAcc  <  1:188980/250‑1 (MQ=255)
                                                                                                                                                                                                                                            tGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGCTTTGCGGGAGCATTTCTAAATGCGTCTGACAGCAAAACAGGTGCCATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGATGCCGTTTCTCTGGCTTGTCTGGGCGATCAATCCCGGTGGAC‑TGGGTGCCGATCCGGTGAAAGATATTCAGCCTTTTAATGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTa        >  2:82000/1‑241 (MQ=255)
                                                                                                                                                                                                                                            tGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGAC‑TGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTa        <  1:82000/241‑1 (MQ=255)
                                                                                                                                                                                                                                                       |                                                                                                                                                                                                                                            
TGCGACTGATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTC‑CGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGAC‑TGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACC  >  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: ‑