Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,325,644 C→T S33N (AGC→AAC)  yjdN ← PF06983 family protein YjdN

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,325,6440CT98.3% 174.1 / NA 58S33N (AGC→AAC) yjdNPF06983 family protein YjdN
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base T (32/25);  minor base G (0/1);  total (32/26)
Fisher's exact test for biased strand distribution p-value = 4.48e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.85e-01

CATTTCGATTTTTCCGTTAGCGGCAAGATTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGCTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGG  >  NC_000913/4325400‑4325890
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
cATTTCGATTTTTCCGTTAGCGGCAAGATTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGAttt                                                                                                                                                                                                                                                  >  1:171707/1‑251 (MQ=255)
 aTTTCGATTTTGACGTTAGCGGCAAGATAGTCAAACCAGCGTTTTCCTTCTTCCAACTGTTGCGAATCGAGCACCAGCGGAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGACGCTTCCGGAAATGCGGACATTGACGGGGGCGATGGCGGTATCGGGAAGTTTCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGAGCGGATTTTGGAATTTCGCCGAAGTGGAtttt                                                                                                                                                                                                                                                 >  2:287912/1‑251 (MQ=255)
            tCCGTTAGCGGCAAGATTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGt                                                                                                                                                                                                                                       <  1:287194/250‑1 (MQ=255)
                  aGCGGCAAGATTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGTATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCcgcg                                                                                                                                                                                                                                 <  2:112738/250‑1 (MQ=255)
                     ggCAAGATTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTACGGAAGGGCAGTTCTCGGAGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGccc                                                                                                                                                                                                                              >  2:259252/1‑250 (MQ=255)
                                             tCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATc                                                                                                                                                                                                      <  1:351144/250‑1 (MQ=255)
                                                    gcGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCg                                                                                                                                                                                                <  2:186643/248‑1 (MQ=255)
                                                           gTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCCTGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGAACGTTGATAATAGGCAATCGCGTCGGAACAGtt                                                                                                                                                                                        >  2:476255/1‑250 (MQ=255)
                                                            ttGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGCTGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGTTTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTAc                                                                                                                                                                                      >  2:284507/1‑251 (MQ=255)
                                                                 aaTCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGc                                                                                                                                                                                  >  1:409244/1‑250 (MQ=255)
                                                                         aCCAGCGTAAAGCCGGAGTCGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCCATGCGAACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCCGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGCAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAagag                                                                                                                                                                           >  1:306321/1‑249 (MQ=255)
                                                                                 aaaGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGGGGTAGGGAc                                                                                                                                                                 >  2:102536/1‑251 (MQ=255)
                                                                                 aaaGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAAGCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGAc                                                                                                                                                                 <  1:112978/251‑1 (MQ=255)
                                                                                   aGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACtt                                                                                                                                                               <  1:121256/251‑1 (MQ=255)
                                                                                    gCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACtt                                                                                                                                                               >  1:389814/1‑250 (MQ=255)
                                                                                      cGGAGTAGCTGGCGTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAAc                                                                                                                                                            <  2:409244/251‑1 (MQ=255)
                                                                                                    tttCCTGACGGAATGGCATCGCTCATCATGATGTCGCTTCCGGCATTGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGAAGGGACTTAACGGCATGATGGTgcc                                                                                                                                              >  2:339242/1‑249 (MQ=255)
                                                                                                        cTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTc                                                                                                                                           >  1:138255/1‑250 (MQ=255)
                                                                                                                     atCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGCAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGACGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTCCGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAagag                                                                                                                                                                           >  2:345027/1‑205 (MQ=255)
                                                                                                                     atCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGTCAGTTCTCGGCGCTGTCCTGCGCTGATTGTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAagag                                                                                                                                                                           <  1:345027/205‑1 (MQ=255)
                                                                                                                       cGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGCACTCAGTTAAGAGAAGcc                                                                                                                             >  2:208031/1‑249 (MQ=255)
                                                                                                                       cGCTCATCATGATGGCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACTTCAGTTAAGAGAAGCCAg                                                                                                                           <  1:132467/251‑1 (MQ=255)
                                                                                                                          tcatcaTGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGACCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGCCAATCGCGTCGGCACAGGTACCGGCAAAAGAGAGGTAGGGCCTTCACGGCATGTTGGTGACATCATTTAAGAGACGCCAgat                                                                                                                         >  2:406632/1‑248 (MQ=255)
                                                                                                                              catGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGCTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAgtg                                                                                                                    >  1:387448/1‑251 (MQ=255)
                                                                                                                              catGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCAGGAAGGGCAGTTCTCGGCGCTGACCTGCGCTGATTTTGGCATTTCGCCGACGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCACAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGGTAAGAGAAGCCAGTTAAgtg                                                                                                                    >  1:332973/1‑251 (MQ=255)
                                                                                                                                  atgTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGATAAGCCAGTTAAGTGTAGt                                                                                                                >  2:125189/1‑251 (MQ=255)
                                                                                                                                  atgTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGt                                                                                                                >  1:58185/1‑251 (MQ=255)
                                                                                                                                       gCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGACGTTGATTGTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCCATCGCGTCGGAACAGTTACCGGCAAAAGAGCGGTAGGGACTTAACGGAATGATGGTGACCTCAGTTACGAGAAGCCAGTTAAGTGGAGTTCGa                                                                                                            >  2:195327/1‑250 (MQ=255)
                                                                                                                                       gCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGaa                                                                                                           <  1:244594/251‑1 (MQ=255)
                                                                                                                                       gCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGa                                                                                                            >  1:95848/1‑250 (MQ=255)
                                                                                                                                        cTTCCGGCAATACGCCCATTGGCGTGGGCGATGGCGGTCTCGGGAAACTGCGTTCCGGAAGGGCAGTTCTCGGCGCTGCCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCAGGTTAAGAGTAGTTCGaa                                                                                                           >  2:76050/1‑250 (MQ=255)
                                                                                                                                                aagtcgcACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAATGGCAGTTCTCGGCGGTGGCCTGCGCTGAGTTTGGCATTTCGCCGAAGTGGATTTTATAGAGGAGGTCCGCGCCCAACGTACGTTGATAATAGTCAATCGCGTCGGAACAGTTACCGTCAAAAGAGAGGTAGGGACTGAACGGCATGATGGTGACCTCAGTTCAGAGAAGCCAGTTAAGTGGAGTTCGAAATTAATACa                                                                                                  <  2:432326/247‑1 (MQ=255)
                                                                                                                                                aaTGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACa                                                                                                  <  2:385570/251‑1 (MQ=255)
                                                                                                                                                 aTGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAAATGCATTCCGGAAGCGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATCGGCAATCGCGTCGGAACACTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCAAGATGGTGACCTCAGTTAGGAGAAGCCAGTTAAGTGTAGTTCGAAGTTTATAcgg                                                                                                 >  1:259225/1‑249 (MQ=255)
                                                                                                                                                   gcgcACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGAt                                                                                               <  1:63150/251‑1 (MQ=255)
                                                                                                                                                     gcACATCGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTGCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCGCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATga                                                                                             <  1:246380/251‑1 (MQ=255)
                                                                                                                                                       acaTTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATgaga                                                                                           <  1:238196/251‑1 (MQ=255)
                                                                                                                                                         aTTGGCGTGGAAGATGGCAGTTTCGGGAAATTGAACTCCGGAAGGGACGATCTCGGCGATGTCATACGCTGATTTTGGCCTTTCGGCGAAGTTGATTTTATAGAGAAGTTCAGCGCCCAACATACGTTGGTAATAGGCAATTGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGGTAAGGGAAGGCAGATAAGTGTAGTTCGGAATTTATACAGATGAGAgg                                                                                         >  2:339452/1‑251 (MQ=255)
                                                                                                                                                               gTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCt                                                                                   <  2:219523/251‑1 (MQ=255)
                                                                                                                                                                tgagCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTg                                                                                  >  2:400904/4‑251 (MQ=255)
                                                                                                                                                                tGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTg                                                                                  >  1:73798/1‑251 (MQ=255)
                                                                                                                                                                 ggtcGATGGCGGTATCTGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTATGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATATGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGTCTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTg                                                                                  <  1:215926/247‑1 (MQ=255)
                                                                                                                                                                    cGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGATATTTATACAGGTGAGAGGCGGCCTGATa                                                                               >  2:119991/1‑250 (MQ=255)
                                                                                                                                                                     gATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTAGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTCACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTACACAGATGAGAGGCGGCCTGaca                                                                               >  2:163672/1‑247 (MQ=255)
                                                                                                                                                                        ggCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGAcg                                                                          >  1:5096/1‑251 (MQ=255)
                                                                                                                                                                        ggCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGAcg                                                                          <  1:284507/251‑1 (MQ=255)
                                                                                                                                                                            gTATCGGGAAATTGCATTCCGGACGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTAcc                                                                                                                                                                                     >  2:318411/1‑140 (MQ=255)
                                                                                                                                                                            gTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTAcc                                                                                                                                                                                     <  1:318411/140‑1 (MQ=255)
                                                                                                                                                                                             tCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCCGTTAAGTGTCGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGc                                                       >  2:124847/1‑249 (MQ=255)
                                                                                                                                                                                                              cGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGcc                                       <  1:343016/248‑1 (MQ=255)
                                                                                                                                                                                                              cGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGGTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGcc                                       <  2:389814/248‑1 (MQ=255)
                                                                                                                                                                                                                       ccgtcgcTGATTGTGGCAGTTCGCCGAAGTTGAGTTTATAGAGCAGTTCCGCGCCCAACGTACGGTGATAATAGGCAAGCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGGAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGGGCGGCGCATCAGGCAGGCGGCACTGGTGCCGGATGCGGCGTa                           <  2:234702/247‑1 (MQ=255)
                                                                                                                                                                                                                          gcgcTGATTGTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTGGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTTTTGCCGGAGGCGGCGTaaa                         <  1:383493/250‑1 (MQ=255)
                                                                                                                                                                                                                            gcTGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGc                      >  1:45191/1‑251 (MQ=255)
                                                                                                                                                                                                                              tGATTTTGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACATCAGTTAAGAGAAGCCAGTTAAGTGTCGTTCGAAATTTATACAGATACGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCCAGCAGTCGGCAATGTTGCCGGATGCGGCGTAAACGCCt                    >  2:293109/1‑251 (MQ=255)
                                                                                                                                                                                                                                    tGGCATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATcc               <  2:145494/250‑1 (MQ=255)
                                                                                                                                                                                                                                       cATTTCGCCGAAGTTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGGTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGGTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGcc           <  2:367737/251‑1 (MQ=255)
                                                                                                                                                                                                                                                gAAGTTGATTTTATAGAGCAGTCCCGCGCCCAACTTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGTGACTTTTCGTCATGCTGGTGCCCTCAGTTAAGAGAAGCCAGTTAATTGTATTTCGAAATTTATACAGATGACAGGCGTCCTGATCAGACGCGCGAGCGTCGCATCAGGCGGTCTGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAggg  <  1:311617/251‑1 (MQ=37)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                      
CATTTCGATTTTTCCGTTAGCGGCAAGATTGTCAAACCAGCGTTTTCCTTCTTCGACCTGTTGCGAATCGAGCACCAGCGTAAAGCCGGAGTAGCTGGCTTTTCCTGACGGCATGGCATCGCTCATCATGATGTCGCTTCCGGCAATGCGCACATTGGCGTGGGCGATGGCGGTATCGGGAAATTGCATTCCGGAAGGGCAGTTCTCGGCGCTGTCCTGCGCTGATTTTGGCATTTCGCCGAAGCTGATTTTATAGAGCAGTTCCGCGCCCAACGTACGTTGATAATAGGCAATCGCGTCGGAACAGTTACCGGCAAAAGAGAGGTAGGGACTTAACGGCATGATGGTGACCTCAGTTAAGAGAAGCCAGTTAAGTGTAGTTCGAAATTTATACAGATGAGAGGCGGCCTGATAAGACGCGCGAGCGTCGCATCAGGCAGTCGGCACTGTTGCCGGATGCGGCGTAAACGCCTTATCCGGCCTACGGAGGG  >  NC_000913/4325400‑4325890

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: