Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,415,3950GA25.4% 122.8 / 12.2 64R664R (CGC→CGTrecEexonuclease VIII, ds DNA exonuclease, 5' ‑‑> 3' specific
Reads supporting (aligned to +/- strand):  ref base G (13/34);  new base A (16/0);  total (29/35)
Fisher's exact test for biased strand distribution p-value = 1.85e-07
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.89e-01
Rejected as consensus: Frequency below/above cutoff threshold.

GGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATACCGGGACCCGCGTGGTAATTCTCATTCGAAATTCCGTAATAAATACCTGGCTCTATGTCTTCTACGATTACGGGATCTGCGACTTCG  >  NC_000913/1415147‑1415635
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                            
ggACGGCTCCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCg                                                                                                                                                                                                                                                            <  1:294415/251‑1 (MQ=255)
ggACGGCCCCGAAACAAAATTCCTGTTTCAGGATCTTCCCATTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCg                                                                                                                                                                                                                                                            <  2:498803/251‑1 (MQ=255)
ggACGGAACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTGTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCg                                                                                                                                                                                                                                                            <  2:60828/251‑1 (MQ=255)
           ccactaaaTTCCTGTTGCAGGACCTCCCCATAAAATTGAGTATTCAGCGTTGCCGGCGCTCTCCACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATACTTGAGTTCAATTTTCCGGCCTTCTGCCGCAGTGATAACCGTGTTTCCTGTGCTTGCGCATTCCATCAGAAACTCTTTCCCTTCAACTTCTCCGGCGTTTGTACGGCGGGTAAATTCATGTGCTACGATAAAGCGGTTACTGAAtt                                                                                                                                                                                                                                                 <  2:740699/246‑1 (MQ=255)
                  aTTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTTCGGTAAAACGGGTTCTGAAttcttctgg                                                                                                                                                                                                                                         >  1:722073/1‑248 (MQ=255)
                    tCCTGTTTCAGGATCTTCCCAGTAAATTGAGGAGTCAGCGTGTCAGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATACATGATTTCAATGTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTGCGATAAAGCGGTTACTGAATTCTTCCGGt                                                                                                                                                                                                                                        <  1:631270/250‑1 (MQ=255)
                           tCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTa                                                                                                                                                                                                                               <  1:152076/250‑1 (MQ=255)
                            cAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTAcc                                                                                                                                                                                                                             <  1:372757/251‑1 (MQ=255)
                                    tCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGc                                                                                                                                                                                                                         <  1:67586/247‑1 (MQ=255)
                                    tCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCa                                                                                                                                                                                                                        <  1:507118/248‑1 (MQ=255)
                                    tCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTGTGATACATGAGTTCAATTTGCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCGTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGt                                                                                                                                                                                                                     <  2:303573/250‑1 (MQ=255)
                                               tgatgaTTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTGCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGtt                                                                                                                                                                                                        <  1:113102/250‑1 (MQ=255)
                                                 atgatTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCt                                                                                                                                                                                                     <  1:225040/251‑1 (MQ=255)
                                                         gCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTTAATTCTTTCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAAtag                                                                                                                                                                                               >  1:491326/1‑247 (MQ=255)
                                                                         caacaaGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTTAATTTAAGTGCTACGGTAAAGAGGGTACTGGATTTTTTCGGGTCCAGAACCCCGCCGTGGAAAGCAGGTCCTAAAACGAGGGTGTTTGTCtttttg                                                                                                                                                                             >  1:687559/1‑248 (MQ=255)
                                                                           acaaGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTgtggt                                                                                                                                                                          <  1:733447/251‑1 (MQ=255)
                                                                             aaGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTACGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGTCGGTTTAAATACAGTGCTTCGATTAAGCCGTTACTGATTTCTTTCGGTTCTAGTACCCCGCTGTGGGAGGCACTTCCTACATAGAGCACTTTTGTCtttttgtgtg                                                                                                                                                                         >  2:12/1‑244 (MQ=255)
                                                                               gCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATTCATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCTTTTTGCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTGTTTGTCTTTGTGGTGTcc                                                                                                                                                                    <  1:275871/251‑1 (MQ=255)
                                                                                ccATTTCCCCAGCGGCAAAGCCATAACGCTCTGATACATGAGTTCAATTTTCCGGCCGTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGCTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGTCAGTGGAAAGCAGTTCCTAAAACGAGCGTTTTTGTCTTTGTGTTGTCCa                                                                                                                                                                   <  2:408408/251‑1 (MQ=255)
                                                                                ccATTGCCCCAGCGGCAAAGCCATAACGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCa                                                                                                                                                                   <  1:109367/251‑1 (MQ=255)
                                                                                    tGCCCCAGCGGCAAAGCCATAACGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAAATCAGGGGCTACGTTAAAACGGTTTCTGAATTTTTTCCGTTCAAGAAACCAGGAGTGGGAAACAGTTCCTAAAACGGGCGTTTTTGTTTTTGTGGTGTcccaggg                                                                                                                                                               >  1:678136/1‑246 (MQ=255)
                                                                                       cccAGCGGCAAAGCCATAACTCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGc                                                                                                                                                          <  1:331843/251‑1 (MQ=255)
                                                                                       cccAGCGGCAAAGCCATAACGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTTCGGTGCTACGATAAAACGGTTACTGAATTCTTCCGGGTCAAGGACCCGGCAGTGGGAAGCAGGTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACgggggg                                                                                                                                                          >  1:461558/1‑250 (MQ=255)
                                                                                                aaagactgaatgCTTTGATTCATGAGTTCAATTTTCCGGTCTTTTTCCGTAGTGATAACCGTTTTTCCGGTGCTTGCGCAATCCATCAGAAACGCTTTCTCTTCTTCTTTGACGGCGTTTGTTCGGCGGGTAAATTCAGGTGCTACGATAAATCGGTCACTGAATTCTTCCGGTTCAAGTACCCGGCATTTGAAAGCAGTTCCTAAATCGATCGTTTGTGCCTTTGTGGTGTCCCCGGTGGCATTTTTAc                                                                                                                                                 <  2:719968/239‑1 (MQ=255)
                                                                                                aaaGCCATAACGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTTAGGTGCTACGATAAAACGGTTACTGAATTCTTCCGGGTCAAGTAACCGGCAGGGGAAAACAGTTCCTTAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACg                                                                                                                                                >  1:642100/1‑251 (MQ=255)
                                                                                                     cATAACGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTTAATTTAGGTGCTACCGTAAAACGGTTACTAAATTTTTCCGGTTTAAGTACCCGGCAGTGGGAAGCAGTTTCTAAATCGCGCGGTTTTGTCTTTTTGGTTGTCACGGGGGCATTTTTTCGCcac                                                                                                                                            >  1:578454/1‑250 (MQ=255)
                                                                                                      aTAACGCTTTGATAAATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGGTTCAATTAAGGTGGAACATAAAACAGTTTTTGGATTCTTCCGGGTCCAGTACCCCGCCGCGGGAAAGGGTTCCGTCATCCAGCGTTTTTTTCTTTGTGGTTGCCCCGGGGGCGAttttttgccccc                                                                                                                                            >  2:714354/1‑243 (MQ=255)
                                                                                                       tAACGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTTAATTCAGGTGCTACGATAAAGCGGTTACTGTATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCCGTTTCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACaa                                                                                                                                          >  1:595994/1‑250 (MQ=255)
                                                                                                          cGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTTAATTCAGGTGCTACGCTAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTTTCTTTTTGGGGTCCACGGGGGCATTTTTACGCCACAAatat                                                                                                                                      >  1:850282/1‑251 (MQ=255)
                                                                                                            cTTTGATACATGAGTTCAATTTTCCGGCATTCATCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGGTTAATTCAGGTGCTCCGATAAAACGGGTACTGAATTATTCCGGATCAAGAACCCGACCGTGGAAAGCAGTTCCTCAATCGAGCAGTTTTGTCTTTTGTGTGTCACCGGGGGCATTTTTACGCCACAAataaa                                                                                                                                     >  1:751839/1‑248 (MQ=255)
                                                                                                             tttGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGAACGGCGGTTAAATTTCTGGGCTCCGATAAAGAGGGTACTGAATTCTTCCGGTTTAAGTACCCGGGAGGGGAAAGCAGTTACTCAATCGGGGGGTTTTGTCTTTTTGGAGTCCAAGGGGGGCTTTTTTGCCCACAAATATAgg                                                                                                                                   >  1:48459/1‑250 (MQ=255)
                                                                                                                gATACATGAGTTCAATTTTCCGGCCTTCCTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTTCCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGcc                                                                                                                                <  2:15669/251‑1 (MQ=255)
                                                                                                                  tACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTTAATTTCGGGTGTACGCTAAAACGGGTACTGAATTCTTCCGGCTCAAATACCCCGCAAGGGAAAGCAGGTACTAAATCGCGGGTTTTTTTCTTTTTGGTGTCCACCGGGGCAttttttcccccccaaaatattgtcgg                                                                                                                              >  2:299546/1‑230 (MQ=255)
                                                                                                                    cATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGAGCATTCCATCAGAAACGCCTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAAGTCTTCCGGTTCAAGTTCCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAg                                                                                                                            <  2:686052/251‑1 (MQ=255)
                                                                                                                     aTGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGt                                                                                                                           <  1:162650/251‑1 (MQ=255)
                                                                                                                            cAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTTAGGTGCTACGGTAAAACGGATAATGAATTCTTCCGGGTCCAGTAACCGGCAGAGGGAAACAGTTCCTAAAACGAGCGTTTTTGTCTTTGTGGTTGCCACCGGGGGCTTTTTACGCCACAAATTTAGAGCCgggggaacagc                                                                                                                     >  1:838049/1‑241 (MQ=255)
                                                                                                                              aTTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGGTTAATTTCAGGGGTCCGATAAAACGGTTTCTGAATTCTTTCGGTTTCAATACCCGGCAGGGGAAAGGCAGTCCTCAATCCAGCGTTTTTTTCTTTGTAGGGTCCCCAGGGGGATTTTTACGTCCCCAATATAGTGTCGGAggattaacaa                                                                                                                   >  2:186908/1‑241 (MQ=255)
                                                                                                                                    cGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGCAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTTAATTCCAGTGCTACCATCAAGAGGGTACTGAATTTTTTCCGATTCAGGACCCCGCAGCGGCAAAGCGGTTCTCAATCGTACGGTTTTTTCTTTTTGGGGTCCCCGGGGGCCTTTTTTCGCCACAAATATCGTGCCCGAGTGTCCGTAATGTCATc                                                                                                            >  1:483174/1‑251 (MQ=255)
                                                                                                                                     ggCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAACGGGTACTTAATTCTTCCGGGTCAAGTACCCGGCAGGGGAAAACAGTTCCTAAATCGAGCGGTTTTTTCTTTTTGGTGTCCCCGGGGGCCTTTTTACGCCACAAATATAGTGCCCGGGTATTAGCAATGTCATCg                                                                                                           >  1:734974/1‑251 (MQ=255)
                                                                                                                                        cttcttCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTTAATTCTAGTGCTACGATAAAACGGTTACTGAATTCTTCCGGGTCAAGTACCCCGCAGTGGGAAGCAGTTCCTTAATCGAGCGTTTTTGTCTTTTTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCCACAATGTCATCGAGc                                                                                                        >  1:23450/1‑251 (MQ=255)
                                                                                                                                          tcttcCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTTAATTCAGGTGCTACGGTAAAGCGGGTACTGAATTTTTCCGGTTCAAGTACCCGGCAGGGGAAAGCAGTTCCTTAATCGAGCGTTTTTTTCTTTTTGGTGTCCACGGGGGGCTTTTTTCGCCACAAATATAGTGCCg                                                                                                                               >  1:504081/1‑226 (MQ=255)
                                                                                                                                          tcttcCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCg                                                                                                                               <  2:504081/226‑1 (MQ=255)
                                                                                                                                                gCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCCGGTGGTTCGATAAAGCGGTTACTTTATTCTTCCGGTTCCCATAACCCGGGGGGGGAAGAAAGTTCTAAATCGCGCGTTTTTTTCTTTGTGGGGTGCCCGGGGGGCTTTTTTCGCCCCAAAAATAGTGGCGGGGTATCCGCCATGTCATCgggatagattgg                                                                                                >  1:372933/1‑240 (MQ=255)
                                                                                                                                                gCAGTGATAACCGTTTGTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTGCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGGTGTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCt                                                                                                       <  2:372933/244‑1 (MQ=255)
                                                                                                                                                      aTAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAACTTCAGGTGCTACAATAAAAAGGTTTCTGAATTTTTCCAGTTCAAGTACCCCGGAGAGGAAAGAACGTCCTAACACGCGCGTATTTTTCTTTGTGAGGTACCCGGGGGGGATTTTTCAGCACAAATATATTGCCGTAGTAACACCCATATCATCGAGATGAGACATACtaat                                                                                          >  1:316660/1‑248 (MQ=255)
                                                                                                                                                         aCCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATAcc                                                                                       <  1:433324/251‑1 (MQ=255)
                                                                                                                                                          ccGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATACCg                                                                                      <  2:734974/251‑1 (MQ=255)
                                                                                                                                                            gTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCTGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATACCggg                                                                                    <  2:651459/251‑1 (MQ=255)
                                                                                                                                                                     gtgCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTGAATTCATGGGCTTCGGTAAAACGGGTTCTGAATTTCTCCGGTTTAAGTACCCGGCCGTGGGAAGCAGTTCCTAAATCTAGCGTTTTTGTCTTTTTGGGGTCCACGGGGGCCTTTTTACGCCACCAATATAGTGCCGGGGTATCCGCCATGTTATTCAGCTGCGACTTACTGGAACCCGGACCcgcggg                                                                           >  1:238148/1‑249 (MQ=255)
                                                                                                                                                                      tgCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCcaca                                                                                                                                           <  1:138651/186‑1 (MQ=255)
                                                                                                                                                                       gCTTGCTCAGTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTGACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATACCGGGACCCGCGTgg                                                                          <  2:687559/250‑1 (MQ=255)
                                                                                                                                                                       gCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATACCGGGACCCGCGTgg                                                                          <  1:835685/250‑1 (MQ=255)
                                                                                                                                                                       gCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGGGCTCCGCTAAAAAGGGTAATGCATTATTTCTGGTCAAGTACCCGGCAGGGGAAAACAGTACCTACATTGCGCGCTTTTGTCTTTTTGGGGGTCACGGGGGCATTTTTTCGCCACAAAAACAAGGGAGGAGTATCAACAATGTCATCGAGGTGAGACTTaacgataaccggaaccgcggtgg                                                                         >  1:427729/1‑229 (MQ=255)
                                                                                                                                                                          tGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATACCGGGACCCGCGTGGTa                                                                        <  1:246329/249‑1 (MQ=255)
                                                                                                                                                                          gtcgcATTCCATCAGAATCTCTTTTTCTTCTTCTTTTACGGCGTCTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTGACTGAATTCTTCCGGTGCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGTGGGCATTGTTACGCCACAAATATAGTGCCGGAGTCTCAGCAATGTCATCGAGCTTAGACTTACTGATACCGGGACCCGAGCGGTAAt                                                                      <  2:779161/249‑1 (MQ=255)
                                                                                                                                                                             gcATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGGGCTACGATCAAGCGGTTACTGGATTATTCCGGGTCAAGTACCCGGCCGCGGAAAGCAGTTCCTAAATCGAGCGGTTTTGTCTTTGTGGTGTCCCCGGGGGCATTTTTACGCCACAAATATAGTGGCGGAGTATCCGCAATTTCATCGCGCTGAGACTTACTGATACCGGGGCCCGCGAGGTAATtct                                                                   >  1:36827/1‑251 (MQ=255)
                                                                                                                                                                                                ttCTCTTCATCATTGCCGGCGTTTGTACGGCGGTTAAATGCAGGGGCTACGATAAAGCGGTTACTGAAGTCTACCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCGAAATCGAGCGTTTGGGTCTTTGTGGTGTCCACGGGGGCATTTTTACGTCACAAATATAGTGTCGGAGTATCAGCAAGGTCATCGTGCAGAGACTTACTGATTCCGGGGCCGGCGTGGTAATTCTCATTCGAGAGTCCGtaat                                                 <  2:671108/250‑1 (MQ=255)
                                                                                                                                                                                                              tCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTTCCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATACCGGGACCCGCGTGGTAATTCTCATTCGAAATTCCGTAATAAATACCTGGctct                                   <  2:18470/250‑1 (MQ=255)
                                                                                                                                                                                                              tCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAACGGTTTCTGAATTCTTCCCGGTCCAGTACCCGGCAGGGGGAAGCAGTTCCTAAAACGAGCGTTTTTGTCTTTTTGGTGTCCCCGGGGGCCTTTTTTCGCCCCAAATATAGTGCCGGGGTATCAGCAATGTCATCGAGCTGGGACTTACTTATACCCGGACCCCCGTGGGAATTTTTATTCGAAATTCCGTTATAAAAACCTGGctctt                                  >  1:88757/1‑250 (MQ=255)
                                                                                                                                                                                                               ccGGCGTTTGTACTGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATACCGGGACCCGCGTGGTAATTCTCATTCGAAATTCCGTAATAAATACCTGGCTCTAt                                 <  1:867721/251‑1 (MQ=255)
                                                                                                                                                                                                                 ggCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATACCGGGACCCGCGTGGTAATTCTCATTCGAAATTCCGTAATAAATACCTGGCTCTATGt                               <  1:625312/251‑1 (MQ=255)
                                                                                                                                                                                                                   cGTTTGTACGGCGGATAAATTCAGGTGATACGGTAAAGCGGGTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTTAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCCTTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCCTCGAGCTGAGACTTACTGATACCGGGGACCGCGGGGTAATTCTCCTTCGAAATTCCGTAATAAATAACTGGCTCCATGtct                             >  1:684315/1‑251 (MQ=255)
                                                                                                                                                                                                                                             ggTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATACCGGGACCCGCGTGGTAATTCTCATTCGAAATTCCGTAATAAATACCTGGCTCTATGTCTTCTACGATTACGGGATCTGCGAc      >  1:246833/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 cTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATACCGGGACCCGCGTGGTAATTCTCATTCGAAATTCCGTAATAAATACCTGGCTCTATGTCTTCTACGATTACGGGATCTGCGACTTCg  >  1:847674/1‑251 (MQ=255)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                                                                            
GGACGGCACCGACACAAAATTCCTGTTTCAGGATCTTCCCAGTAAATTGATGATTCAGCGTGTCCGGCGCTTTCAACAAGCCATTGCCCCAGCGGCAAAGCCATAACGCTTTGATACATGAGTTCAATTTTCCGGCCTTCTTCCGCAGTGATAACCGTTTTTCCTGTGCTTGCGCATTCCATCAGAAACGCTTTCTCTTCTTCTTTTCCGGCGTTTGTACGGCGGTTAAATTCAGGTGCTACGATAAAGCGGTTACTGAATTCTTCCGGTTCAAGTACCCGGCAGTGGAAAGCAGTTCCTAAATCGAGCGTTTTTGTCTTTGTGGTGTCCACGGGGGCATTTTTACGCCACAAATATAGTGCCGGAGTATCAGCAATGTCATCGAGCTGAGACTTACTGATACCGGGACCCGCGTGGTAATTCTCATTCGAAATTCCGTAATAAATACCTGGCTCTATGTCTTCTACGATTACGGGATCTGCGACTTCG  >  NC_000913/1415147‑1415635

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