Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH11,165,4440AC23.2% 1046.1 / 77.5 444T71P (ACG→CCG) yfiFmethyltransferase
Reads supporting (aligned to +/- strand):  ref base A (58/280);  new base C (102/0);  total (164/280)
Fisher's exact test for biased strand distribution p-value = 0.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 449 total aligned reads displayed.

ATGATTCTGATAAACGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGCACTTTTCC  >  DH1/1165295‑1165594
                                                                                                                                                     |                                                                                                                                                        
      cTGATAAACGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCGGGTTc                                                                                                                                                  >  8:685094/1‑151 (MQ=255)
         aTAAACGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCcg                                                                                                                                               <  4:411766/151‑1 (MQ=255)
           aaaCGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGGCCGGGTTGCcgcg                                                                                                                                             >  2:408642/1‑151 (MQ=255)
           aaaCGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCGGTTTcccccc                                                                                                                                             >  4:865217/1‑148 (MQ=255)
           aaaCGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCcgcg                                                                                                                                             <  7:1697703/151‑1 (MQ=255)
             aCGTACCCAGAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGAAGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCcgcgcg                                                                                                                                           <  2:129122/151‑1 (MQ=255)
              cGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGGCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCgcgcgc                                                                                                                                          <  4:2102021/151‑1 (MQ=255)
                   ccACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGt                                                                                                                                     <  4:1901641/151‑1 (MQ=255)
                   ccACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGt                                                                                                                                     <  7:718877/151‑1 (MQ=255)
                    cacaACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCGGGTTCCCGCGCGCCcgggg                                                                                                                                    >  3:1441172/1‑149 (MQ=255)
                    cacaACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGtg                                                                                                                                    <  5:15840/151‑1 (MQ=255)
                     acaACCGGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGtga                                                                                                                                   <  5:2140915/151‑1 (MQ=255)
                     acaACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTGCCCCCCGCCgggggg                                                                                                                                   >  1:101294/1‑148 (MQ=255)
                     acaACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGtga                                                                                                                                   <  1:1602959/151‑1 (MQ=255)
                         cccGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGGGGCCGGGGTCCCGCGGGCCgcgggggggg                                                                                                                               >  6:1967685/1‑142 (MQ=255)
                           cgcgTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCCCCCCGTTTCCCCCCCGCCGGGggatgcgcc                                                                                                                             >  3:1319413/1‑147 (MQ=255)
                             cgTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCGGGTTGCCGCGCGCCGGGGGATGAGACGc                                                                                                                           >  5:349902/1‑151 (MQ=255)
                              gTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCGGGTTACCCCGCGCCcgggggggccccgcc                                                                                                                          >  3:587153/1‑139 (MQ=255)
                                aCCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGGGGAGCGGTTGCCGGCGCCCgggggggggggcgccgg                                                                                                                         >  8:1570565/1‑138 (MQ=255)
                                aCCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCgg                                                                                                                        <  6:1158811/151‑1 (MQ=255)
                                  cGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGaa                                                                                                                     <  4:1769619/151‑1 (MQ=255)
                                     gAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAg                                                                                                                  <  6:1764499/151‑1 (MQ=255)
                                     gAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTAGCCGTCGCCCCGCCCGCGATGACAAACAGACTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAg                                                                                                                  <  2:152731/151‑1 (MQ=255)
                                      aaaGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAgg                                                                                                                 <  1:1175321/151‑1 (MQ=255)
                                         gggggCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCg                                                                                                              <  1:1738397/151‑1 (MQ=255)
                                         gggggCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCg                                                                                                              <  1:2263045/151‑1 (MQ=255)
                                             gcgggcgTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCa                                                                                                          <  4:2247455/151‑1 (MQ=255)
                                                ggcgTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGGGCCCCGGTTGCCGCCCGCCCGGGGGTGAGACCCCGGCAAAAGGcccagacgg                                                                                                        >  6:476533/1‑144 (MQ=255)
                                                     ccAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCa                                                                                                  <  8:2072883/151‑1 (MQ=255)
                                                     ccAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACCGTTTACCCCGCGCCCGGTGATGAGACGCCGGCAAAGGCCCCTCACGgggggc                                                                                                  >  8:1720104/1‑149 (MQ=255)
                                                      cAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGGCCGGGTTGCCGCGCGCCCGGGGGTGCGACGCCCGAAAACGCCgcgccccgggggcg                                                                                                 >  4:773976/1‑138 (MQ=255)
                                                       aGGAAAATCTCGTTCTGACGGTGGCCGTCGCCACGCCCGCGATGACAAACAGAGTAAGCCCCGTGACCGCAAGTGGGAAGATTCGCCCTCGCGCCCGGTTTACAGGGCGAACGGGGATGAGACGGCGGAAAAAGACGCTAACGGTGGCaa                                                                                                 >  6:371368/1‑149 (MQ=255)
                                                       aGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATc                                                                                                <  3:532294/151‑1 (MQ=255)
                                                          aaaaTCTCGTGCTTACGGTGGCCGTCGCCCAGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGATACGCCGGAAAAGGCCGATCACGGTGGCATCAg                                                                                              <  6:983644/150‑1 (MQ=255)
                                                           aaaTCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTg                                                                                            <  6:792406/151‑1 (MQ=255)
                                                           aaaTCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTg                                                                                            <  4:644987/151‑1 (MQ=255)
                                                               ctcGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCGGGTTGCCGCGCGCCGGGGGGTGAGACGCCCGCAAAGGCCGCTCCCGGTGGGCTCAgggggt                                                                                         >  1:1390195/1‑148 (MQ=255)
                                                                 cGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGt                                                                                          <  4:1666163/147‑1 (MQ=255)
                                                                 cGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTa                                                                                         <  4:1003623/148‑1 (MQ=255)
                                                                       gACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCCCACCGGTTGCCGCGCGCCGGGGGGTGTGACGCCGGCAAAGGCCGCTCACCGGGGGCTCACTGGTAAAAGgtttt                                                                                >  4:1979466/1‑150 (MQ=255)
                                                                        aCGGTGGCCTTAGCCCCGCCCGAGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACTGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTTGTAAAAGTTTTAt                                                                               <  6:1473685/151‑1 (MQ=255)
                                                                        aCGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTAt                                                                               <  4:554705/151‑1 (MQ=255)
                                                                         cGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGGGGCCGGGTTGGCGCGCGGCGGGGGCTGGGACGCCGGGCAAGGGCGCTCGCGGGGGGGTGAGGGGTTAAAGTtttttt                                                                              >  8:1310579/1‑151 (MQ=255)
                                                                         cGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTAtt                                                                              <  5:1816551/151‑1 (MQ=255)
                                                                          ggTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCCGGTTACCGCGCGCCCGGGGCTGAGACCCCCGCAAACGCCGCTCACCGGGGCAGGAGGGGTGAAAATttttttt                                                                             >  6:1253484/1‑150 (MQ=255)
                                                                          ggTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTg                                                                             <  8:22810/151‑1 (MQ=255)
                                                                           gTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGa                                                                            <  3:774063/151‑1 (MQ=255)
                                                                            gggCCGGCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGAt                                                                           <  8:287611/150‑1 (MQ=255)
                                                                               ccGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCg                                                                        <  2:920372/151‑1 (MQ=255)
                                                                               ccGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCg                                                                        <  3:396612/151‑1 (MQ=255)
                                                                                   cgccccgcccGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGaa                                                                     <  2:1629349/150‑1 (MQ=255)
                                                                                   cgccccgcccGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTACCGCGCGCCGGGGGGTGAGCCGCCGGAAAAGGCCGCTCACGGTGGCCTCAGTGGTAAAATTTTTTTTGCTCCggc                                                                      >  3:178170/1‑148 (MQ=255)
                                                                                    gccccgcccGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGGGCCCCGGTTGCCGCGCGCCGGGTGGTGAGCCGCCCGCAAAGGCCGCGCACGGGGGGATCAGTGGTTAAACTTTTTTTTATCCGGaatt                                                                   >  5:1128128/1‑149 (MQ=255)
                                                                                    gccccgcccGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCCGTTTCCCCCGCGCCGGGGGGTGAGCCGCCGGGAAAGGCCGCTCCCGGGGGCCTCAGTGGTAAAAGTTTTATTTATCCCGAAgt                                                                   >  1:514053/1‑151 (MQ=255)
                                                                                    gccccgcccGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAgt                                                                   <  4:1317331/151‑1 (MQ=255)
                                                                                      cccgcccgCGATGACAAACAGAGTCAGCCGCGTGACCGCAAGTGGGAAGATTCGCCGTCGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGGTTTATTGATCCGGAAgtgt                                                                 <  7:1306505/151‑1 (MQ=255)
                                                                                      cccgcccgCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAgtgt                                                                 <  3:1590804/151‑1 (MQ=255)
                                                                                       ccgcccgcGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCGGGTTCCCGCGCGCCGGGGGGTGAGGCGCCGGGAAAGGCCGCGCCCCGGGGCATCAGTGGTCAAAGTTTTATTGAGACGGTCGTGtt                                                                >  5:756861/1‑151 (MQ=255)
                                                                                        cgcccgcGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTg                                                               <  2:222924/151‑1 (MQ=255)
                                                                                          cccgcGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGcg                                                             <  8:143139/151‑1 (MQ=255)
                                                                                          cccgcGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGcg                                                             <  3:2231546/151‑1 (MQ=255)
                                                                                           ccgcGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGGCCCGTTTGCCCCGCGCCGGGGGGTGCGCCGCCCGCAAAGGCCGATCACCGTGGCATCAGTGGTAAAACTTTTATTTATCCGGCAGTGTTGcgg                                                            >  1:2008573/1‑150 (MQ=255)
                                                                                           ccgcGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGcgt                                                            <  1:1270898/151‑1 (MQ=255)
                                                                                            cgcgATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGcgtc                                                           <  7:1920926/151‑1 (MQ=255)
                                                                                             gcgATGACAAACAGAGTCAGCCCCGGGACCGCAAGTGGGAAGATTCGCCGGGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGGGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGcgtcg                                                          <  6:231365/151‑1 (MQ=255)
                                                                                                aTGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCa                                                       <  3:67749/151‑1 (MQ=255)
                                                                                                   aCAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACCGTTTACCGCCCGCCGGGTGAAGAGACCCCGGAAAAAGCCGCAAACCGTGGGATCACTGGGTAAAATTTTAATTATGCCGAAGTGTTTCGGGGTCAGcc                                                    >  6:1017979/1‑150 (MQ=255)
                                                                                                    cAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGt                                                   <  5:2251957/151‑1 (MQ=255)
                                                                                                     aaaCAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGGCCGGTTTGCCGCCCGCCGGGTGAGGAGCCCCCGGCAAACGCGGCTCACGGTGGGAACAGTGGTTAAAATTTTTTTGATGCGGAAGAGTTTCGTCGTCAGCGTg                                                  >  4:330126/1‑151 (MQ=255)
                                                                                                     aaaCAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTg                                                  <  3:2049599/151‑1 (MQ=255)
                                                                                                         agagTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGGCAGCGTGCGga                                              <  2:2005906/151‑1 (MQ=255)
                                                                                                             tCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGGGCCCGGTTTGCCCCGCGCCGGGGGGTGAGGCGCCCGCAAAGGCCGCTCACGGTGGGCTGAGTGGTTAAAATTTTTTTTATCCGGCAGTGTTTGGTCGTGAGCGTGGGgcaga                                           >  1:968590/1‑150 (MQ=255)
                                                                                                               aGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAAc                                        <  1:4206/151‑1 (MQ=255)
                                                                                                                 ccccGTGACCGCAAGTTGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAAccc                                      <  6:1559354/151‑1 (MQ=255)
                                                                                                                     gTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCGGTTTTCCGCGCGCCGGGTGGTGAGGCGGCGGATAAGGCGGCTGACGGTGGGGTGAGTGGTTAAATTTTTTTTTGTTCGGAAGTGTTGGGTCGTGAGCGTGGGGTAGCAAACcgcgg                                  >  5:217656/1‑150 (MQ=255)
                                                                                                                       gACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTGCCGCGCGGCGGGTGGTGAGACGCCGGGAAAGGCCGCTCCCGGGGGCATCAGTGGGTAAAGTTTTTTTGATACGGCAGTGTTGCGTCGGGAGCGTGGGGCAGCAACCCGCGTCt                                >  3:873364/1‑151 (MQ=255)
                                                                                                                          cGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTACCGCGCGCCCGGTGATTAGACGCCGGCAAAGGCCGCTGACGGGGGGATGCGTGGTTAAAGTTTTTTTTATGCGGCAGTGTTTCGTCGTGAGCGTGGGGGACCAACCcgcggcgtcc                             >  2:816308/1‑145 (MQ=255)
                                                                                                                            cAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGc                           <  5:1720152/151‑1 (MQ=255)
                                                                                                                            cAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGc                           <  5:519663/151‑1 (MQ=255)
                                                                                                                            cAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGc                           <  7:582629/151‑1 (MQ=255)
                                                                                                                             aaGTGGGAAGATTCGCCGTGGCGCCCGGTTTGCCGCGCGCCGGGTGGTGAGCCGCCGGCAAAGGCCGGTGACGGTGGGCTGAGTGGGTAAAGTTTTATTGTTGCGGAAGTGTTGGGGGGTGAGCGGGGGGCAGAAACCcgcgggttcgggg                          >  3:1243693/1‑142 (MQ=255)
                                                                                                                               gTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGaa                        <  1:749338/151‑1 (MQ=255)
                                                                                                                                 gggAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGaaaa                      <  7:1268544/151‑1 (MQ=255)
                                                                                                                                 gggAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGaaaa                      <  5:1429524/151‑1 (MQ=255)
                                                                                                                                 gggAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGaaa                       <  8:553924/150‑1 (MQ=255)
                                                                                                                                  ggAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAACGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAAt                     <  7:378550/151‑1 (MQ=255)
                                                                                                                                   gAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATg                    <  5:97342/151‑1 (MQ=255)
                                                                                                                                    aaGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGc                   <  5:448567/151‑1 (MQ=255)
                                                                                                                                    aaGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGc                   <  6:2112626/151‑1 (MQ=255)
                                                                                                                                      gATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCAt                 <  4:1474013/151‑1 (MQ=255)
                                                                                                                                      gATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCAt                 <  2:1696245/151‑1 (MQ=255)
                                                                                                                                           gCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATTAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTTATCCGGAAGTTTTGCGTCGTCAGCGTGGGGAAGAAACCCGCGTCTACGGCGAAAAATCATGta              >  4:1152043/1‑148 (MQ=255)
                                                                                                                                              gTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGcc         >  2:712477/1‑150 (MQ=255)
                                                                                                                                              gTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGCa         <  5:1065952/151‑1 (MQ=255)
                                                                                                                                                ggCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGCAc        <  4:587490/150‑1 (MQ=255)
                                                                                                                                                ggCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGCACt       >  2:187283/1‑151 (MQ=255)
                                                                                                                                                  cgcACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGCACttt     <  8:2019494/151‑1 (MQ=255)
                                                                                                                                                  cgcACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGCACttt     <  2:1159593/151‑1 (MQ=255)
                                                                                                                                                   gcACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGCACtttt    <  8:1029002/151‑1 (MQ=255)
                                                                                                                                                   gcACGGGTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGCACtttt    <  5:67242/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: