Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,838,222 A→T 7.5% intergenic (‑117/+32) hydN ← / ← ascG formate dehydrogenase‑H, [4Fe‑4S] ferredoxin subunit/asc operon transcriptional repressor; prpBC operon repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,838,2220AT 7.5% 157.4 / 3.7 53intergenic (‑117/+32)hydN/ascGformate dehydrogenase‑H, [4Fe‑4S] ferredoxin subunit/asc operon transcriptional repressor; prpBC operon repressor
Reads supporting (aligned to +/- strand):  ref base A (23/26);  new base T (2/2);  total (25/28)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.03e-01

ACACTGACGGCAAACTGTCGCCGTGGAAATGTTCACACCTTTAATGACATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCGAATCCGATAACTGACACCTGCTCTGGCACCGCTACGCCGCGCTCGTGTAACGCTTTCATCGCACCTATCGCCATATCGT  >  NC_000913/2837923‑2838519
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                         
tcacTTACGGCAAACTGTCGCCGTTGAAATGTTCACACCTTTAATGACATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGac                                                                                                                                                                                                                                                                                                          <  2:411948/300‑1 (MQ=255)
 cacTGACGGCAAACTGTCGCCGTGGAAATGTTTACACCTTTAATGACATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCAAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATATCATCGaaa                                                                                                                                                                                                                                                                                                         >  1:411385/1‑299 (MQ=255)
  ccTGACGGCAAACTGTCGCCGTGGAAATGTTCACACCTTTAATGACATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACt                                                                                                                                                                                                                                                                                                       <  2:339874/300‑1 (MQ=255)
          cAAACTGTCGCCGTGGAAATGTTCACACCTTTAATGACATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGTACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGTCACTAGTgtct                                                                                                                                                                                                                                                                                                >  1:54336/1‑298 (MQ=255)
           aaaCTGTCGCCGTGTAAATGTTTACTCCGTTAATGTCTTGGATACGCGGTAAAAAAGTTGCCGGGGTCAGTGTTGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGa                                                                                                                                                                                                                                                                                               <  1:371061/301‑1 (MQ=255)
             ctacacgacgctcttccgatctcacCTTTAATGACATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGTCACTAGTGTCGATg                                                                                                                                                                                                                                                                                             <  1:19272/279‑1 (MQ=255)
               tggtgCCGTTGAAATGTTCACACCTTTAATGACATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACGCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGAc                                                                                                                                                                                                                                                                                           <  1:575382/297‑1 (MQ=255)
                  cGCCGTGGAAATGTTCACACCTTTAATGACATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGAAACTAGTAACGATGACAtg                                                                                                                                                                                                                                                                                        >  2:322800/1‑301 (MQ=255)
                   gCCGTGGAAATGTTCACACCTTTAATGACATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACAtgt                                                                                                                                                                                                                                                                                       >  2:605419/1‑301 (MQ=255)
                       gggAAATGTTCACACCTTTAATGACATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGtgatg                                                                                                                                                                                                                                                                                   <  1:343453/300‑1 (MQ=255)
                        tgAAATGTTCACACCTTTAATGACATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGTACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGTCACTAGTGTCGATGACAtcgttttt                                                                                                                                                                                                                                                                                  <  2:54404/300‑8 (MQ=255)
                           aaTGTTCACACCTTTAATGACATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGaaaa                                                                                                                                                                                                                                                                               >  2:590971/1‑301 (MQ=255)
                                tCACACCTTTAATGACATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGttt                                                                                                                                                                                                                                                                          >  2:501552/1‑301 (MQ=255)
                                   cacCTTTAATGACATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGTCACTAGTGTCGATGAGAtcggaagagcacacgtctg                                                                                                                                                                                                                                                                       >  2:19323/1‑283 (MQ=255)
                                             aatagggttACGCGGTAAATAATGTTCCGGGGTCATCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGAACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGc                                                                                                                                                                                                                                                             <  1:467193/293‑1 (MQ=255)
                                                 tGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGCGCAAtg                                                                                                                                                                                                                                                         >  1:417078/1‑301 (MQ=255)
                                                       cgcgGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAAtgagtgaa                                                                                                                                                                                                                                                   >  2:63860/1‑300 (MQ=255)
                                                          ggTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATAATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTcg                                                                                                                                                                                                                                                >  2:550960/1‑301 (MQ=255)
                                                                aaaGTTTCCGGGGTCAGCGTTGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACAATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGcggcgg                                                                                                                                                                                                                                          <  1:544852/301‑1 (MQ=255)
                                                                           ggTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCAGATCAGGTTTcc                                                                                                                                                                                                                               >  2:626290/1‑301 (MQ=255)
                                                                                   aggCACCGTTCTGATTTTTCTGATGGGAAACCACGCACTCTACTTCACAGGTACGGCAACCAATACAGTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAAGGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAAGTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAgg                                                                                                                                                                                                                       <  1:625413/299‑1 (MQ=255)
                                                                                           tCCTGTTTTTCCTGTTGGGAAACCCCGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGGGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAAGTACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTcggc                                                                                                                                                                                                               <  1:500905/301‑1 (MQ=255)
                                                                                                 ttttCCTGATAGAAAACCAAGCACGCTACTTAACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAAGGATGAAACGGGTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTgag                                                                                                                                                                                                         <  1:322396/301‑1 (MQ=255)
                                                                                                          tgaGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCg                                                                                                                                                                                                <  2:211104/301‑1 (MQ=255)
                                                                                                                  cacgcacgCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATTAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTTATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAg                                                                                                                                                                                        <  1:604658/301‑1 (MQ=255)
                                                                                                                                  caGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGt                                                                                                                                                                        >  2:536565/1‑301 (MQ=255)
                                                                                                                                    ggTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTcc                                                                                                                                                                      >  2:317902/1‑301 (MQ=255)
                                                                                                                                                    cATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAata                                                                                                                                                      >  2:141996/1‑300 (MQ=255)
                                                                                                                                                        aattCGCGTCTGCAATGATGTAACGGTTCATTAAATTCTCCAGCAATTACAGTTAATTCGCCGATACATTCTCAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCAtct                                                                                                                                                  <  1:535892/298‑1 (MQ=255)
                                                                                                                                                          cTCGCGTCAGCAAGGATGAAACGGTTCATCAAATTTTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCa                                                                                                                                                <  1:63792/301‑1 (MQ=255)
                                                                                                                                                                 aaGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCg                                                                                                                                         <  1:141871/300‑1 (MQ=255)
                                                                                                                                                                   gCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATTTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGAAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGATTTCGGCGGTGAGAAATCCCCGCCATCGCGCATAAAAACCAGCCACCCAAAAATTTCCTACATCTTCTCAATCCCCGGa                                                                                                                                       >  1:321358/1‑301 (MQ=255)
                                                                                                                                                                           aaaaCGGTTCATCAAATTTTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTAc                                                                                                                               <  2:417641/300‑1 (MQ=255)
                                                                                                                                                                                    catcaAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGgagaa                                                                                                                      >  2:636763/1‑300 (MQ=255)
                                                                                                                                                                                      tcaAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGAAAAgcg                                                                                                                    >  2:228632/1‑301 (MQ=255)
                                                                                                                                                                                                         cAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGc                                                                                                 <  2:480754/301‑1 (MQ=255)
                                                                                                                                                                                                               ttGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGc                                                                                           <  2:31107/300‑1 (MQ=255)
                                                                                                                                                                                                               atgtggcGATAAATTTATAGATCATGCCAGTTTTTTATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCAGAGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGGCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTACTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGc                                                                                           <  1:439465/295‑1 (MQ=255)
                                                                                                                                                                                                               aTGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGc                                                                                           >  2:424931/1‑301 (MQ=255)
                                                                                                                                                                                                                tGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCCGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGACGGAGCAATGCGTGAGTCTCGTCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGCAATCCCCGCCATCGAGCATAAACATCAGCCGTCCAATAATTTCCTGAATCATCACAGTCACCGAAATTTTTACGCTGGAGCACGCCGGAAACGTGTAAGGGGCAATAGca                                                                                          >  2:186386/1‑300 (MQ=255)
                                                                                                                                                                                                                         ctcgctcttccgatctccAGTTTTTAATTTTCTGTTTTTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCg                                                                                 <  1:3140/285‑1 (MQ=255)
                                                                                                                                                                                                                            ttCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGACTCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGAGTAGGGGGCAATAGCGAAATCATCGAAt                                                                              >  1:599412/1‑301 (MQ=255)
                                                                                                                                                                                                                              cacaAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCGAATcc                                                                            >  2:330562/1‑301 (MQ=255)
                                                                                                                                                                                                                               agaaaTCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCGAATCCg                                                                           <  1:424368/299‑1 (MQ=255)
                                                                                                                                                                                                                                caAATCATGCCCGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGAAATCATCAAATCCgc                                                                          >  1:593883/1‑300 (MQ=255)
                                                                                                                                                                                                                                        gCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCGAATCCGATAACTGac                                                                  <  1:282334/301‑1 (MQ=255)
                                                                                                                                                                                                                                         ccAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCGagatcggaagagcaca                                                                 >  2:3140/1‑286 (MQ=255)
                                                                                                                                                                                                                                                         ttttttttAAGGTAATTAATTTCTTTAATGCAGGAAAAACGATGTCATCGAAACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTTAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCAGCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCGAATCCGATAACTGACACCTGCTCTGGCACCGc                                                 <  1:590210/296‑1 (MQ=255)
                                                                                                                                                                                                                                                                 aaGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCGTCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCGAATCCGATAACTGACACCTGCTCTGGCACCGCTAcgccgc                                         >  2:29151/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                        ttaattTTTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCGAATCCGATAACTGACACCTGCTCTGGCACCGCTACGCCGCGCTCgtg                                  <  2:321762/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                   aaTGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCGAATCCGATAACTGACACCTGCTCTGGCACCGCTACGCCGCGCTCGTGTAACGCTTTCa                       <  2:600173/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                   aaTGCAGGAAAAAAGTGGGTGTCGACACTAGTGAGGATGACAGGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCGAATCCGATAACTGACACCTGCTCTGGCACCGCTACGCCGCGCTCGTGTAACGCTTTCa                       <  1:317498/301‑1 (MQ=255)
                                                                                                                                                                                                                                                                                          gAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCGGAAGGTTTTCGGCGGTGAGACATCCCCGCCATCGAGCATAAAGATCAGCGGTCCAATAATTGCGTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCGAATCCGATAACTGACACCTGCTCTGGCACCGCTACGCCGCGCTAGTGTAACGCTTTCATCGCAcc                >  1:249465/1‑301 (MQ=255)
                                                                                                                                                                                                                                                                                                    ttcttCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCGAATCCGATAACTGACACCTGCTCTGGCACCGCTACGCCGCGCTCGTGTAACGCTTTCATCGCACCTATCGCCata      <  2:198866/297‑1 (MQ=255)
                                                                                                                                                                                                                                                                                                        tCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCGAATCCGATAACTGACACCTGCTCTGGCACCGCTACGCCGCGCCCGGTTAACGCTTTCATCGCAcccacccccaaaaacc  >  1:198623/1‑287 (MQ=255)
                                                                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                                                                         
