Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH11,165,4440AC22.7% 646.0 / 48.2 279T71P (ACG→CCG) yfiFmethyltransferase
Reads supporting (aligned to +/- strand):  ref base A (43/172);  new base C (63/0);  total (106/173)
Fisher's exact test for biased strand distribution p-value = 4.10e-34
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 280 total aligned reads displayed.

ATGATTCTGATAAACGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGCACTTTTCC  >  DH1/1165295‑1165594
                                                                                                                                                       |                                                                                                                                                       
 tgatTCTGATAAACGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGGGCAcc                                                                                                                                                       >  1:650240/1‑150 (MQ=255)
   atTCTGATAAACGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGt                                                                                                                                                     <  3:62990/151‑1 (MQ=255)
    ttCTGATAAACGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGGCGCCCGGt                                                                                                                                                     >  4:898859/1‑151 (MQ=255)
     tCTGATAAACGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGtt                                                                                                                                                    <  6:454987/150‑1 (MQ=255)
        gATAAACGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGttt                                                                                                                                                   <  7:1353699/148‑1 (MQ=255)
          tAAACGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGGGCCCCGGTTGCcgc                                                                                                                                              >  4:1414889/1‑151 (MQ=255)
                   ccACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGt                                                                                                                                     <  8:886960/151‑1 (MQ=255)
                    cacaACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGtg                                                                                                                                    <  3:69301/151‑1 (MQ=255)
                     acaACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTAACGCGCGCCGCGtgc                                                                                                                                   >  6:611770/1‑150 (MQ=255)
                       aaCCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGtgatg                                                                                                                                 <  3:128017/151‑1 (MQ=255)
                       aaCCCGCGTACCGGGAAAGGGGGCGGGCCTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCCGACGCGTTAGCGCGGGCCgggggggg                                                                                                                                 >  2:505846/1‑146 (MQ=255)
                        aCCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCGGTTTCCCGCGCGCCCGGtgatg                                                                                                                                 >  4:346559/1‑150 (MQ=255)
                         cccGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATgag                                                                                                                               <  1:211453/151‑1 (MQ=255)
                          ccGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCCCACCGGTTCCCCCGCGCCGGGGGATgagc                                                                                                                              >  4:1124968/1‑150 (MQ=255)
                           cgcgTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGGGCCCGCGGTGCCGCCCGCCggggggtgggcc                                                                                                                             >  5:395169/1‑142 (MQ=255)
                           cgcgTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCGGGTTGCCGCGCGCCggggggtcagcc                                                                                                                             >  6:1407092/1‑142 (MQ=255)
                              gTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACCGGTTCCCCCGCCGCGGGtgctgccccgcc                                                                                                                          >  5:107535/1‑144 (MQ=255)
                                aCCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGGCCGGTTTGCCCCGCGCCGGGGGCTGCGCCGCCgg                                                                                                                        >  1:1279783/1‑151 (MQ=255)
                                 ccGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGa                                                                                                                       <  8:438936/151‑1 (MQ=255)
                                 ccGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGa                                                                                                                       <  7:309499/151‑1 (MQ=255)
                                    ggAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGaaaa                                                                                                                    <  2:934227/151‑1 (MQ=255)
                                       aaGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGc                                                                                                                 <  5:1123182/151‑1 (MQ=255)
                                         gggggCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCg                                                                                                               <  7:941200/151‑1 (MQ=255)
                                          ggggCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCCGGTTTCCCCCCGCCCGGGGAGGAGACCCCCGCAAACGcccc                                                                                                              >  8:93940/1‑149 (MQ=255)
                                            ggcgggcgTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATc                                                                                                            <  2:146135/151‑1 (MQ=255)
                                              cgggcgTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGGACGGTTTCCCGCGCGGCGGGTGCTGAGCCGCCGGAACAGGcccggcac                                                                                                          >  4:762912/1‑145 (MQ=255)
                                               gggcgTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACg                                                                                                         <  1:1405535/151‑1 (MQ=255)
                                               gggcgTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACg                                                                                                         <  2:1382439/151‑1 (MQ=255)
                                                 gcgTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGt                                                                                                       <  4:411354/151‑1 (MQ=255)
                                                 gcgTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGt                                                                                                       <  1:82208/151‑1 (MQ=255)
                                                  cgTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGGCCGGGTTGCCCCCCGCCGGGGGGTGCGCCGCCCGCCAAGGCCgcgccccggg                                                                                                      >  1:220723/1‑142 (MQ=255)
                                                  cgTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTg                                                                                                      <  7:668500/151‑1 (MQ=255)
                                                   gTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCGGGGTACCGCGCCCCGGGGGGGGAGCCCCCGGcaaacggcgcgcccggggg                                                                                                     >  8:496502/1‑136 (MQ=255)
                                                   gTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCCGGTTGCCGCGCGCCGGGGGCTGAGCCGCCCGCAAAGGCCGCTCCCGgggg                                                                                                     >  5:1213632/1‑151 (MQ=255)
                                                   gTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGGCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTgg                                                                                                     <  1:1096967/151‑1 (MQ=255)
                                                   gTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTgg                                                                                                     <  5:791052/151‑1 (MQ=255)
                                                    tCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCCGTTTACCGCGCGCCCGGTGATGAGACCCCGGAAAAGGCCCCTCAAGGGGGc                                                                                                    >  6:884180/1‑151 (MQ=255)
                                                    tCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGGGATGAGACGCCGGAAAAGGCCGCTCCAGGTGGc                                                                                                    >  4:554788/1‑151 (MQ=255)
                                                      cAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCCGTTTACCGCGCGCCGGGGGATGAGACGCCGGAAAAGGCCGATCACGGTGGGAt                                                                                                  >  6:550937/1‑151 (MQ=255)
                                                        ggAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCa                                                                                                <  5:1255718/151‑1 (MQ=255)
                                                          aaaaTCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCGCGTTGACCCGCGCCCGGGGGTGTGGCGCCGGCAAACGCCGCACACGGGGGcctgacg                                                                                              >  8:752851/1‑145 (MQ=255)
                                                          aaaaTCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGt                                                                                              <  1:291567/151‑1 (MQ=255)
                                                           aaaTCTCGTGCTGACGGTGGCCGGCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTg                                                                                             <  4:1339625/151‑1 (MQ=255)
                                                               gtcGTGCTGACGGTGGCCGTCGCCGCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAACGCCGATCACGGTGGCATCAGTGGTa                                                                                         <  3:486756/149‑1 (MQ=255)
                                                                 cGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTaaaa                                                                                      <  8:418769/151‑1 (MQ=255)
                                                                    gCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGtt                                                                                   <  1:1178527/151‑1 (MQ=255)
                                                                      tGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGGCCCCCGGTTACCGCGCGCCCGGGGGTGAGACCCCGGGAAAAGCCGCGCACCGGGGGCTCAGGGGGAAAAAtttt                                                                                 >  4:679028/1‑151 (MQ=255)
                                                                       gACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGtttt                                                                                 <  5:674569/150‑1 (MQ=255)
                                                                        aCGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGGGCACCGGTTACCCCGCCCCCGGGGAATACACGCCGGCAAAAGCAGAAAACCGGGGCCTCAGGGGGAAAAAGTTTAt                                                                               >  8:221408/1‑151 (MQ=255)
                                                                          ggTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTg                                                                             <  7:1385100/151‑1 (MQ=255)
                                                                          ggTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTg                                                                             <  2:841331/151‑1 (MQ=255)
                                                                           gTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGa                                                                            <  8:1419486/151‑1 (MQ=255)
                                                                           gTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGa                                                                            <  7:357542/151‑1 (MQ=255)
                                                                           gTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACAGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGa                                                                            <  8:1237078/151‑1 (MQ=255)
                                                                             ggCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCCGGTTTCCCCGCGGCCGGGGGTGAGACCCCGGCCAAGGCCGCGCACGGGGGGCTGCGTGGTTAAAAGTttttttatc                                                                           >  8:1110572/1‑148 (MQ=255)
                                                                             ggCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATc                                                                          <  1:933544/151‑1 (MQ=255)
                                                                               ccGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGGGCCCGGGTTGCCGCGCGCCGGGGGGTGGGGCGCCCGCAAACGCCGCGCACGGGGGGGTCAGGGGTTAAAGTtttttttttgcg                                                                        >  1:397870/1‑145 (MQ=255)
                                                                               ccGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCg                                                                        <  1:1060821/151‑1 (MQ=255)
                                                                                cGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCGGTTTCCCGCGCGCCGGGGGATGAGCCGCCGGCAAAGGCCGCTCACGGGGGGCTCAGTGGTAAAAGTTTTTTTGATCCgg                                                                       >  3:802981/1‑151 (MQ=255)
                                                                                 gTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCCGGTTACCGCGCGCCGGGGGGTGAGACGCCGGGAAACGCCCCTCACGGGGGCCTCAGTGGGAAAACTTTTATTGATGCggg                                                                      >  7:1217889/1‑150 (MQ=255)
                                                                                   cgcccagcccGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGTAAGATTCTCATGTGCGAACTGTTTCCCGCGCGCCGGGTGATGAGACGCCTGAAAAGGCCGATCACGGTTGCATCAGTTGTAAAAGTTTTATTGATCCGGaa                                                                     <  6:353130/150‑1 (MQ=255)
                                                                                    gccccgcccGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGGTGAGACGCCGGCAAAGGCCGCTCACCGGGGCATCAGTGGTAAAAGGTTTATTGATCCGGTAgt                                                                   >  3:838093/1‑151 (MQ=255)
                                                                                    gccccgcccGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAgt                                                                   <  3:326227/151‑1 (MQ=255)
                                                                                    gccccgcccGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTACCCCGCGCCGGGGGATGAGACGCCGGAAAAGGCCGCTCACGGTGGCATCACGGGTTAAACTTTTATTGATCCGGaact                                                                   >  1:778983/1‑149 (MQ=255)
                                                                                       ccgcccgcGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGtt                                                                <  5:360659/151‑1 (MQ=255)
                                                                                         gcccgcGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGc                                                              <  3:167608/151‑1 (MQ=255)
                                                                                         gcccgcGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCCCCCCGTTTACCGCCCCCCCCGGGATAAAACCCCCGCAAAAGACCATACAGGTGGCATCAGTTGTTAAAATTTTTTTGATGGGGAAATGtttg                                                              >  2:663462/1‑149 (MQ=255)
                                                                                          cccgcGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGcg                                                             <  7:1148456/151‑1 (MQ=255)
                                                                                            cgcgATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCGGGTTGCCCCGCGCCGGGGGATGGGACGCCGGCAAAGGCCGCTCCCGGTGGGCTGAGGGGGAAAAGTTTTTTTGATCCGGAAGTGtttgggg                                                           >  3:1284624/1‑146 (MQ=255)
                                                                                              cgATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCGGTTTTCCGCGCGCCGGGGGATGATCCGCCGGGGAAGGCCGCTCACGGGGGGCTGAGTGGGTAAAGTTTTTTTTGTTCGGTAGTGTTTggtcgg                                                         >  5:142497/1‑150 (MQ=255)
                                                                                                aTGACAAACAGAGTCAGCCTCGTGACCGCAAGTGGGAAGATTCGCCGTGGCTCTGGGTTTCGGGCGCGTGTGGTGATGAGACGCCGGATTAGGCCGATCTCGGTGGGTTCAGTGGTAAAATTTTTATTGATCCGGAAGTGTTGCGTCGTCa                                                       <  4:1074324/151‑1 (MQ=255)
                                                                                                aTGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGGGCACCGTTTGACCCGCGCCCGGGGGTTAGACCCCCGCAAAAGCCGCGCACCGGGGCATAAGAGGTTAAAATTTTATTGTTACGGGAGGGTTggggtgta                                                        >  6:255399/1‑143 (MQ=255)
                                                                                                  gACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGCTTTATTGATCCGGAAGTGTTGCGTCGTCAGc                                                     <  3:1026451/151‑1 (MQ=255)
                                                                                                   aCAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGCTTTACCGCGCGCCCGGTGGTGAGACGCCGGCAAAAGGCGATCACCGTGGCATCACTGGTAAAAATTTTTTTTATACGGAACTGTTTGGTCGTCAGCg                                                    >  2:1151249/1‑151 (MQ=255)
                                                                                                     aaaCAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTg                                                  <  6:693595/151‑1 (MQ=255)
                                                                                                        cAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCgg                                               <  4:152351/151‑1 (MQ=255)
                                                                                                            gTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGgaaga                                           <  5:255807/151‑1 (MQ=255)
                                                                                                            gTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGgaaga                                           <  5:619013/151‑1 (MQ=255)
                                                                                                             tCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCCGGGGATGAGACCCCGGAAAAGGCCGCTCACGGTGGCCTCAGTGGGTAAAGTTTTATTTAAGCGGAAGTGTTGCGGCGGCAGCGTGCGgcagc                                           >  6:125945/1‑149 (MQ=255)
                                                                                                              cAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACCGTTTGCCGCGCGCCGGGGGATGAGACGCCGGAAAAGGCCGCTCACGGGGGGATCAGTGGTTAAATTTTTTTTTGTGCGGGAGTGTTTGGTGGTGAGCGTGGGggagaa                                          >  4:660860/1‑150 (MQ=255)
                                                                                                                 ccccGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAAccc                                      <  3:1167574/151‑1 (MQ=255)
                                                                                                                  cccGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCGGGTTGCCCCGCGCCGGGTGGTGCGACGCCGGCAAAAGGCGCGGACCGGGGGGTGACTGGTAAAAGTTTTATTGATCCGGGAGTGTTGGGGCGGCAGCGTGCGgattcaaaccc                                     >  1:855105/1‑142 (MQ=255)
                                                                                                                   ccGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCCGGTTTGCCCGGGCCCGGGGCGGGGACGCCGGATAAGGCCGCTGACGGTGGCCGCAGGTGTGAAAATTTTTTTTTTCCGGGAGTTTTTAGGCGTGCCCGTGCGGCAGAAAccccc                                    >  3:899190/1‑149 (MQ=255)
                                                                                                                    cGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCGGGTTGCCGCGCGCCCGGTGGTGAGCCGCCGGAAAAAGCCGCTCACGGGGGGCTCAGTGGGAAAAGTTTTTTTTATCCGGGAGTGGTTGGTGGTGAGCGGGGGGAAGCAACCcccg                                   >  1:340390/1‑151 (MQ=255)
                                                                                                                    cGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCcgcg                                   <  5:1027366/151‑1 (MQ=255)
                                                                                                                     gTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCGGTTTCCCGCGCGCCCGGTGCTGAGACGCCCGAAAAGGCCGATGACGGGGGCATGACTGGTAAAAATTTTATTGATACGGAAGTGTTTCGTCGTCAGCGTGGGGGATAAAACcgcgg                                  >  4:1253870/1‑150 (MQ=255)
                                                                                                                       gACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATGACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCt                                <  8:461381/151‑1 (MQ=255)
                                                                                                                       gACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCt                                <  1:1026313/151‑1 (MQ=255)
                                                                                                                          cGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCC‑‑AAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGc                           <  2:351654/151‑1 (MQ=255)
                                                                                                                               gTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGaa                        <  4:732749/151‑1 (MQ=255)
                                                                                                                                  ggAAGATTCGCCGTGGCGCACGGTTTCGCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGGTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAAt                     <  1:545661/151‑1 (MQ=255)
                                                                                                                                     aGATTCGCCGTGGCGCCCCGTTTACCCCGCGCCGGGGGATGAGACCCCCGAAAAGGCCGATCACCGTGGCATCAGTGGTAAAAGTTTTATTTATACCGAAGTGTTGCGGCGTCAGCGTGCGCAGAAACCCCCCGCTACGGCGCAAAAGCAt                 >  4:1057853/1‑151 (MQ=255)
                                                                                                                                     aGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCa                  <  6:1480180/151‑1 (MQ=255)
                                                                                                                                            ccGTGGCGCCCCGTTTACCGCGCCCCGGGTGATGAGACCCCCGAAAACGCCGATCACGGTGGCCTCAGTGGTAAAAGTTTTATTTATCCGGAAATGTTTCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGGATGTGAgg           >  8:795062/1‑151 (MQ=255)
                                                                                                                                             cGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGc          <  7:508873/151‑1 (MQ=255)
                                                                                                                                              gTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTCCGGCGAAAATGCATGTCAGGcc         >  5:163145/1‑150 (MQ=255)
                                                                                                                                                  gcgcACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGCACtt      <  5:986595/151‑1 (MQ=255)
                                                                                                                                                    gcACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGCACtttt    >  2:751770/1‑151 (MQ=255)
                                                                                                                                                    gcACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGCACtttt    <  2:1228160/151‑1 (MQ=255)
                                                                                                                                                      cACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGCACTTTTc   <  1:1458662/151‑1 (MQ=255)
                                                                                                                                                       |                                                                                                                                                       
ATGATTCTGATAAACGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGCACTTTTCC  >  DH1/1165295‑1165594

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