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

GATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGAC‑CTTGTTAATCACCCCT  >  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)
                                                     aGCTGACCCCCGAGGGTCCACCAACCACCTGGAATCGGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGAGTTATTGGTTCAGGCGGCATCCGCGATGTACAGCGCCAGCCAACGCGACTGTTTAATGGTTACGCCGACCAGGTGGCAGCGCCGGAGCGTGGCCTGGATTTGCGTGAGAATTGCTAAATGCGTCTGACAGC‑a                                                                                                                                                                                 <  2:307930/251‑1 (MQ=255)
                                                         gACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAAc                                                                                                                                                                             >  1:16160/1‑251 (MQ=255)
                                                                              caccTGGAATCTGGCAGCGCCTGACGAATACGGTTTGTACGCCATCGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGttt                                                                                                                                                        <  2:315532/251‑1 (MQ=255)
                                                                                          ggCAGCGCCTTACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTGTGCCTGCATCTTg                                                                                                                                            <  2:270571/251‑1 (MQ=255)
                                                                                          ggCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCAt                                                                                                                                                <  1:447884/247‑1 (MQ=255)
                                                                                           gCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGc                                                                                                                                           >  1:118691/1‑251 (MQ=255)
                                                                                            cAGCGCCTGACGAATACGGTTTTTTCTCCAACGTTAATCCGTATGTTGATCACCCGCGCTTGTCACAGGCTACCGAACGATTTATTGGTGCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGGATCGCTGTATCGTGGCCTGGAGTTGCGGGAGAATTACTAAATGCGTCTTACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGcc                                                                                                                                          <  2:124440/251‑1 (MQ=255)
                                                                                            cAGCGCCTGACGAATACGGTTTTTACGCCATCGTTAAATCGTATGTTGATCCCCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGGACAGCGACAGCCATCGCTTCTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGcc                                                                                                                                          <  2:66322/251‑1 (MQ=255)
                                                                                                    gACGAATACAGGTCTTACGCCAACGTTCATCTGAATGTAGCTCACCCGCGATGGTCACAGGCAACCGAACGATTTATTGGTTCCGGCGGCATCATCGATGTACAGCGCCAGCCAACGCTACTGTTTAATCGTTCCGCCGACCAGCTGGCCTAGCTGTCTCGTGGGCTGGCTTTGCGGGAGAATTTCTAAATGCGTCTGAAAGCCAAAAC‑GGTGACATAGCCGACAGTTTGCGTGCATTTTGCCGGAttgct                                                                                                                                  >  2:306711/1‑249 (MQ=255)
                                                                                                        aagaCGGTCTTTACGCCTACGTTAATCCGTATGTTGATCACCCGCGCTTGTCACAGGCTACCGAACGATTTAGTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACTCTACTGTGTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCg                                                                                                                              <  2:152473/248‑1 (MQ=255)
                                                                                                         aTACGGTTTTTTCGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCCAAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGAGTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGt                                                                                                                             <  2:165360/251‑1 (MQ=255)
                                                                                                                ttttACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTctct                                                                                                                       <  1:432365/250‑1 (MQ=255)
                                                                                                                    aagcccaCGTTAATCCGTATGTTGATCACCCGCGCTCGTGACAGGCTACCGAACTAGTGATTGGTTTAGGCTGCATTCGCGATCGACCGCGCTAGCTAACGCTACTGTTTAATGTTTACTCCGACCATGTTGCATCGCTGTATCGTGGCCTGGATATGCTGGATAAGTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCTTTTTCtctgtctg                                                                                                                  <  2:164484/245‑1 (MQ=255)
                                                                                                                                 tCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGACGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATCCGTCTGCCAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGCCTGGTCTGGGCGatcc                                                                                                     >  2:391041/1‑250 (MQ=255)
                                                                                                                                       tGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCtt                                                                                                                                             >  1:439363/1‑205 (MQ=255)
                                                                                                                                       tGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCtt                                                                                                                                             <  2:439363/205‑1 (MQ=255)
                                                                                                                                             tCACCCGCGCTAGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATACTCGATGTACAGAGCCCGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATGGCTGTATCGTGACCTGGATTTGCGGGAGAATTTATAAATGCGTCTGACAGA‑AAAACAGGGGACATGGCTGCAAGTTTGGCTGCATCTTGCCGGCTTGTTGCCGTTTCTCTGGCGGGTCTGGGCGatcaatcccgaggggt                                                                                         >  2:416644/1‑242 (MQ=255)
                                                                                                                                              cACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACt                                                                                        >  2:8358/1‑251 (MQ=255)
                                                                                                                                                 ccGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑ACAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTATCTGGCTGGTCTGGGCGATCAATCACGGTGGACTggg                                                                                     >  1:436761/1‑251 (MQ=255)
                                                                                                                                                        ggTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGAt                                                                              <  1:328395/251‑1 (MQ=255)
                                                                                                                                                         gTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATc                                                                             <  1:110849/251‑1 (MQ=255)
                                                                                                                                                                      acaCGATTTATTGGTTCAAGCGGCTTCCTCGATGCACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTTCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGGTTGTTGCCGTTTCTCTGGCTGGTCTGGGCGTTCAATCACGGGGGACTgggctccggaccggggaactatat                                                                >  2:40237/3‑230 (MQ=255)
                                                                                                                                                                        aCGATTTATTGGTTCAGGCGGCATCCGCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTc                                                              <  2:393623/251‑1 (MQ=255)
                                                                                                                                                                               aTTGGTTCAGGCGGCATCCTCGATGGACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCTGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGTGTGCCGATCCGGTGAAAGATATTCAGCAtt                                                        <  1:416644/250‑1 (MQ=255)
                                                                                                                                                                                        ggtggcAGCCTCGATGGTCAGGGCCAGCCAACGCTTCTGTTTAATGGTTACGCCGACCATGTGGCATCGGTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGTCTGCATCTTGCCGGATTGTTGCCGTGTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTgggct                                              <  2:332013/248‑4 (MQ=255)
                                                                                                                                                                                               tCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCgcactgc                                        >  1:447449/1‑250 (MQ=255)
                                                                                                                                                                                                ccTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCgcactgcac                                      >  1:401399/1‑251 (MQ=255)
                                                                                                                                                                                                   cGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAAGGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCgcactgcact                                     <  1:94494/249‑1 (MQ=255)
                                                                                                                                                                                                       gTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAAt                                <  1:109532/250‑1 (MQ=255)
                                                                                                                                                                                                          cAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAAttttt                            <  1:217224/251‑1 (MQ=255)
                                                                                                                                                                                                                   ccaACGCTACTGTTTAATGGTTACGCCGGCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGAc                   <  2:16160/251‑1 (MQ=255)
                                                                                                                                                                                                                     aaCGCTACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGAC‑Ct                >  1:179909/1‑251 (MQ=255)
                                                                                                                                                                                                                          tACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTg                                                                                            >  1:59890/1‑171 (MQ=255)
                                                                                                                                                                                                                          tACTGTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTg                                                                                            <  2:59890/171‑1 (MQ=255)
                                                                                                                                                                                                                              gTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGAC‑CTTGTTAAt         <  1:293177/249‑1 (MQ=255)
                                                                                                                                                                                                                              gTTTAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGAC‑CTTGTTAATc        <  1:164125/250‑1 (MQ=255)
                                                                                                                                                                                                                               tttAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCCATGCACTGATAttttt                            >  2:323106/1‑230 (MQ=255)
                                                                                                                                                                                                                               tttAATGGTTACGCCGACCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAAttttt                            <  1:323106/230‑1 (MQ=255)
                                                                                                                                                                                                                                   aTGGTTACGCCGACAAGGTGGAATCGCTGTATAGAGGCCTAGATTTGCGGGAGAATTTCTAAAAGCGTCTGACAGC‑AATACAGATGACATGGCTGAAAGTTTGCCTGCATCTTGACGGATTGTTGCCGTTTCTCTGGCTGGTCTTGTCGATCACTCACGGTTGACTGTGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCAATACCCTGAATTTTGTG‑TGGGGGCACTTGTTAATCACCCct  >  2:280192/1‑251 (MQ=255)
                                                                                                                                                                                                                                               |                                                                                                                                                                                                                                                
GATTGTGCCGTGGAAATATGGCTTTAAAGGGATTAAATCGATAGTCAGTATTAAGCTGACCCGCGAGCGTCCGCCAACCACCTGGAATCTGGCAGCGCCTGACGAATACGGTTTTTACGCCAACGTTAATCCGTATGTTGATCACCCGCGCTGGTCACAGGCTACCGAACGATTTATTGGTTCAGGCGGCATCCTCGATGTACAGCGCCAGCCAACGCTACTGTTTAATGGTTACGCCGCCCAGGTGGCATCGCTGTATCGTGGCCTGGATTTGCGGGAGAATTTCTAAATGCGTCTGACAGC‑AAAACAGGTGACATGGCTGAAAGTTTGCCTGCATCTTGCCGGATTGTTGCCGTTTCTCTGGCTGGTCTGGGCGATCAATCACGGTGGACTGGGTGCCGATCCGGTGAAAGATATTCAGCATTTTACTGGTCGCACTGCACTGAAATTTTTGCTGGCGAC‑CTTGTTAATCACCCCT  >  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: ‑