ACACTGACGGCAAACTGTCGCCGTGGAAATGTTCACACCTTTAATGACATGGATACGCGGTAAAAAAGTTTCCGGGGTCAGCGATGCACAGTCCTGATTTTCCTGATGAGAAACCACGCACGCTACTTCACAGGTACGGCAACCAATACATTTACTCGCGTCAGCAATGATGAAACGGTTCATCAAATTCTCCAGCAATGACAGTTAATGCGCCGATACATTCACAAATCATGCCAGTTTTTAATTTACTGTTATTTAAGGAAATTAATTTCTGTAATGCAGGAAAAACGATGTCATCGACACTAGTGACGATGACATGTGATGACAATGTTTATCGCGAAGGAGCAATGAGTGAGTCGCGGCGGATCAGTTTTCCGCTGAAGGTTTTCGGCGGTGAGAAATCCCCGCCATCGAGCATAAAAATCAGCCGTCCAATAATTTCCTGAATCATCTCAGTCACCGGAATTTTTACGCTGGAGAGCGCCGGAACGGTGTAGGGGGCAATAGCGATATCATCGAATCCGATAACTGACACCTGCTCTGGCACCGCTACGCCGCGCTCGTGTAACGCTTTCATCGCACCTATCGCCATATCGT  >  NC_000913/2837923‑2838519

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 9 ≤ ATCG/ATCG < 20 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: