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,4330CA97.9% 131.6 / NA 49A319D (GCC→GAC) msrPprotein‑L‑methionine sulfoxide reductase catalytic subunit MsrP
Reads supporting (aligned to +/- strand):  ref base C (1/0);  new base A (21/26);  total (22/27)
Fisher's exact test for biased strand distribution p-value = 4.58e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCCT  >  NC_000913/2040194‑2040671
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                
gATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGc                                                                                                                                                                                                                                         <  1:143376/249‑1 (MQ=255)
    gtgCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGATGGCATCGCt                                                                                                                                                                                                                                   >  1:288894/1‑251 (MQ=255)
      gCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGt                                                                                                                                                                                                                                 >  1:358732/1‑251 (MQ=255)
          gtgcaaTATGGTTTTAAAGGGAATAAATTGATATACAGTTTTAAGCTTTACCGCGAGCGTCCGCCAAACACCTGGAATCTGGCAGCGCCTGACGATTACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGCACGATTTATTGGTTCAGGCGGCATTCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCg                                                                                                                                                                                                                             <  2:172650/247‑1 (MQ=255)
                  tGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTgg                                                                                                                                                                                                                     >  1:387494/1‑251 (MQ=255)
                    gCTTTGAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGACCGATTTATTGGTTCAGGCGGCCTCCTCGATGTACAGCGACAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGGAGCATCGCTGTATCGTGGCCTGGAt                                                                                                                                                                                                                   >  2:100375/1‑251 (MQ=255)
                       ttcaaGGGTTTACATCGATAGTCAGCATTACGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGACAGCGCCAGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGCTCACCCGCGCTGGTCACAGGCTACCGCACGATTTATTGGTTCAGGCAGCATCCTCGATGTACAGCGCCAGCCAACGCAACTGTTTAATGGTTACGCTGACCAGGTGGCAACGCTGTATCGTGGCCTGGATTTg                                                                                                                                                                                                                >  2:169654/4‑251 (MQ=255)
                       ttAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGAttt                                                                                                                                                                                                                 <  1:298769/250‑1 (MQ=255)
                                 tAAATCGATAGTCAGTATTAAGCTGACCCGCGTGCGTCCCCCAACCACATGCAATCTGGCAGCGCCTGCCAAATACGGTTTTTACGCCAACGTTACTACGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAAAGATTTATTGGTTCAGACGGCCTCCTCGATGTACAGCGCCACCCCACGCTACTGCTTAATGGTTACGCAGCACAGGTGGCATCGCGGTATCATGGCATTGCTTTGCGGGAGAAtt                                                                                                                                                                                                      >  1:271482/1‑251 (MQ=255)
                                  aaaTCGATAGTCAGTATTAAGCTGACCCGCGCGCGTCCGCAAACCACCTGGAATCTGGCAGCGCCTGACGCATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGTCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAAttt                                                                                                                                                                                                     >  1:166789/1‑251 (MQ=255)
                                                     aGCTGACCCCCGAGGGTCCACCAACCACCTGGAATCGGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGAGTTATTGGTTCAGGCGGCATCCGCGATGTACAGCGCCAGCCAACGCGACTGTTTAATGGTTACGCCGACCAGGTGGCAGCGCCGGAGCGTGGCCTGGATTTGCGTGAGAATTGCTAAATGCGTCTGACAGCa                                                                                                                                                                                 <  2:307930/251‑1 (MQ=255)
                                                         gACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAAc                                                                                                                                                                             >  1:16160/1‑251 (MQ=255)
                                                                              caccTGGAATCTGGCAGCGCCTGACGAATACGGTTTGTACGCCATCGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGttt                                                                                                                                                        <  2:315532/251‑1 (MQ=255)
                                                                                          ggCAGCGCCTTACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTGTGCCTGCATCTTg                                                                                                                                            <  2:270571/251‑1 (MQ=255)
                                                                                          ggCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCAt                                                                                                                                                <  1:447884/247‑1 (MQ=255)
                                                                                           gCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGc                                                                                                                                           >  1:118691/1‑251 (MQ=255)
                                                                                            cAGCGCCTGACGAATACGGTTTTTTCTCCAACGTTAATCCGTATGTTGATCACCCGCGCTTGTCACAGGCTACCGAACGATTTATTGGTGCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGGATCGCTGTATCGTGGCCTGGAGTTGCGGGAGAATTACTAAATGCGTCTTACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGcc                                                                                                                                          <  2:124440/251‑1 (MQ=255)
                                                                                            cAGCGCCTGACGAATACGGTTTTTACGCCATCGTTAAATCGTATGTTGATCCCCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGGACAGCGACAGCCATCGCTTCTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGcc                                                                                                                                          <  2:66322/251‑1 (MQ=255)
                                                                                                    gACGAATACAGGTCTTACGCCAACGTTCATCTGAATGTAGCTCACCCGCGATGGTCACAGGCAACCGAACGATTTATTGGTTCCGGCGGCATCATCGATGTACAGCGCCAGCCAACGCTACTGTTTAATCGTTCCGCCGACCAGCTGGCCTAGCTGTCTCGTGGGCTGGCTTTGCGGGAGAATTTCTAAATGCGTCTGAAAGCCAAAACGGTGACATAGCCGACAGTTTGCGTGCATTTTGCCGGAttgct                                                                                                                                  >  2:306711/1‑249 (MQ=255)
                                                                                                        aagaCGGTCTTTACGCCTACGTTAATCCGTATGTTGATCACCCGCGCTTGTCACAGGCTACCGAACGATTTAGTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACTCTACTGTGTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCg                                                                                                                              <  2:152473/248‑1 (MQ=255)
                                                                                                         aTACGGTTTTTTCGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCCAAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGAGTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGt                                                                                                                             <  2:165360/251‑1 (MQ=255)
                                                                                                                ttttACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTctct                                                                                                                       <  1:432365/250‑1 (MQ=255)
                                                                                                                    aagcccaCGTTAATCCGTATGTTGATCACCCGCGCTCGTGACAGGCTACCGAACTAGTGATTGGTTTAGGCTGCATTCGCGATCGACCGCGCTAGCTAACGCTACTGTTTAATGTTTACTCCGACCATGTTGCATCGCTGTATCGTGGCCTGGATATGCTGGATAAGTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCTTTTTCtctgtctg                                                                                                                  <  2:164484/245‑1 (MQ=255)
                                                                                                                                 tCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGACGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATCCGTCTGCCAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGCCTGGTCTGGGCGatcc                                                                                                     >  2:391041/1‑250 (MQ=255)
                                                                                                                                       tGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCtt                                                                                                                                             >  1:439363/1‑205 (MQ=255)
                                                                                                                                       tGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCtt                                                                                                                                             <  2:439363/205‑1 (MQ=255)
                                                                                                                                             tCACCCGCGCTAGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATACTCGATGTACAGAGCCCGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATGGCTGTATCGTGACCTGGATTTGCGGGAGAATTTATAAATGCGTCTGACAGAAAAACAGGGGACATGGCTGCAAGTTTGGCTGCATCTTGCCGGCTTGTTGCCGTTTCTCTGGCGGGTCTGGGCGatcaatcccgaggggt                                                                                         >  2:416644/1‑242 (MQ=255)
                                                                                                                                              cACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACt                                                                                        >  2:8358/1‑251 (MQ=255)
                                                                                                                                                 ccGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCACAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTATCTGGCTGGTCTGGGCGATCAATCACGGTGGACTggg                                                                                     >  1:436761/1‑251 (MQ=255)
                                                                                                                                                        ggTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGAt                                                                              <  1:328395/251‑1 (MQ=255)
                                                                                                                                                         gTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATc                                                                             <  1:110849/251‑1 (MQ=255)
                                                                                                                                                                      acaCGATTTATTGGTTCAAGCGGCTTCCTCGATGCACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTTCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGGTTGTTGCCGTTTCTCTGGCTGGTCTGGGCGTTCAATCACGGGGGACTgggctccggaccggggaactatat                                                                >  2:40237/3‑230 (MQ=255)
                                                                                                                                                                        aCGATTTATTGGTTCAGGCGGCATCCGCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTc                                                              <  2:393623/251‑1 (MQ=255)
                                                                                                                                                                               aTTGGTTCAGGCGGCATCCTCGATGGACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCTGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGTGTGCCGATCCGGTGAAAGATATTCAGCAtt                                                        <  1:416644/250‑1 (MQ=255)
                                                                                                                                                                                        ggtggcAGCCTCGATGGTCAGGGCCAGCCAACGCTTCTGTTTAATGGTTACGCCGACCATGTGGCATCGGTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGTCTGCATCTTGCCGGATTGTTGCCGTGTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTgggct                                              <  2:332013/248‑4 (MQ=255)
                                                                                                                                                                                               tCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCgcactgc                                        >  1:447449/1‑250 (MQ=255)
                                                                                                                                                                                                ccTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCgcactgcac                                      >  1:401399/1‑251 (MQ=255)
                                                                                                                                                                                                   cGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAAGGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCgcactgcact                                     <  1:94494/249‑1 (MQ=255)
                                                                                                                                                                                                       gTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAAt                                <  1:109532/250‑1 (MQ=255)
                                                                                                                                                                                                          cAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAAttttt                            <  1:217224/251‑1 (MQ=255)
                                                                                                                                                                                                                   ccaACGCTACTGTTTAATGGTTACGCCGGCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGAc                   <  2:16160/251‑1 (MQ=255)
                                                                                                                                                                                                                     aaCGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCt                >  1:179909/1‑251 (MQ=255)
                                                                                                                                                                                                                          tACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTg                                                                                            >  1:59890/1‑171 (MQ=255)
                                                                                                                                                                                                                          tACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTg                                                                                            <  2:59890/171‑1 (MQ=255)
                                                                                                                                                                                                                              gTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAAt         <  1:293177/249‑1 (MQ=255)
                                                                                                                                                                                                                              gTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATc        <  1:164125/250‑1 (MQ=255)
                                                                                                                                                                                                                               tttAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCCATGCACTGATAttttt                            >  2:323106/1‑230 (MQ=255)
                                                                                                                                                                                                                               tttAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAAttttt                            <  1:323106/230‑1 (MQ=255)
                                                                                                                                                                                                                                   aTGGTTACGCCGACAAGGTGGAATCGCTGTATAGAGGCCTAGATTTGCGGGAGAATTTCTAAAAGCGTCTGACAGCAATACAGATGACATGGCTGAAAGTTTGCCTGCATCTTGACGGATTGTTGCCGTTTCTCTGGCTGGTCTTGTCGATCACTCACGGTTGACTGTGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCAATACCCTGAATTTTGTGTGGGGGCACTTGTTAATCACCCct  >  2:280192/1‑251 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                
GATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGCAAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGACCTTGTTAATCACCCCT  >  NC_000913/2040194‑2040671

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