Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH1325,4730AG53.6% 125.9 / inf 628G8G (GGT→GGC) rtcRFis family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (133/157);  new base G (130/205);  total (263/365)
Fisher's exact test for biased strand distribution p-value = 8.79e-02
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.71e-03

Only 100 of 644 total aligned reads displayed.

GTTTCAAATAGCGAGCGCGAGCGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAACCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTAAAACAATGAATTAAAAAATTAAGCTGCTTATTT  >  DH1/325324‑325623
                                                                                                                                                     |                                                                                                                                                      
   tCAAATAGCGAGCGCGAGCGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCa                                                                                                                                                     <  5:1328999/150‑1 (MQ=255)
     aaaTAGCGAGCGCGAGCGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCaaa                                                                                                                                                   >  5:368201/1‑150 (MQ=255)
     aaaTAGCGAGCGCGAGCGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAACCaaaa                                                                                                                                                  <  2:1157503/151‑1 (MQ=255)
      aaTAGCGAGCGCGAGCGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCaaa                                                                                                                                                   >  3:1986907/1‑149 (MQ=255)
       aTAGCGAGCGCGAGCGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGc                                                                                                                                                <  5:1532952/151‑1 (MQ=255)
        tAGCGAGCGCGAGCGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAACCAAAAGcc                                                                                                                                               >  6:1950982/1‑151 (MQ=255)
          gcgagcgcgagcGGGCGTCGTGCAACAATTCCATTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAACCAAAAGCCAc                                                                                                                                             <  1:1727257/151‑1 (MQ=255)
           cgagcgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACt                                                                                                                                            >  3:1091777/1‑151 (MQ=255)
            gagcgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTg                                                                                                                                           >  7:1975160/1‑151 (MQ=255)
            gagcgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTg                                                                                                                                           >  1:355129/1‑151 (MQ=255)
             agcgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGt                                                                                                                                          >  1:1401058/1‑151 (MQ=255)
              gcgcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGtt                                                                                                                                         <  3:387092/151‑1 (MQ=255)
                gcgagcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGtttt                                                                                                                                       <  7:1757022/151‑1 (MQ=255)
                 cgagcgGGCGTCGTGCAACAAATCAAGTAGATAGATGACCAACCATTCTTGCTGGCATAAAAATATTATCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGtttt                                                                                                                                       >  8:1936214/1‑150 (MQ=255)
                    gcgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGc                                                                                                                                   >  4:456643/1‑151 (MQ=255)
                     cgGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAACCAAAAGCCACTGTTTTACGc                                                                                                                                   >  8:1209035/1‑150 (MQ=255)
                      gggCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCAt                                                                                                                                 >  5:2031326/1‑151 (MQ=255)
                       ggCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAACCAAAAGCCACTGTTTTACGCATc                                                                                                                                >  1:923624/1‑151 (MQ=255)
                          gtcgtGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGGCGCCATTTTGACCAGCGCTGACTGCAGCGCCCGGCATAATCCACTACGGTTCCGACAAAGCCAAAATCCACTGTTTTACGCATCtt                                                                                                                              >  6:32206/1‑150 (MQ=255)
                            cgtGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAg                                                                                                                            <  1:123663/150‑1 (MQ=255)
                             gtGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAACCAAAAGCCACTGTTTTAc                                                                                                                                     >  8:642015/1‑140 (MQ=255)
                              tGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAACCAAAAGCCACTGTTTTACGCATCTTAGat                                                                                                                          >  6:1521755/1‑150 (MQ=255)
                               gCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGatat                                                                                                                        >  3:844182/1‑151 (MQ=255)
                                caacaaTTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATc                                                                                                                       >  2:1070272/1‑151 (MQ=255)
                                 aacaaTTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATcc                                                                                                                      >  2:496138/1‑151 (MQ=255)
                                 aacaaTTCCAGTCGATCGATGACCAACGATTCTTGCTGGCAAAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATcc                                                                                                                      <  4:1091767/151‑1 (MQ=255)
                                  acaaTTCCAGTCGATCGATGATCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAACCAAAAGCCACTGTTTTACGCATCTTAGATATCCt                                                                                                                     >  3:157101/1‑151 (MQ=255)
                                  acaaTTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCt                                                                                                                     <  4:1838848/151‑1 (MQ=255)
                                  acaaTTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAACCAAAAGCCACTGTTTTACGCATCTTAGATATCCt                                                                                                                     <  8:1062578/151‑1 (MQ=255)
                                   caaTTCCAGTCGATCGATGACCAACGATTCTTGCTGTCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCtt                                                                                                                    <  7:638296/151‑1 (MQ=255)
                                   caaTTACAGTCGATCGATGACCAAAGATTCTTGCTGGCATAAACAGAGTGTCGGACGGCATTATGACCAGCGCTGACTGCCGCGGCAGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTCTGTTACGCATCTTAGATATACtt                                                                                                                    >  2:1549629/1‑151 (MQ=255)
                                     aTTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTtat                                                                                                                  <  7:141088/151‑1 (MQ=255)
                                      ttCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAACCAAAAGCCACTGTTTTACGCATCTTAGATATCCTtat                                                                                                                  >  4:1621524/1‑150 (MQ=255)
                                          aGTCGATCGATGACAAACGATTCTTGCTGGCTTAAACAGAGTGTCGGACGCCATTTTGACCAGAGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAACACCAAAAGCCACTGTTTTACGCATCTTAGATATACTTATaaaac                                                                                                             >  7:1687433/1‑150 (MQ=255)
                                           gTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATaaa                                                                                                               >  5:890006/1‑148 (MQ=255)
                                             cgattgatGACGAACGATTCTTGCTGGCATAAACAGAGTGTCTGTCGCCATTTTGAGCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAACCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGAc                                                                                                           <  4:282648/150‑1 (MQ=255)
                                              gatcgatGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAACCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACg                                                                                                          >  8:2019961/1‑150 (MQ=255)
                                                tcgatGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGat                                                                                                        >  7:317194/1‑150 (MQ=255)
                                                  gatGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATAtt                                                                                                     >  1:1085723/1‑151 (MQ=255)
                                                  gatGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATAtt                                                                                                     <  1:100171/151‑1 (MQ=255)
                                                  gatGACCAACGAATCTTGCTGGCATAAACAGACTGTCGCACGCCATTTTGACCAGCGCTGACTGCCCCGCCCGGCATAATCCAGAACGGTACCGCCAAAACCAAAAGCCACTCTTTTACGCATCTTACATATCCTAATAAAAGACCATAtt                                                                                                     >  8:767911/1‑151 (MQ=255)
                                                   atGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATAttt                                                                                                    <  8:1430112/151‑1 (MQ=255)
                                                       ccAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCAtt                                                                                                >  5:518313/1‑151 (MQ=255)
                                                        cAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCAtt                                                                                                >  4:929700/1‑150 (MQ=255)
                                                           cGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAACAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCg                                                                                            <  2:322441/151‑1 (MQ=255)
                                                           cGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCg                                                                                            <  3:43348/151‑1 (MQ=255)
                                                              ttCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAACCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGt                                                                                          >  6:1161828/1‑150 (MQ=255)
                                                               tCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAACCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTa                                                                                         >  4:1232921/1‑150 (MQ=255)
                                                                 ttGCTGGCATAAACAGAGTCTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAACCAAAAGCCACTCTTTTACGCATATTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAAt                                                                                      >  2:751281/1‑151 (MQ=255)
                                                                  tGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTaaa                                                                                       <  8:1813961/149‑1 (MQ=255)
                                                                  tGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTaaa                                                                                       <  1:1549963/149‑1 (MQ=255)
                                                                  tGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAAtt                                                                                     >  7:1268974/1‑151 (MQ=255)
                                                                  tGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAACCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAAtt                                                                                     >  3:941456/1‑151 (MQ=255)
                                                                   gCTGGCATAAACAGAGTGTGGGCCGGCATTTTGACCAGCGCTGACTGCCGGGCCCGGCATAATCCAGAACGGTACCGACAAAGCCAAACGCCAGTGTTTTACGCATCTTATATATCCTTATAAAACACGATATTTCATTATCGGTAAATTg                                                                                    <  1:589686/151‑1 (MQ=255)
                                                                   gCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTg                                                                                    <  2:69365/151‑1 (MQ=255)
                                                                   gCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTg                                                                                    >  2:900528/1‑151 (MQ=255)
                                                                     tGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAACCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGc                                                                                  >  5:1468277/1‑151 (MQ=255)
                                                                      ggCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTgg                                                                                   <  1:22819/149‑1 (MQ=255)
                                                                       gCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAACCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCtt                                                                                >  5:1073292/1‑151 (MQ=255)
                                                                           aaaCAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAACCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGc                                                                             <  7:822596/150‑1 (MQ=255)
                                                                            aaCAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAACCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCgg                                                                           <  7:482602/151‑1 (MQ=255)
                                                                                gagTGTCGGACGCCATTTTGACCAGCGGTGACTGCCGCGCCCGGGATAATCCAGTACGGTACCGAGAAAAGCAAAAGCCTCTGTTTTACGGATCTTAGATATACTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGaaa                                                                       <  3:1191923/151‑1 (MQ=255)
                                                                                cagTGTCGGACGCCATTTTGACCAGCGCGGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCACAAGCCACTGATTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGaaa                                                                       <  6:1755819/150‑1 (MQ=255)
                                                                                       ggACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCgg                                                                <  8:500045/151‑1 (MQ=255)
                                                                                       ggACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAACCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCgg                                                                <  8:920768/151‑1 (MQ=255)
                                                                                          cGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAACCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTa                                                             <  1:1494134/151‑1 (MQ=255)
                                                                                            ccATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAACCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGCCGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACg                                                           >  6:1396090/1‑151 (MQ=255)
                                                                                            ccATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAACCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACg                                                           >  3:1380294/1‑151 (MQ=255)
                                                                                                tttGACCAGCGCTGACTGCCGCGCCCGGCATAATCTAGTACGGTACCGACAAAACCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGAc                                                        <  2:1402556/150‑1 (MQ=255)
                                                                                                tttGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGAc                                                        <  5:1782739/150‑1 (MQ=255)
                                                                                                      cAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAACCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGa                                                 <  7:74012/151‑1 (MQ=255)
                                                                                                        gcgcTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAACCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGata                                               <  6:886052/151‑1 (MQ=255)
                                                                                                         cgcTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGatat                                              <  3:299727/151‑1 (MQ=255)
                                                                                                         cgcTGACTGCCGCGCCCGGCATAATCCAGTACGCTACCGACAAAACCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGatat                                              <  6:571097/151‑1 (MQ=255)
                                                                                                            tGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAACCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAAt                                           >  5:1254027/1‑151 (MQ=255)
                                                                                                             gACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATaa                                            <  3:1659571/149‑1 (MQ=255)
                                                                                                             gACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACCGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATc                                          >  7:1549850/1‑151 (MQ=255)
                                                                                                                tGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTACGGTTACGGACGTAAAGATATAATCttt                                       >  6:1002830/1‑151 (MQ=255)
                                                                                                                 gCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCtttt                                      <  5:1891696/151‑1 (MQ=255)
                                                                                                                 gCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCttt                                       <  3:548514/150‑1 (MQ=255)
                                                                                                                   cgcgCCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCtttttt                                    <  3:1827845/151‑1 (MQ=255)
                                                                                                                     cgcCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTaa                                  <  4:368517/151‑1 (MQ=255)
                                                                                                                     cgcCCGGCATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTaa                                  <  2:1871318/151‑1 (MQ=255)
                                                                                                                        cctgaaTAAACCATTAAGGTACCGAAAAAACCAAAAGCCACTGTTTTACCCATCTTAGATATCCTTATAAAAGACGAAATTTCATTATCTGTAAATTTGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTaaaaa                               <  5:1672080/146‑2 (MQ=255)
                                                                                                                        ccGGCATAATCCAGTACGGTACCGACAAAACCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTAAAAc                               >  5:705372/1‑151 (MQ=255)
                                                                                                                          ggCATAATCCAGTACGGTACCGACAAAACCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTAAAACaa                             <  2:1757823/151‑1 (MQ=255)
                                                                                                                            cATAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTAAAACAATg                           <  2:540834/151‑1 (MQ=255)
                                                                                                                              tAATCCAGTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTAAAACAATGa                          >  2:966247/1‑150 (MQ=255)
                                                                                                                                  ccAGTACGGTACCGACAAAACCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTAAAACAATGAATTaa                     <  4:1379730/151‑1 (MQ=255)
                                                                                                                                     gTACGGTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTAAAACAATGAATTaaaaa                  <  8:336327/151‑1 (MQ=255)
                                                                                                                                       aCGGTACCGACAAAACCAAAATCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTAAAACAATGAATTAAAAAAt                <  8:15979/151‑1 (MQ=255)
                                                                                                                                       aCGGTACCGACAAAACCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTAAAACAATGAATTAAAAAAt                >  3:686389/1‑151 (MQ=255)
                                                                                                                                        cGGTACCGACAAAACCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTAAAACAATGAATTAAAAAAtt               <  1:823048/151‑1 (MQ=255)
                                                                                                                                         ggTACCGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTAAAACAATGAATTAAAAAATTa              <  6:367514/151‑1 (MQ=255)
                                                                                                                                            aCCGACAAAACCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTAAAACAATGAATTAAAAAATTa              >  8:1376774/1‑148 (MQ=255)
                                                                                                                                              cGACAAAGCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTAAAACAATGAATTAAAAAATTAAgctg         <  4:849219/151‑1 (MQ=255)
                                                                                                                                                aCAAAACCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTAAAACAATGAATTAAAAAATTAAgctgct       >  7:967190/1‑151 (MQ=255)
                                                                                                                                                  aaagccaaaagccaCTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTAAAACAATGAATTAAAAAATTAAGCTGCTTa     <  4:631739/151‑1 (MQ=255)
                                                                                                                                                  aaaaCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTAAAACAATGAATTAAAAAATTAAGCTGCTTa     <  6:709406/151‑1 (MQ=255)
                                                                                                                                                  aaaaCCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATAGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTAAAACAATGAATTAAAAAATTAAGCTGCtt      >  7:1056831/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                      
GTTTCAAATAGCGAGCGCGAGCGGGCGTCGTGCAACAATTCCAGTCGATCGATGACCAACGATTCTTGCTGGCATAAACAGAGTGTCGGACGCCATTTTGACCAGCGCTGACTGCCGCGCCCGGCATAATCCAGTACGGTACCGACAAAACCAAAAGCCACTGTTTTACGCATCTTAGATATCCTTATAAAAGACGATATTTCATTATCGGTAAATTGGCTTTGCGGGAAATCTCCGGTTACGGACGTAAAGATATAATCTTTTTTAAAACAATGAATTAAAAAATTAAGCTGCTTATTT  >  DH1/325324‑325623

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