Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH11,165,4440AC26.2% 388.7 / 34.6 169T71P (ACG→CCG) yfiFmethyltransferase
Reads supporting (aligned to +/- strand):  ref base A (16/108);  new base C (44/0);  total (60/109)
Fisher's exact test for biased strand distribution p-value = 2.37e-27
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 173 total aligned reads displayed.

TCTGATAAACGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGCACTTTTCC  >  DH1/1165300‑1165594
                                                                                                                                                 |                                                                                                                                                           
tCTGATAAACGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGtt                                                                                                                                                       <  4:741093/150‑1 (MQ=255)
  tGATAAACGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCGGgttgc                                                                                                                                                    >  1:410691/1‑149 (MQ=255)
   gATAAACGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGttt                                                                                                                                                      <  1:203714/148‑1 (MQ=255)
    aTAAACGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCcg                                                                                                                                                  <  5:350047/151‑1 (MQ=255)
       aaCGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGGGCCGGGGTTGCcgccc                                                                                                                                               >  6:765735/1‑149 (MQ=255)
       aaCGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCGGGTTCCcgcgg                                                                                                                                               >  3:189449/1‑150 (MQ=255)
       aaCGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCGGGTTCCcgcgc                                                                                                                                              >  1:574053/1‑151 (MQ=255)
       aaCGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCcgcgc                                                                                                                                              <  8:648367/151‑1 (MQ=255)
         cGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCCGGTTGCCgcgcgc                                                                                                                                            >  3:724218/1‑151 (MQ=255)
            aCCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGGCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGGCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCgg                                                                                                                                         <  5:804311/151‑1 (MQ=255)
             cccACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAATCAGAGTCAGCCCCGTGACCGCAAGTCGGAAGATTCGCCGTGGGGGAAGGGTTGCCGCGCGCCggg                                                                                                                                        >  6:1065570/1‑151 (MQ=255)
              ccACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACTGTCGCCGTCGCCCCGCACGCGATGACAAACAGAGTCAGCCCCGTGACAGCAAGTGGGAAGATTAGCCGTGGGGCACAGTTTCCAgcgggcgccgt                                                                                                                                       >  6:894536/1‑146 (MQ=255)
                 caACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGCCAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGGCCCGGTTCCCGCGCGGCggggg                                                                                                                                      >  8:899742/1‑147 (MQ=255)
                 caACCCGCGTACAGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCCGGACCGCAAGTGCGAAGCTTCGCCGAGGGCCCCGGGGTCCcgcgcgggcggggcc                                                                                                                                    >  6:811621/1‑142 (MQ=255)
                   aCCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGtgatga                                                                                                                                  <  3:84877/151‑1 (MQ=255)
                      cgcgTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCGGGTTGCCGCCCGCCGGGGGGTGAGAc                                                                                                                               >  1:371122/1‑151 (MQ=255)
                      cgcgTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCCGGTTGCCGCGCGCCggggggtgggcc                                                                                                                               >  5:184808/1‑142 (MQ=255)
                           aCCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGGGGCCCGGGTTCCGCGCGCCGCGGGGTGAGCCCCCgg                                                                                                                         >  7:19068/1‑151 (MQ=255)
                             cGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGGCCGGTTTGCCGCGCGCCGGGGGATGAGCCGCCggc                                                                                                                        >  1:145822/1‑149 (MQ=255)
                             cGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCGGGTTCCCGCGCGCCCGGGGCTGAGACCCCggc                                                                                                                        >  7:992768/1‑149 (MQ=255)
                             cGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGaa                                                                                                                       <  3:307861/151‑1 (MQ=255)
                              gggAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGaaa                                                                                                                      <  3:50154/151‑1 (MQ=255)
                                gAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAg                                                                                                                    <  7:407260/151‑1 (MQ=255)
                                gAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGGCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAg                                                                                                                    <  5:603152/151‑1 (MQ=255)
                                  aaGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTGCCCCGCGCCGGGTGCTGAGCCGCCGGAAAAGGc                                                                                                                  >  2:89416/1‑151 (MQ=255)
                                   aGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGcc                                                                                                                 <  3:1055015/151‑1 (MQ=255)
                                     ggggCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGGCCGCGTTTCCCCCCCCCCGGGGGTGAGCCGCCGCGAAAGGCCgc                                                                                                               >  5:292315/1‑150 (MQ=255)
                                      gggcgggcGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGAt                                                                                                              <  4:287785/151‑1 (MQ=255)
                                         cgggagTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCGGGTTGCCGCGCGCCGGGGGGTGAGCCGCCGGGAAAgggccggccc                                                                                                           >  6:1021137/1‑143 (MQ=255)
                                          gggcgTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACg                                                                                                         <  5:711227/151‑1 (MQ=255)
                                          gggcgTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACg                                                                                                         <  8:734448/151‑1 (MQ=255)
                                              gTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCCGGTTACCGCGCGCCGGGGGGTGCGGCGCCGGCAAAGGccccgcacgggtg                                                                                                      >  1:85761/1‑140 (MQ=255)
                                              gTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTgg                                                                                                     <  5:981326/151‑1 (MQ=255)
                                                ccAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGGGGCCGGTTTACCGCGCGCCGGGGGGTGAGTCGCCGGCAAAGGCCGCGCACGGTgggc                                                                                                   >  5:1007453/1‑149 (MQ=255)
                                                ccAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCGGGTTCCCGCGCGCCGGGGGTTGGGCCGCCGGCAAAGGCCGCTCACCGTGGcc                                                                                                   >  3:578002/1‑150 (MQ=255)
                                                ccAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCa                                                                                                   <  6:118949/151‑1 (MQ=255)
                                                    gAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAg                                                                                               <  1:91533/151‑1 (MQ=255)
                                                       aaTCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAACGGCCGATCACGGTGGCATCAGTgg                                                                                            <  5:46741/151‑1 (MQ=255)
                                                       aaTCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATAAGTgg                                                                                            <  5:1066105/151‑1 (MQ=255)
                                                          ctcGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCGGGTTGCCGCGCGCCGGGGGGTGGGCCGCCGGAAAAGGCCGCGCACCGTGGGCTCAGTggga                                                                                          >  5:443838/1‑148 (MQ=255)
                                                            cGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTa                                                                                          <  2:1053685/148‑1 (MQ=255)
                                                            cGTGCTGACGCTGGCCGGCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGGCGGGTGATGAGACGCCGGAAAAGTCCGATCACGGTGGCATCAGTGGTa                                                                                          <  3:112887/148‑1 (MQ=255)
                                                             gTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCCGGTTGCCCCGCGCCCGGGGGTGTGACGCCCGGAAAAGGCGCTCACGGGGGGCGCAGTgggtaa                                                                                        >  1:399197/1‑145 (MQ=255)
                                                             gTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAg                                                                                      <  3:976714/151‑1 (MQ=255)
                                                                  gACGGTGGCCGTCGCCCCGCCGTCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGttt                                                                                   <  3:338450/149‑1 (MQ=255)
                                                                  gACGGTGGCCGTCGCCCCGCCCGCGATGACAAAGAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCTCGGTTTCCCGCGCGCCGGGAGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGtttt                                                                                  <  6:215520/150‑1 (MQ=255)
                                                                  gACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTa                                                                                 <  4:167079/151‑1 (MQ=255)
                                                                   aCGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTAt                                                                                <  8:80856/151‑1 (MQ=255)
                                                                    gggTGGCCGTCGCCGCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCTGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTAtt                                                                               <  7:679506/150‑1 (MQ=255)
                                                                    cGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTAtt                                                                               <  8:88530/151‑1 (MQ=255)
                                                                      gTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGGGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGa                                                                             <  4:587377/151‑1 (MQ=255)
                                                                       gggCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGAt                                                                            <  1:516044/150‑1 (MQ=255)
                                                                        ggCCGTCGCCCCGCCCGCGATGACAATCAGAGTCAGCGCGGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCGCGCGCGCCGGGTGATGAGACGCCGGAATAGGGCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGAt                                                                            <  2:69122/150‑1 (MQ=255)
                                                                        ggCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGATGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGAt                                                                            <  8:1049176/150‑1 (MQ=255)
                                                                            gTCGCCCCGCCAGCTATTACAAACAGAGTCAGCCACGTGACAGCAAGTGGGAAGATTCGCCGTGGCGGACGTTTTCCCGGGCGCCGTGTGATGAGACGCCGGAAAAGGCCGATAACGGTTGCATCTGTGGTAAAAGTTTTATTGATCCGGa                                                                       <  3:1033775/151‑1 (MQ=255)
                                                                             gcgccccgcccGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGaa                                                                      <  6:224631/150‑1 (MQ=255)
                                                                              cgccccgcccGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGaa                                                                      <  1:491499/150‑1 (MQ=255)
                                                                              cgccccgcccGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAg                                                                     <  6:589020/151‑1 (MQ=255)
                                                                              cgccccgcccGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAg                                                                     <  3:174549/151‑1 (MQ=255)
                                                                              cgccccgcccGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAg                                                                     <  8:302670/151‑1 (MQ=255)
                                                                                ccccgcccGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCGGGTTCCCGCGCGCCGGGGGATGAGCCGCCGGAAAACGGCGCGCCCGGTGGGCTCAGTGGGTAAAGTTTTTTTGATACGGCAgtg                                                                   >  5:1082290/1‑151 (MQ=255)
                                                                                ccccgcccGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAgtg                                                                   <  6:955331/151‑1 (MQ=255)
                                                                                  ccgcccgcGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCGGGTTACCGCCCGCCCGGGGATGAGACGCCGGAAAAGGGCCAACAACGTGGCATACGTGGTAAAAATTTTATTGGTCCGGCAGTGtt                                                                 >  6:270823/1‑151 (MQ=255)
                                                                                  ccgcccgcGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGtt                                                                 <  4:190554/151‑1 (MQ=255)
                                                                                    gcccgcGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGc                                                               <  7:859452/151‑1 (MQ=255)
                                                                                    gcccgcGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGc                                                               <  4:812729/151‑1 (MQ=255)
                                                                                    gcccgcGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGc                                                               <  6:548626/151‑1 (MQ=255)
                                                                                     cccgcGATGACAAACAGAGTCAGCCCCGTGAGCGGATGTGGGAAGATTCGCCGTGGCGCACGGGTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTTTTGATCCGGAAGTGTTGcg                                                              <  4:37109/151‑1 (MQ=255)
                                                                                         ggATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGcgtcgt                                                          <  8:935213/150‑1 (MQ=255)
                                                                                           aTGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCAGATCACTGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCa                                                        <  7:262603/151‑1 (MQ=255)
                                                                                            tGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACCGTTTGCCCCGCGCCGGGGGATGAGGAGCCGGCAAAGGCCGCTGACGGGGGCCTGAGGGGGTAAAGTTTTTTTGATCCCGCAGTGTTgcggcggc                                                         >  3:444400/1‑147 (MQ=255)
                                                                                             gACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAg                                                       <  1:16655/150‑1 (MQ=255)
                                                                                                     gagTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGgaa                                             <  5:558715/151‑1 (MQ=255)
                                                                                                       gTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTGCCGCCCGCCCGGGGATGAGACCCCGGAAAAGGCCGCTCACCGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGGAGTGTTTCGTCGTCAGCGTGCggg                                              >  4:553937/1‑147 (MQ=255)
                                                                                                          aGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAAc                                        <  2:339200/151‑1 (MQ=255)
                                                                                                             cccGTGACCGCAAGTGGGAAGATTCGCCGTGGCACACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCcg                                     >  5:532722/1‑151 (MQ=255)
                                                                                                               cGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCcgcg                                   <  4:866532/151‑1 (MQ=255)
                                                                                                               cGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCcgcg                                   <  3:100639/151‑1 (MQ=255)
                                                                                                                gTGACCGCAAGTGGGAAGATTCGCCGTGGCGCCCGGGTTGCCCCGCGCCGGGGGGTGAGACCCCGGCAAAGGGCGCTCACCGTGGGCTCAGTGGTTAAAATTTTTTTTATACGGGAGTGTTGGGGGGGCAGCGGGCGGGAGACACCcgcgg                                  >  7:538594/1‑150 (MQ=255)
                                                                                                                   aCCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTa                               <  2:1040364/151‑1 (MQ=255)
                                                                                                                   aCCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTa                               <  2:1035448/151‑1 (MQ=255)
                                                                                                                   aCCGCAAGTGGGAAGATTCGCCGTGGCCCCCGGTTTCCCCCGCCGCCGGGGAGTAGGCGGCGGGAAAGGGCGATACACGTGGCATGAGTGGGTAAAGGTTTTTTTATACCGACGGGTTGGGTCGTGCGCGTGCGGGACGAAAccccggctc                                >  3:972316/1‑145 (MQ=255)
                                                                                                                    ccGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTAc                              <  2:28570/151‑1 (MQ=255)
                                                                                                                       cAAGTGGGAAGATTCGCCGTGGCGCCCGGTTTCCCGCGCGCCGGGTGATGAGCCGCCGGGAAAGGCCGCTCACGGTGGCCTGAGTGGGAAAACTTTTATTGATCCGGAAGTGTTGGGTCGTCACCGTGCGGGAGAAACCCGCGGCTAAGGc                           >  3:607174/1‑151 (MQ=255)
                                                                                                                        aaGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCg                          <  4:187186/151‑1 (MQ=255)
                                                                                                                          gTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGaa                        <  1:724272/151‑1 (MQ=255)
                                                                                                                            gggAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGaaaa                      <  8:509102/151‑1 (MQ=255)
                                                                                                                            gggAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGaaa                       <  2:675559/150‑1 (MQ=255)
                                                                                                                               aaGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGGGATGAGACGCCGGAAAAGGCCGATCAAGGGGGCATCAGTGGTATAACTTTTATTGATACGGAAGTGTTTGGTCGTGAGCGTGGGGGAGAATCCCGCGTCTACGGCGCGAATGc                   >  2:558836/1‑151 (MQ=255)
                                                                                                                                 gATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCAt                 <  7:59472/151‑1 (MQ=255)
                                                                                                                                 gATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCAt                 <  6:1033167/151‑1 (MQ=255)
                                                                                                                                  aTTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCATg                <  5:270003/151‑1 (MQ=255)
                                                                                                                                     cGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCa             <  7:353040/151‑1 (MQ=255)
                                                                                                                                       ccGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAgg           >  4:529668/1‑151 (MQ=255)
                                                                                                                                         gTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGCa         <  7:451499/151‑1 (MQ=255)
                                                                                                                                             gcgcACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGTCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGCACtt      <  7:553228/151‑1 (MQ=255)
                                                                                                                                                cACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGCACTTTTc   >  8:459018/1‑151 (MQ=255)
                                                                                                                                                cACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGCACTTTTc   >  2:505740/1‑151 (MQ=255)
                                                                                                                                                 aCGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGCACTTTTcc  <  1:118377/151‑1 (MQ=255)
                                                                                                                                                 aCGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGCACTTTTcc  >  5:367448/1‑151 (MQ=255)
                                                                                                                                                 |                                                                                                                                                           
TCTGATAAACGTACCCACAACCCGCGTACCGGGAAAGGGGGCGGGCGTCCAGGAAAATCTCGTGCTGACGGTGGCCGTCGCCCCGCCCGCGATGACAAACAGAGTCAGCCCCGTGACCGCAAGTGGGAAGATTCGCCGTGGCGCACGGTTTCCCGCGCGCCGGGTGATGAGACGCCGGAAAAGGCCGATCACGGTGGCATCAGTGGTAAAAGTTTTATTGATCCGGAAGTGTTGCGTCGTCAGCGTGCGGAAGAAACCCGCGTCTACGGCGAAAATGCATGTCAGGCACTTTTCC  >  DH1/1165300‑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: