Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,762,3080AG32.4% 111.0 / 11.3 75G42G (GGA→GGGrhsArhs element protein RhsA
Reads supporting (aligned to +/- strand):  ref base A (17/33);  new base G (24/0);  total (42/33)
Fisher's exact test for biased strand distribution p-value = 8.87e-09
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CATGGTGATATAAATATTTTCCTAATTATTTCACTCTGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAAC  >  NC_000913/3762059‑3762509
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                              
cATGGTGATATAAATATTTTCCTAATTATTTCACTCTGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGa                                                                                                                                                                                                                <  2:317068/250‑1 (MQ=255)
  tGGTGATATAAATATTTTCCTAATTATTTCACTCTGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGGATTGTTCAGGGGTCAACCGGGGTGCGCATTGGTGCCCCCACCGGCGGGGCCAGTTCGGTGTGCCCCCGCGGCGTg                                                                                                                                                                                                             >  2:588387/1‑251 (MQ=18)
   ggTGATATAAATATTTTCCTAATTATTTCACTCTGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTg                                                                                                                                                                                                             <  2:160934/250‑1 (MQ=255)
      gATATAAATATTTTCCTAATTATTTCACTCTGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGGGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGGGTGACgg                                                                                                                                                                                                         >  1:674932/1‑250 (MQ=255)
                 tttCCTAATTATTTCACTCTGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGCCGGGGGGCGCATTGGTGCCCCCACCGGCGCGGCCTGTTCGGGGGGCCCCGGCGGGGGGACGTCCGGCCACCCg                                                                                                                                                                                              >  1:305910/1‑251 (MQ=17)
                  ttCCTAATTATTTCACTCTGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCg                                                                                                                                                                                              <  2:80460/250‑1 (MQ=255)
                  ttCCTAATTATTTCACTCTGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGGAGCCTTGTTCCGGGTTTACCCGGGGGGCGCATTGGTGCCCCCCCCGGCGCGGCCTTTTCGGTGGTGCCCGGCGGGGTGGCGTCCGGCCATCCgg                                                                                                                                                                                             >  2:538577/1‑251 (MQ=17)
                       aaTTATTTCACTCTGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGAGACGTCCGGCCATCCGGTCAAt                                                                                                                                                                                        >  1:728402/1‑251 (MQ=255)
                        attattTCACTCTGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGGGCCCCCACCGGCGTGGCCTGGTCGGGGTGCCCCCGCGGGGGGGCGTCCGGCCCTCCGGTCAATc                                                                                                                                                                                       >  1:435728/1‑251 (MQ=17)
                           attTCACTCTGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCCGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCCGCGGAGTGACGTCCCGCCATCCGGGCAATcccc                                                                                                                                                                                    >  1:101414/1‑251 (MQ=255)
                            tttacctctGATGGATATCTCACTTCAGGCTTTCTTATCAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGCGCGGAAAACCGGCAGCGCGTCAGGGCGACATGCCGCCGAATGGCGGTCGCATTTGTCCGGGTTCAGCCCGGGAGCGCCTTGGGGCCCCCCCCGGCTTGGCCTGTTCCGTGGTCCCCTCCGGGGTGCCGTCCGGcgccgcgggcaacacccc                                                                                                                                                                                   >  2:116241/6‑234 (MQ=9)
                               cACTCTGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCGTTGTTCAGGGGTCAGCCGGGGGGGGGATTGGTGCCCCCCCCCGCGTGGGCAGTTCGGGGTGCCCCGGCGGGGAGACGGCCCGCCATcacggcaaacccctccc                                                                                                                                                                                >  2:170991/1‑235 (MQ=14)
                                    tGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCCGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCCGCCGGGGTGCGCATTGGTGCCCCCCCCCGCGGGGCCTGGTCGGGGTGGCCCGGCGGGGGGAGGTACGGCCATCCGGGGCACCCCCTGCTCgggg                                                                                                                                                                          >  1:222078/1‑249 (MQ=14)
                                            atCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGCTCAGCCGGGGGGCGCATTGGTGCCCCCACCGGCGTGGCCCGTTTGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCCGCCCGGTGCAAAgggc                                                                                                                                                                 >  1:641523/1‑249 (MQ=14)
                                              ctcACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAACATTGTTCAGGGTTCAGCCGGGGGGCGCAATGGTGCCCCCCCCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGGGTGACGGCCGGCCATCCGGTCCATcccccgcccggggccaagggccc                                                                                                                                                                 >  1:553398/1‑232 (MQ=14)
                                               tcACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGAAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGCCGGGGGGCGCATTGGGGCCCCCACCGGCGTGGCCTGTTCGGGGGGCCCCCGCGGGGTGACGGCCGGCCCTCCCGGCAATCCCCTGCCCGGTGCAAAgggccct                                                                                                                                                              >  1:751784/1‑246 (MQ=14)
                                                cACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCTGTGCAAAGGTc                                                                                                                                                                 <  1:556283/247‑1 (MQ=255)
                                                     cAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCTGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCGGCTCGGTGCAAAGGTCCTTCCCGGt                                                                                                                                                        <  1:538577/251‑1 (MQ=255)
                                                           ttCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAAc                                                                                                                                                   <  1:580788/250‑1 (MQ=255)
                                                               tataAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCg                                                                                                                                                 <  2:72286/248‑1 (MQ=255)
                                                               tataAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGAc                                                                                                                                               <  1:305177/250‑1 (MQ=255)
                                                                ataAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGAATTGTTCAGGGTACAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTACGGGTGCCCCGGGGGGGGGACGGCCGGCCACCCGGTCAATCCCCTGCTCGGTGCCAAAGGACATTCCCgggagaccggcac                                                                                                                                              >  1:363994/1‑240 (MQ=14)
                                                                  aaaTCTGGAGGGTTTCGCCGGTCAGCAGTCAAAAAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGGGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGCCAATCCCCTGCTCGGGGCAAATGTCCTTCACGGTGAAACCGACATc                                                                                                                                            <  2:519655/250‑1 (MQ=14)
                                                                      cggTAGGTTTTCGCCAGTCAGCAGAAAAATATCCCGTTAAGACGAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTGAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGccc                                                                                                                                        <  2:5470/248‑1 (MQ=14)
                                                                                 cGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCCGCGGGGTGGCGTCCGGCCCTCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCCGTGAAAACGACATCCCCCCgcccggccccc                                                                                                                           >  1:717186/1‑249 (MQ=14)
                                                                                   ccTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTACCAATGGTCCGGGTTCAGCCGGGGGGCGCATTTGGGGCCCCACCGGCGGGGCCCGTTCGGGGTGCCCCCGCGGGGGGCCGTCCGGCCATCCGGTCAATCCCCTGCCCGGGGACAAGGGCCTTCCCGGTGAAACCCACACCGCCCTGCCccgcccgccg                                                                                                                         >  2:268841/1‑249 (MQ=14)
                                                                                      gTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGGCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCg                                                                                                                    <  2:499965/251‑1 (MQ=255)
                                                                                      gTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCCGCGGGGGGACGTTCGGCCCTCCGGTCAATCCCCCGCCCGGGGCAAAGGTCCTTCCCCGGGAAACCCACATCGCCCTgccccgccccccgccc                                                                                                                      >  1:377968/1‑244 (MQ=14)
                                                                                       gcagcagACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGt                                                                                                                   <  2:20496/250‑1 (MQ=255)
                                                                                           cagACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCAt                                                                                                               <  1:365445/251‑1 (MQ=255)
                                                                                                aaaTAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGTGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAAACCCCCGCTCGGGGGAAAAGGCCCTCCCGGGGAAACCGACATCGCCCCGCCCCGCCCGCTGCCGTtcctccccc                                                                                                          >  1:408793/1‑248 (MQ=14)
                                                                                                 aaTAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCtctc                                                                                                         <  2:435728/251‑1 (MQ=255)
                                                                                                  aTAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTcc                                                                                                        <  1:208939/251‑1 (MQ=255)
                                                                                                   aaaCCCGATAAAACAAGGAGGAGAAATGAGCGGAAAACCGGCAGCGCGGCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGTTGCCGTTCATCCTCTccc                                                                                                       <  2:641523/250‑1 (MQ=255)
                                                                                                     aCCCGATAAAACAAGGAGGAGAAAGGAGCGGAAAACCGGCAGCGCGGCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGc                                                                                                     <  2:553398/251‑1 (MQ=255)
                                                                                                       ccGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGGGCGCATTGGTGCCCCCACCGGCGTGGCCCGTTCGGTGGGCCCCCGCGGGGGGACGGCCGGCCATCCCGTCAATCCCCCGCTCGGTGCAAAGGGCCCTTCCGGTGAAAACCACATCCCCCTGCCCCGCCCCCCGCCGtttcttccctcccggcc                                                                                                   >  1:376916/1‑236 (MQ=14)
                                                                                                         gATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCt                                                                                                 <  1:573641/251‑1 (MQ=255)
                                                                                                          aTAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTc                                                                                                                                                                                           <  2:217794/165‑1 (MQ=255)
                                                                                                          aTAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGGGCGCATTGGTGCCCCCCCCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGGGGGGCGGCCGGCCATCCgggc                                                                                                                                                                                           >  1:217794/1‑163 (MQ=14)
                                                                                                           aaaaaCAAGGATGAGAAAGTAGCGGACAACCGGCAGCGCGTCAGGGCGAAGTGACGCAGTATGGCGGTAGCATTGTTCAGGGTACAGCCGGGGTGCGCATTGGTGCCCCTACCGGCGTGGCCTTTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGGCAATCCCCTGCTCGGTGCTAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCAGGCCCGCTGCCGTTCATCCTCTCCCGCACCTAc                                                                                               <  2:102179/250‑1 (MQ=14)
                                                                                                                  aGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGGGCGCATTGGTGCCCCCACCGGCGTGGCCCGTTCGGGGTGGCCCGGCGGGGGGACGTCCGGCCCTCCGGTCAATCCCCCGCCCGGGGCCAAGGGCCTTCCCGGGGAAACCCACCTCGCCCTGCCCCGGCCCCCGCCGTTCATCCCCTCCCCCACCTCCACCAGtt                                                                                        >  1:111948/1‑251 (MQ=14)
                                                                                                                  aGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCAATGTTCCGGGGTCAGCCGGGGGGCGCATTGGTGCCCCCCCCGGCGTGGGCTGTTTGGGGTGCCCCCGCGGGGGGACGTCCCGCCACCCGGGCAATCCCCCGCTCGGGGGCAAGGGCCTTCCCGGGGAAAACGACCTCCCCCCCCCCGGCCCCCCGCCGTTCCtccctccccccccccccaacagct                                                                                       >  2:323196/1‑235 (MQ=9)
                                                                                                                   ggATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGGGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGAGTGACGTCCGGGCATCCGGGCCATCCCCTGCTCGGTGCAAAGGGCCTTCCCGGTGAAACCGACATCGCCCCGCCCGGCCCGCTGCCGGTCATCCTCCCCCGCACCCAcagcaggt                                                                                        >  1:409917/1‑248 (MQ=14)
                                                                                                                    gATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTTTTCCGGGTTCAGCCCGGGGGCGCATTGGTGCCCCCACCGGCGTGGCCCGTTCGGGGTGCCCCGGGGGGGTGGCGCCCGGCCCTCCCGTCAATCCCCTGCTCCGGGCCAAAGGCCCTCCCGGGGAAACCGACATCCCCCCCCCCGGCCCGCGCCCCTTCATCCTCTCCCCCAcccacctccctttc                                                                                      >  2:552481/1‑239 (MQ=9)
                                                                                                                     aTGAGAAAGGAGCGGAAAAACGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGGGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGGCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTAcc                                                                                     <  2:30431/251‑1 (MQ=255)
                                                                                                                        agaAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGa                                                                                  <  1:70382/251‑1 (MQ=38)
                                                                                                                         gaAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTTGGGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGAGGGACGTCCGGCCATCCGGGCAATCCCCTTCTCGGTGCAAAAGTCCTTCCCGGTGAAACCCACATCCCCCTTCCCCGCCCGCTGCCGTTCATCCTCTcccccccccacagccggtacccggc                                                                                 >  1:215591/1‑229 (MQ=11)
                                                                                                                          aaaTGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGGTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGGGGGACGTCCGGCCCTCCGGGCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCCGCCGCTCATCCTCTCCCCCCCCTACAGCCGTTACCGGAc                                                                                 >  1:757492/1‑250 (MQ=11)
                                                                                                                             tGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCCGGGGTCAGCCGGGGTGCGCATTGGTGCCCCCCCCCGCGGGGCCTGGTTGGGGTGCCCCGGCGGGGGGACGTACGGCCATCCCGTTAATCCCCCGCTCCGGGCAAAGGTCCTTCCCCGTGAAACCGACATCGCCCTGCCCCGCCCGCCGCCCGTTCTCCTCTCCCCCCCCTACAGCAGTTACCgggaaca                                                                              >  2:694122/1‑245 (MQ=9)
                                                                                                                                gCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTAGGGCGGTAGCATTGTTCAGGGGTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGCCGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGATACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACg                                                                          <  2:561719/251‑1 (MQ=255)
                                                                                                                                 cGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGc                                                                         <  2:717186/251‑1 (MQ=37)
                                                                                                                                    aaaaCCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCcg                                                                      <  1:323196/251‑1 (MQ=37)
                                                                                                                                      aaCCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGGGGGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGGGCAAAAGTCCTTCCCGGGGAAACCGCCATCGCCCTGCCCGGCCCGCCGCCGTTCATCCTCTCCCCCCCCTACAGCCGTTACCGGGCAAAAACGCCcgcc                                                                    >  1:414442/1‑250 (MQ=11)
                                                                                                                                       aCCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCCTTGTTCAAGGTTCAGCCGGGGTGGGCATTTGTGCCCCCACCGGCGTGGCCCGTTCGGTGGGCCCCGGCGGAGTGACGTCCGGCCCTCCCGTCAATCCCCCGCTCGGGGCCAAGGGGCCTCCCCGTGAAACCGACATCGCCCTGCCCGGCCCCCCGCCGCTCCTCCCCTCCCCCACCTACACCAGTTACCCGGCAAAAAAGCCCgcgc                                                                   >  1:765179/1‑251 (MQ=9)
                                                                                                                                          ggCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCCGCGGGGGGACGGCCGGCCATCCGGTCAATCCCCTGCTCGGGGCAAAGGTCCTTCCCGGGGAAACCGACATCGCCCTGCCCGGCCCGCCGCCGTTCATCCTCTCCCGCACCCACAGCAGGTACCGGACAAAAACGCCCGCCCCgg                                                                >  1:693156/1‑251 (MQ=11)
                                                                                                                                           gCAGCGCGTCAGGGCGACATGACGTAGTATGGCGGTAGCATGGTGCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTGCGGGGTGCCCCGGCGGATTGACGTCCGGCCATCCGGGCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCAGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGt                                                               <  2:440627/251‑1 (MQ=11)
                                                                                                                                            cAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCCCCCGGGGGGCGCCTTTGTTCCCCCACCGGCGTGGCCTGTTCGGTGTTCCCCGGCGGAGGGACGGCCGGCCCTCCCGTCAAACCCCCTCTCGGTGCAAAGGTCCTTCCCGGGGAAACCGACATCGCCCCTCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAACAGTTACCGGACAAAAACGCCCGCGCCGGTg                                                              >  1:364454/1‑251 (MQ=11)
                                                                                                                                             aGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGTGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTgg                                                             <  2:728402/251‑1 (MQ=25)
                                                                                                                                             aGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCCGTGCAAAGGGCCTTCCCGGTGAAACCGACATCGCCCTGCCCCGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCCCCggggg                                                             >  1:670674/1‑248 (MQ=255)
                                                                                                                                                         cGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAACCGGGGTGCGCCTTGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAAACCCCTGCTCCGTGCAAAGGGCCTTCCCGGCGAAACCGACCTCGCCCTGCCCCGCCCCCTGCCGGTCCTCCTCTCCCCCACCCACAGCAGTTACCGGACACAAACGCCCGCCCCGGTGGGGAGCCTCGGcc                                                 >  1:56667/1‑251 (MQ=11)
                                                                                                                                                           aCATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGGGTGCCCCGGCGGGGTGACGTCCGGCCCTCCGGTCAATCCCCTGCTCGGGGCAAAGGGCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCCTCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGGGGGGAGCCTCGGcccc                                               >  1:500254/1‑251 (MQ=11)
                                                                                                                                                               gACGCAGTATGGCGTTAGCATTGTTCAGGGTGCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTGCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCt                                           <  2:674932/251‑1 (MQ=18)
                                                                                                                                                               gACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCCGGGTGCGCATTGGAGGCCCCACCGGCGAGGCCTGTTCGGTGTGCCCCGGCGGAGGGACGTCCGGCCCTCCGGTCAATCCCCTGCTCGGGGCCAAGGGCCTTCCCGGGGAAACCGAACTCGCCCTGCCCCGCCCGCTGCCGTTCATCCTCTCCCGCACCCACAGCCGTTACCGGAAAAAAACCCCCGCGCCGGGGGGGGGGCTCGGcacccgct                                           >  1:436684/1‑244 (MQ=11)
                                                                                                                                                                  gCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGGTCGGGGGGGCCCCGCGGGGAGACGGCCGGCCCTCCGGGCAATCCCCCGCTCGGGGCCAAGGTCCTTCCCGGGGAAACCGACATTGCCCCGCCCGGCCCGCTGGCCGTCCTCCCCCCCCGCACCCACCACCGGTACCCGCCACAAACGCACGCCCCGGTGGGGAGCCCCGGCCCCCGCTggg                                        >  2:540059/1‑250 (MQ=2)
                                                                                                                                                                  gCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGGGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGGCAATCCCCCGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCCCCCTGCCCCGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCCGTTACCGGACAAAAACGCCCGCGCCGGTGGGGGGCCTCGGCCCCGGCTggg                                        >  1:165729/1‑250 (MQ=11)
                                                                                                                                                                          gCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGCGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCt                                <  2:414442/251‑1 (MQ=18)
                                                                                                                                                                           cGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTg                               <  1:750842/251‑1 (MQ=33)
                                                                                                                                                                               aGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGa                           <  1:540059/251‑1 (MQ=33)
                                                                                                                                                                                  aTTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGGTCATCCTCTCCCGCCCCTACCGCAGTTCCCGGACCAAAACGCCCCCGCCGGCGGGGAGCCTCGGCCCCGGCTGGGAAAACCCTGCgggtt                         >  2:275840/1‑247 (MQ=11)
                                                                                                                                                                                      ttCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGc                    >  1:1823/1‑251 (MQ=33)
                                                                                                                                                                                       tcgtgtttCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCTTTTGCGGTGGGCCCCGGCGGAGCGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCAGCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGc                    <  2:765179/244‑1 (MQ=11)
                                                                                                                                                                                       gcAGGGTTCAGCCGGGGTGTGCATTGGTGCCCCCCCCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGc                    <  2:446816/249‑1 (MQ=18)
                                                                                                                                                                                              tCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGc             <  1:694122/250‑1 (MQ=33)
                                                                                                                                                                                                 gCCGGGGGGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGGCGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACCTCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAACACGCCCCCGCCGGCGGGGAGCCTCCGCCCCGGCTGGAAAAAGCCAGCGGatcaccggtttcagccacg         >  2:664028/1‑234 (MQ=11)
                                                                                                                                                                                                     gggTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCCTACAGCCGCGCGAt     >  2:101227/1‑251 (MQ=18)
                                                                                                                                                                                                         gcgcATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAac  >  1:489425/1‑250 (MQ=31)
                                                                                                                                                                                                                                                           |                                                                                                                                                                                                              
CATGGTGATATAAATATTTTCCTAATTATTTCACTCTGATGGATATCTCACTTCAGGCTTTCTTATAAATCTGTAGGGTTTCGCCTGTCAGCAGACAAATAACCCGATAAAACAAGGATGAGAAATGAGCGGAAAACCGGCAGCGCGTCAGGGCGACATGACGCAGTATGGCGGTAGCATTGTTCAGGGTTCAGCCGGGGTGCGCATTGGTGCCCCCACCGGCGTGGCCTGTTCGGTGTGCCCCGGCGGAGTGACGTCCGGCCATCCGGTCAATCCCCTGCTCGGTGCAAAGGTCCTTCCCGGTGAAACCGACATCGCCCTGCCCGGCCCGCTGCCGTTCATCCTCTCCCGCACCTACAGCAGTTACCGGACAAAAACGCCCGCGCCGGTGGGGAGCCTCGGCCCCGGCTGGAAAATGCCTGCGGATATCCGCTTACAGCTGCGCGATAAC  >  NC_000913/3762059‑3762509

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