Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*1428,8820AG67.9% 53.3 / 33.0 56intergenic (‑270/+522)unknown/unknownhypothetical protein/LSU rRNA ## 23S rRNA, large subunit ribosomal RNA
Reads supporting (aligned to +/- strand):  ref base A (11/7);  new base G (11/27);  total (22/34)
Fisher's exact test for biased strand distribution p-value = 3.88e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

AATACCAAATCAGCGACCTTATTGCCCATTATTCGAGAGCAATTCAAGCCTGATAGCACCGTTTATCCGTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGAGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTTTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACAC  >  1/428643‑429046
                                                                                                                                                                                                                                               |                                                                                                                                                                    
ggtccacatatcctttcttcgataaattatcctgaatcatgtcttccttttcctgattccttatatCCGTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTAt                                                                                                                                                           >  1:389451‑M2/67‑251 (MQ=255)
          tcctttcttcgataaattatcctgaatcatgtcttccttttcctgattccttatatCCGTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCg                                                                                                                                                   >  1:776683‑M2/57‑249 (MQ=255)
                   cgataaattatcctgaatcatgtcttccttttcctgattccttatatCCGTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCtt                                                                                                                                        >  1:887480‑M2/48‑251 (MQ=255)
                         attatcctgaatcatgtcttccttttcctgattccttatatCCGTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTc                                                                                                                                  <  1:803722‑M2/210‑1 (MQ=255)
                            atcctgaatcatgtcttccttttcctgattccttatatCCGTTCTAAGCCTTCCGTAAATGGCTTGTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTa                                                                                                                                               <  1:454353‑M2/197‑1 (MQ=255)
                             tcctgaatcatgtcttccttttcctgattccttatatCCGTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATa                                                                                                                              >  1:344126‑M2/38‑251 (MQ=255)
                                 gaatcatgtcttccttttcctgattccttatatCCGTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGct                                                                                                                                           <  1:637616‑M2/201‑1 (MQ=255)
                                            tccttttcctgattccttatatCCGTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGt                                                                                                               <  1:200055‑M2/229‑1 (MQ=255)
                                              cttttcctgattccttatatCCGTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGgttg                                                                                                             <  1:540410‑M2/231‑2 (MQ=255)
                                                   cctgattccttatatCCGTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAagctctttgattccgtacgctatgcgtatggaatcaaag                                                                                                        <  1:416073‑M2/236‑20 (MQ=255)
                                                          acctatatCCGTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGAGGAtttttt                                                                                                                                                             <  1:147992‑M2/183‑1 (MQ=255)
                                                          acctatatCCGTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGAGGAtttt                                                                                                                                                               <  1:796573‑M2/181‑1 (MQ=255)
                                                          acctatatCCGTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGAGGAtt                                                                                                                                                                 <  1:167894‑M2/179‑1 (MQ=255)
                                                          acctatatCCGTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGAGGATTTTTTATGt                                                                                                                                                         >  1:580187‑M2/9‑195 (MQ=255)
                                                          acctatatCCGTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGAGGATTTTTTATGt                                                                                                                                                         <  1:77317‑M2/187‑1 (MQ=255)
                                                          acctatatCCGTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGAGGATTTTTTATGt                                                                                                                                                         >  1:520300‑M2/9‑195 (MQ=255)
                                                          acctatatCCGTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGAGGATTTTTTATGTTg                                                                                                                                                       >  1:136013‑M2/9‑197 (MQ=255)
                                                          acctatatCCGTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGAGGATTTTTTATGTTGACCGCTTAAGctc                                                                                                                                          >  1:358920‑M2/9‑210 (MQ=255)
                                                          acctatatCCGTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGAGGATTTTTTATGTTGACCGCTTAAGct                                                                                                                                           >  1:349061‑M2/9‑209 (MQ=255)
                                                          acctatatCCGTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGAGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATa                                                                                                                              >  1:673076‑M2/9‑222 (MQ=255)
                                                           cctatatCCGTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGAGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTacggttttgttt                                                                                                        <  1:551573‑M2/236‑5 (MQ=255)
                                                            ttatatCCGTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACaa                                                                                               <  1:398054‑M2/245‑1 (MQ=255)
                                                               tatCCGTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGAGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCAAATGCATAGc                                                                                                                      >  1:740698‑M2/4‑225 (MQ=255)
                                                                    gTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAg                                                                                        >  1:481458/1‑250 (MQ=255)
                                                                    gTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGAGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCg                                                                                                                     >  1:148355/1‑221 (MQ=32)
                                                                    gTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGAGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGttttgttt                                                                                                        <  1:171621/234‑5 (MQ=11)
                                                                     ttCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCAt                                                                                                                         >  1:402851/1‑216 (MQ=17)
                                                                     ttCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTta                                                                                                          <  1:555390/231‑1 (MQ=17)
                                                                     ttCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCt                                                                                      <  1:683253/251‑1 (MQ=255)
                                                                              ttCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACaaaa                                                                             >  1:423367/1‑251 (MQ=255)
                                                                                             ttcttcTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGAGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGc                                                                                                                           >  1:733964/1‑190 (MQ=31)
                                                                                               cttctATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACg                                                                                                                            <  1:823261/187‑1 (MQ=14)
                                                                                               cttctATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTca                                                             >  1:697884/1‑250 (MQ=255)
                                                                                                  ctATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATc                                                         <  1:383538/251‑1 (MQ=255)
                                                                                                  ctATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATc                                                         <  1:431610/251‑1 (MQ=255)
                                                                                                     tatTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATCAgg                                                      <  1:750447/251‑1 (MQ=255)
                                                                                                              gTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGAGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGc                                                                                                                      <  1:671009/178‑1 (MQ=32)
                                                                                                                  aGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAgg                                         >  1:605527/1‑251 (MQ=255)
                                                                                                                    aCGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCTTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGc                                        <  1:447879/250‑1 (MQ=255)
                                                                                                                    aCGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGc                                        <  1:707975/250‑1 (MQ=255)
                                                                                                                     cGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATc                                                         <  1:839331/232‑1 (MQ=255)
                                                                                                                       aTGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCAt                                                                                                                         <  1:889998/166‑1 (MQ=17)
                                                                                                                       aTGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATc                                                         <  1:864807/230‑1 (MQ=255)
                                                                                                                         gtgATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGt                                                                                                               <  1:406391/174‑1 (MQ=17)
                                                                                                                                ttACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGAGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGt                                                                                                               >  1:605431/1‑167 (MQ=32)
                                                                                                                                     gTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGt                                                                                                               <  1:11793/162‑1 (MQ=17)
                                                                                                                                         cATTTGGTATGTGATAGACTTGAATCGTCAGTGGCTTTACTTACCATTGACTTATCCGCCTGTATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGGCAACGTCTGGCTATAACAAAAAGCCCCGATGCCCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGGCCTACTCTcaca                  <  1:300667/251‑1 (MQ=14)
                                                                                                                                               gTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGa            >  1:320201/1‑251 (MQ=37)
                                                                                                                                               gTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGAGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCAt                                                                                                                         >  1:286565/1‑142 (MQ=32)
                                                                                                                                                   gtgATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGc                                                                                       >  1:214452/1‑172 (MQ=255)
                                                                                                                                                   gtgATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGAcc                            <  1:255498/231‑1 (MQ=25)
                                                                                                                                                   gtgATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATgggg             <  1:101300/246‑1 (MQ=25)
                                                                                                                                                    tgATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGccc       <  1:739401/251‑1 (MQ=21)
                                                                                                                                                    tgATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCCCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGccc       <  1:677476/251‑1 (MQ=21)
                                                                                                                                                          aCTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCacac  <  1:854493/250‑1 (MQ=17)
                                                                                                                                                          aCTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGGGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTGTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCacac  <  1:856424/250‑1 (MQ=17)
                                                                                                                                                                                                                                               |                                                                                                                                                                    
AATACCAAATCAGCGACCTTATTGCCCATTATTCGAGAGCAATTCAAGCCTGATAGCACCGTTTATCCGTTCTAAGCCTTCCGTAAATGGCTTTTCTTCTATATTTCGGGGTGAAGACGATGTGATACTTACAGTTCCATTTGGTATGTGATAGACTTGAATCGTCATTGGCTTTACTTACCATTGACTTATCCTCCTATATCTTGAATTTGGGTTGTCAGCCTTATTCTTATATAGTGAGGATTTTTTATGTTGACCGCTTAAGCTCTTTGATTCCATACGCATAGCGTACGGTTTTTATAGCTAGACAACGTCTAGCTATAACAAAAAGCCCCGATGCTCTCACATCAGGGCTCCCTATCAGGCTTGGCGGTGACCTACTCTCACATGGGGATGCCCCACAC  >  1/428643‑429046

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