Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,472,5080AG78.0% 69.8 / 12.2 41V9A (GTC→GCC) pyrLpyrBI operon leader peptide
Reads supporting (aligned to +/- strand):  ref base A (5/4);  new base G (16/16);  total (21/20)
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.47e-01
Rejected as consensus: Frequency below/above cutoff threshold.

CAATGACTTTGTGCTTCAACAGCTCTGGTTGCGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGC  >  NC_000913/4472260‑4472751
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   
cAATTACTTTGTGCTTCAACAGCTCTGGTTGCGGGTTTGTTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCGGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGTTAAGACa                                                                                                                                                                                                                                                   <  2:404851/251‑1 (MQ=255)
           tgCTTCAACAGCTCTGGTTGCGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGACAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGGCAAAATGTCGAac                                                                                                                                                                                                                                        >  1:555268/1‑251 (MQ=255)
            gCTTCAACAGCTCTGGTTGCGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGGTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCGTGTTTCAGACGCGGGAAGGCAAAATGTCGAaca                                                                                                                                                                                                                                       <  2:420089/251‑1 (MQ=255)
                          ggTTGCGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGGCAAAATGTCGAACACACTGAACCATAc                                                                                                                                                                                                                          <  1:404255/250‑1 (MQ=255)
                                  gTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGGCAAAATGTCGAACACACTGAACCATACAtcctcctg                                                                                                                                                                                                                 <  2:427494/251‑3 (MQ=255)
                                    ttGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGGCAAAATGTCGAACACACTGAACCATACATCCTCCCgg                                                                                                                                                                                                                <  2:569488/250‑1 (MQ=255)
                                     tGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGc                                                                                                                                                                                                               <  1:113631/250‑1 (MQ=255)
                                                  gccgcTGTCGCCAGCACCAGATTAAGGTCATCGCGGATAAGGTCGTTTATGGAAATGATATGTTTCTTAAAAAGCGGATTAGCCATCTTTTATCTCCTGCCGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGAGCAACGGGAAGAAAAACGGCAGGCCAGCGTCTGTGTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAGCCATACACCCCCCCGGCAAATTGTCCGGcgt                                                                                                                                                                                                 <  2:245873/251‑3 (MQ=255)
                                                          cGCCAGCACCAGATTAAGGGCATCGCGAATAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGGCAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTc                                                                                                                                                                                          <  2:434733/250‑1 (MQ=255)
                                                            ccAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGGCAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATc                                                                                                                                                                                       <  1:150461/251‑1 (MQ=255)
                                                                             tCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGGCAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGATAGAAgcggcgaa                                                                                                                                                                        >  2:280876/1‑247 (MQ=255)
                                                                                       tAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGACATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGGCAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGc                                                                                                                                                            >  2:256917/1‑251 (MQ=255)
                                                                                           gTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGGCAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAata                                                                                                                                                        <  1:43637/251‑1 (MQ=255)
                                                                                                  aTGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGGCAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTc                                                                                                                                                 >  2:178847/1‑251 (MQ=255)
                                                                                                               gTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGTATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGGCAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGTTGCAATATTTTGTCAGAGCAAACGGtt                                                                                                                                    >  2:254038/1‑251 (MQ=255)
                                                                                                                  tCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGACTGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGGCAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTAt                                                                                                                                 >  2:85372/1‑251 (MQ=255)
                                                                                                                    tGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAACAAAAAGACCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCCGCGACTTTTTTCAGACGCGGTAAGGCAAAATGTCGCACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCCTTTGATGCAATATTTTGTCAGCGCACACGGTTTAttt                                                                                                                               >  1:608785/1‑251 (MQ=255)
                                                                                                                    tGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGGCAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTAttt                                                                                                                               >  2:316318/1‑251 (MQ=255)
                                                                                                                    tGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTAttt                                                                                                                               >  2:476481/1‑251 (MQ=255)
                                                                                                                        ataGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGGCAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGaa                                                                                                                            <  2:247637/250‑1 (MQ=255)
                                                                                                                            cGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGGCAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCTGCATTATTCTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTTAATTaaa                                                                                                                       <  1:64895/251‑1 (MQ=255)
                                                                                                                                           ttATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGGCAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATg                                                                                                         >  2:276602/1‑250 (MQ=255)
                                                                                                                                                     aCGCCTGGTAAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACCCACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTAttt                                                                                              <  1:342127/251‑1 (MQ=255)
                                                                                                                                                         cTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCCGGCCAGCGTCTTTTTTCAGACGCGGTAAGGCAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGAtg                                                                                          >  1:405058/1‑251 (MQ=255)
                                                                                                                                                         cTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGGCAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTCTACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCACTATTTTGTCAGCGCAACCGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGAtg                                                                                          >  2:249857/1‑251 (MQ=255)
                                                                                                                                                         cTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGGCAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATCTTTTGTCAGCGCAAAAGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGAtg                                                                                          >  2:138610/1‑251 (MQ=255)
                                                                                                                                                         cTGGGCAAAAAAAAGCACCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGAGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGAtg                                                                                          >  1:461391/1‑251 (MQ=255)
                                                                                                                                                               aaaaaaaaGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGCAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTc                                                                                    >  1:570029/1‑251 (MQ=255)
                                                                                                                                                                         cccTCGATTGCGGGGCTGGGACTGGGTGATCAACGGGCAGAAAACCGGCAGGCCAGAGTCTTTTTTCAGCCGCGGTAAGGCAAACTGTCGAACACACTGCAACATACATCCTCCAGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCCCAAGCCTTTGATGCAATATTTTGTCTGCGACAACGGTTTATTTGAATTAAAAGTCAAGGAATATTCATTTTTATTTGATGTGATTCTGCAGGGGa                                                                           >  2:588499/1‑250 (MQ=255)
                                                                                                                                                                         cccTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGGAAGGCAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGCAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGa                                                                           >  2:119831/1‑250 (MQ=255)
                                                                                                                                                                               aTTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGGCAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTcc                                                                    <  1:102127/251‑1 (MQ=255)
                                                                                                                                                                                          gggAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGGCAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCaa                                                         <  1:178847/251‑1 (MQ=255)
                                                                                                                                                                                                  gTGAGCAACGGGAAGAAAAACGGCAGGCCAGCGTCTGTATTCAGACGCGGTAAGGCAAAATGGCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGGCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACGGTCCTTCGGTATCAATACTTGCa                                                 <  2:555268/251‑1 (MQ=255)
                                                                                                                                                                                                         aCGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGGCAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTc                                          <  1:283309/251‑1 (MQ=255)
                                                                                                                                                                                                             gaagaaAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGGCAAAATGTCGACCACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGGCCTGTCCTTCGGTATCAATAATTGCAGACATTCCCTg                                      >  2:101180/1‑251 (MQ=255)
                                                                                                                                                                                                             gaagaaAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCTg                                      >  1:24349/1‑251 (MQ=255)
                                                                                                                                                                                                                   aaaCGGCAGGCCTGCGTCTTTTTTAAGTCGCGTTAATGCAAAATGGCGAACACACTGAACAATACATCCTCCCGGTAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGTAAACGGTTTATTGGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGCGAACTGTCCTTCTGTAGCAATAATTGCATACATTCCCTGcttttt                                <  2:350341/251‑1 (MQ=255)
                                                                                                                                                                                                                               aGCGTCTTTTTTCAGACGCGGTAAGGCAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCTGCTTTTTCCTTtttcactt                    >  1:378520/1‑251 (MQ=255)
                                                                                                                                                                                                                                          tCAGACGCGGTAAGGCAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATc          <  1:256917/250‑1 (MQ=255)
                                                                                                                                                                                                                                                gcgGTAAGGCAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATg   <  1:488756/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 cgGTACGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGc  >  2:282652/1‑251 (MQ=255)
                                                                                                                                                                                                                                                        |                                                                                                                                                                                                                                                   
CAATGACTTTGTGCTTCAACAGCTCTGGTTGCGGGTTTGCTTTCAGTTTCGCCGCTGTCGCCAGCACCAGATTAAGGTCATCGCGACTAAGGTCGTTTATGGAAATGATATGTTTCTGATATAGCGGATTAGCCATCTTTTATCTCCTGACGCCTGGGCAAAAAAAAGCCCCTCGATTGAGGGGCTGGGAATGGGTGATCAACGGGAAGAAAAACGGCAGGCCAGCGTCTTTTTTCAGACGCGGTAAGACAAAATGTCGAACACACTGAACCATACATCCTCCCGGCAAATTGTCCGGCATTATACTCATCGTCAGAAGCGGCGCAAGCATTTGATGCAATATTTTGTCAGCGCAAACGGTTTATTTGAATTAAAAGTCAAGGTATATGCATTTTTATTTGATGTGATTCTGCAGGGGAACTGTCCTTCGGTATCAATAATTGCAGACATTCCCTGCTTTTTCCTTTTTCACTTTCACGCAATCAGATATGC  >  NC_000913/4472260‑4472751

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