Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,875,269 T→G 18.5% T51P (ACC→CCC)  cysN ← sulfate adenylyltransferase, subunit 1

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,875,2690TG18.5% 119.2 / 2.5 63T51P (ACC→CCC) cysNsulfate adenylyltransferase, subunit 1
Reads supporting (aligned to +/- strand):  ref base T (16/32);  new base G (13/2);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 3.49e-04
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CATATTGCGGGTGTACTGCTCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCA  >  NC_000913/2875048‑2875496
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                   
cATATTGCGGGTGTACTGCTCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTg                                                                                                                                                                                                                               <  2:248806/228‑1 (MQ=255)
        gggtgttacgCTCGGGCCCTGGGGTGTCGGGGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAga                                                                                                                                                                                                                        <  1:150681/218‑1 (MQ=255)
        gggTGTATCGCCCGGGCACTGGGGTGTCGGCTGTAATAAATTTACGCTTCGCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTgcagca                                                                                                                                                                                                                          <  1:290159/225‑1 (MQ=255)
          gtgtACTGCTCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCCGCAgacga                                                                                                                                                                                                                     >  1:440305/1‑228 (MQ=255)
          gtgtACTGCTCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGGATCGTGACGCAgccgg                                                                                                                                                                                                                     >  1:55207/1‑227 (MQ=255)
                  cTCCTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCAACGTCAATGGTGATGCCCTGTTCGCGCACAGCTTGCCGGCCGTCAACCAGCAGAGCCAGACCCAGCTTTTCGCCCTGGGAGCCGGGACCATTATTGTCGTATTCTAGCGATCAGACCTGACCTTCGCAGATTAGCCGGGGATCGTTCGGCAggcgaacccaaca                                                                                                                                                                                                             >  2:510627/1‑220 (MQ=255)
                       gCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTTGCGGGGAACGTGCAGCACAAGACCAAACaaaaaaa                                                                                                                                                                                                        >  2:94296/1‑222 (MQ=255)
                       gCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGGGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGGTTCGTGCAGCAGACGGCACATCAGAGGAc                                                                                                                                                                                                        >  2:381005/1‑228 (MQ=255)
                                  tCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGGATCGGGCCACACACGAACAATACGAATACtttttctcgtc                                                                                                                                                                                             >  2:557422/1‑221 (MQ=255)
                                    ggggATAATAAATTTCCGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATc                                                                                                                                                                                           <  1:59702/225‑1 (MQ=255)
                                    ggCGTTTATAAATTTCCGCTGCCGGGTAGAGAAAGTGCGGTGGGCCCCGTCAAGGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTGCTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATc                                                                                                                                                                                           <  2:123377/228‑1 (MQ=255)
                                    ggCGATAATTAATTTCCCCTTCTCGGTAGAGAAAGTGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATc                                                                                                                                                                                           <  2:440305/228‑1 (MQ=255)
                                    ggCGATAATCAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATc                                                                                                                                                                                           <  1:403031/228‑1 (MQ=255)
                                     gggATAATTAATGTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTTCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATGTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATc                                                                                                                                                                                           <  1:110139/225‑1 (MQ=255)
                                     gCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTGAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATc                                                                                                                                                                                           <  1:260694/227‑1 (MQ=255)
                                                          tggggaggaaaaTCGCGGTAGGCCACGTCAAAGGGGATGCCCTGTTCGCGCCCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGaaaa                                                                                                                                                                     <  2:469719/219‑1 (MQ=255)
                                                                 agaAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTAATTTTGCCGTCATCGACGCTACCCCAGGCCCGCAAACgcagaa                                                                                                                                                              >  1:271937/1‑226 (MQ=255)
                                                                   aaaTAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGGATCGTGCAGCAGACGACCAATCAGAGTACTTTTTCCGTCATCGGCGCTACCACAGGTCAGAAAACGCCGCAgg                                                                                                                                                            >  1:469719/1‑228 (MQ=255)
                                                                      tAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGGATCGTGCAGCAGACGACCAATCATAGTACTTTTGCCGTCATCGACGCTTCCACAGGTCAGAAAACGCCGAAgggt                                                                                                                                                          >  1:649393/1‑225 (MQ=255)
                                                                      tAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGGATCGTGCAGCAGACGACCAATCAGAGTACCTTTTCCGTCATTGCCGCTACCCCAAGACCGGAAACgccgcacgtt                                                                                                                                                          >  1:136968/1‑223 (MQ=255)
                                                                          ggTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCAAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTGCTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGAGTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTa                                                                                                                                                      <  1:409245/227‑1 (MQ=255)
                                                                            tAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTgcagca                                                                                                                                                                                                                          <  1:453708/157‑1 (MQ=255)
                                                                             aGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCCGCAGACGACCAATCAGAGGACTTTTTCCGTCCTCGACGCTACCCCAGGTCAGAAAACGCAGCAGGCTTTTTtgtt                                                                                                                                                  >  1:376657/1‑228 (MQ=255)
                                                                                            ggaggCCCTTGTCGCGCTCAGCTTGGCGCCCGTCCCCCAGCAGGGCCAGATCCAGCTTGTCGCCCTGGGTGCCGTGACGCTTACTGGCGTTATGCAGCGATGAGAGCTGAGCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAatcatc                                                                                                                                    <  2:376657/223‑1 (MQ=255)
                                                                                              aTGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGGACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGAAGCAGGCTTTTATGTTGTTGCGCAATCATCCa                                                                                                                                  >  1:266303/1‑227 (MQ=255)
                                                                                              aTGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGGACTTTTGCCGGCATCGACCCCACCACAAGTCAAAAAACACAGCAGGCTTTTTTGTTTTTGCGCAACCATACAg                                                                                                                                 >  2:591776/1‑228 (MQ=255)
                                                                                               tGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGCGCTGATCTTCGTAGATTTGGCGGGTATCGTGCCGCAGACGCCCAATCCGAGTACTTTTTCCGTCACTGACGCTACCACCGGCCAGAAAACGCAACCGGCTTTTATGTTTTTGCGACATCCTCCAg                                                                                                                                 >  2:340679/1‑227 (MQ=255)
                                                                                                    tGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTTCCGTCATCGACGCTAcc                                                                                                                                                                                  >  2:307739/1‑173 (MQ=255)
                                                                                                     gTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGGGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAACAGGGTTTTATGTTGTTGCGCAATCACCCAGTGTTCGa                                                                                                                          >  1:464050/1‑228 (MQ=255)
                                                                                                     gTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACca                                                                                                                                                                                 <  1:307739/173‑1 (MQ=255)
                                                                                                      ttCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGCCAGCAAACGCAGCAGGCTTGTATGTTGTTGCGCAATCATCCCAGCtttcac                                                                                                                         >  1:458435/1‑224 (MQ=255)
                                                                                                                gCTTGCAGGCCTGCCCCCGGCAGGACCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTGACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCa                                                                                                               <  2:8576/228‑1 (MQ=255)
                                                                                                                 cTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAatcatc                                                                                                                                    <  1:298157/206‑1 (MQ=255)
                                                                                                                         ccGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGTGGATCGTGCAGCAGACGGCCAACCAGAGTACTTTTTCCGTCATCGACGCTACCCCAGGTCCGAAAAAACAGCAGGGTTTTTTGTTTGTTCGCAATCCTCCAGGCTCTCAAGCCCCCTTCATTgggagtt                                                                                                      >  2:81268/1‑223 (MQ=255)
                                                                                                                         ccGTCCACCAGCAGACCCAGTTCCATCTTTCCGCCGGGGGGGCCGTGACTGTTACTGTCTTTAGTCAGCGATGAGACCTGATCTGTGTAGATTTGGCGGGGATCGTGCAGCAGACGCCCAATCAGAGTACTTTTGCCGTCAGCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGTGCAATCATCCAGGCTTCGAGGCCGCGTTCATTGGCGAtt                                                                                                      <  2:410154/228‑1 (MQ=255)
                                                                                                                           gTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGGATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTCCATTGGCGATttg                                                                                                    >  1:195346/1‑228 (MQ=255)
                                                                                                                            tCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGTCGGGGATCGTGCAGCAGACGGCCAATCAGGATACTTTTTCCGGCATCCACGCTTCCACAGGGCAGAAAACGCAGCAGGCTTTTATTTTGTTGTGCAACAATCCAGGCTTCGACGCCGCCTTCATTGGCGTTttgt                                                                                                   >  2:456534/1‑228 (MQ=255)
                                                                                                                                 cggcagGGCCAGGTCCAGCTTTTCGCCCTGGGTGCCGTGACGTTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAGCGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGc                                                                                              <  1:636276/226‑1 (MQ=255)
                                                                                                                                    cagAGCCAGATCCGGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCaa                                                                                            <  1:381005/227‑1 (MQ=255)
                                                                                                                                      ggggcAGATCCAGGTTTTCGCCCTGGGGGCCGTGACCGTTACGGTCGTTATGCAGCGATGAGGGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTg                                                                                         <  2:95125/224‑1 (MQ=255)
                                                                                                                                       agCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCCTCGACGCTACCACAGGTCCGAAAAAGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCCAgggg                                                                                        >  1:190975/1‑225 (MQ=255)
                                                                                                                                              tCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATAATCCAGGCTTCGCCGCCGCCTTCATTGGCGATTTGTTGTGCAAGt                                                                                          >  1:95125/1‑219 (MQ=255)
                                                                                                                                                   cTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGGTCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACCGGTCAGAAAACCCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACCCCGCCCTCCTTGGCGATTTGTTTTGGAAATgggggggtacactt                                                                            >  2:67516/1‑218 (MQ=255)
                                                                                                                                                        cGCCCTGGGTGCCTGGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTaaaaa                                                                       <  1:340679/228‑1 (MQ=255)
                                                                                                                                                         gCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTaaaaa                                                                       <  1:462783/227‑1 (MQ=255)
                                                                                                                                                             tGGTTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTTCAGCAGACGACCAATCATAGTACTTTTGCCGTCATCGACGCTACCACAGGGCAGAAAACGCAGAAGGCTTTTTTGTTGGTGCGCAATCATACAGGCTTCGACACCCCCTTCCTTAGCGATTTGTTTTGCAAGTGGGTGGTGCATCTACAAAAAAccc                                                                  >  2:600721/1‑228 (MQ=255)
                                                                                                                                                                gTGCCGTGACTGTGACTGTCGTTATGCAGAGATGAGGGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTGCTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCCTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCt                                                                <  1:294564/227‑1 (MQ=255)
                                                                                                                                                                         cggttccttttgTTATGCAGCGGTGAGCGCTGATTTTCGGAGGTTTGGCGGGTCTCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTTGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGtttttt                                                      <  2:266303/217‑1 (MQ=255)
                                                                                                                                                                            tttcggTCGTTTAGCAGCGATGAGGGGTGATTTTCTGAGATTTGGCGGGTATCGTGCAGCAGACGACCAAGCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCa                                                    <  2:562486/222‑1 (MQ=255)
                                                                                                                                                                             ttctttggTTTTGCGGGGATGGGAGCTGATCTTCGTTGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAg                                                   <  2:297041/220‑1 (MQ=255)
                                                                                                                                                                             tACTGTCGTTATGCAGCGATGAGAGCTGATTTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGc                                                  <  1:94296/228‑1 (MQ=255)
                                                                                                                                                                             tACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTCTGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTTTTCATCTTACAAATACCCCTGACGTTTTTTCAg                                                   >  1:608548/1‑227 (MQ=255)
                                                                                                                                                                              aCTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTgcagc                                                                                                                                                                                                                           >  2:353397/1‑58 (MQ=255)
                                                                                                                                                                               cTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTgcagca                                                                                                                                                                                                                          <  1:353397/58‑1 (MQ=255)
                                                                                                                                                                                             cGATGAGAGCTGATCTTCGTAGATTTGGCGGGGATCGTGGAGCAGACGACCAATCAGAGTACTTTTGCCGTAATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACCCCCCCTTCATTGGCGATTTTTTGTGCCAGGGCGGGGTTCACCTTAAAAATACCCCCGAAATTTTTTAACGTCCAttaaaccccccg                                  >  2:193030/1‑217 (MQ=255)
                                                                                                                                                                                                 tatggttgcttttCGTTGATTTGGCGGGTTGCGTCCGCCAAACCACCAAACAGAGTCCTTTGGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGTTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGCCCCCGCCTGGTc                              <  2:271937/217‑1 (MQ=255)
                                                                                                                                                                                                    accttgtcttctGAGATTTGGCGGGGATCGTGCAGCAGACGACCAATCCGAGTCCTTTTGCCGTCATCGACGCTCCCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATCCCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTcgcg                           <  1:207677/222‑1 (MQ=255)
                                                                                                                                                                                                              cGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCCCAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCCTCCAGGCTTCAACGCCGCCTTCATTGGGGTTTTTTTGGGCAATACCGGGGTTCTCTTTTAAAATACCCCTGACGTTTTTTTAGCTCCCAAACCCCCGCCTGTGCGCggcgaacaccc                 >  2:607970/1‑219 (MQ=255)
                                                                                                                                                                                                              cGTAGATTTGGCGGGGATCGTGCAGCAGACGACCAATCAGGGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGCAGCCGCCTTCCTTGGCGATTTGTTTGTCAAGTGGGGGGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTCGTCGCGGGCAATCAcc                 >  2:655658/1‑227 (MQ=255)
                                                                                                                                                                                                               gTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACGTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCAcg                 <  1:165933/227‑1 (MQ=255)
                                                                                                                                                                                                               gTAGATTTGGCGGGTATCGGGCAGCAGACGACCAATCAGAGTTCTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATCCCCCTGACGtttttt                                                      <  1:645734/190‑1 (MQ=255)
                                                                                                                                                                                                                      tGGGGGGTGACGGGCAGCAGGCGGCCAATCAGGGTTATTTTGCCGGCACGGACGCTACCACAGGGCAGAAAACGCAGCAGGCTTTTTTGTTTGTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTg         <  2:464050/228‑1 (MQ=255)
                                                                                                                                                                                                                         cGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGt                               <  1:516842/203‑1 (MQ=255)
                                                                                                                                                                                                                             tATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGGTTTGTTGTGCAAGTGCGGGGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCCCCTGGTCGCGGTCAATCACCCGGCCCTGACGTTCa  >  1:41388/1‑228 (MQ=255)
                                                                                                                                                                                                                             |                                                                                                                                                                                                                                   
CATATTGCGGGTGTACTGCTCGTGCCCTGGGGTGTCGGCGATAATAAATTTACGCTTCTCGGTAGAGAAATAGCGGTAGGCCACGTCAATGGTGATGCCCTGTTCGCGCTCAGCTTGCAGGCCGTCCACCAGCAGAGCCAGATCCAGCTTTTCGCCCTGGGTGCCGTGACGCTTACTGTCGTTATGCAGCGATGAGAGCTGATCTTCGTAGATTTGGCGGGTATCGTGCAGCAGACGACCAATCAGAGTACTTTTGCCGTCATCGACGCTACCACAGGTCAGAAAACGCAGCAGGCTTTTATGTTGTTGCGCAATCATCCAGGCTTCGACGCCGCCTTCATTGGCGATTTGTTGTGCAAGTGCGGTGTTCATCTTAAAAATACCCCTGACGTTTTTTCAGCTCCATAGACCCCGCCTGGTCGCGGTCAATCACGCGGCCCTGACGTTCA  >  NC_000913/2875048‑2875496

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