Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH12,326,2870GT31.5% 668.9 / 126.0 407intergenic (+207/+200)yddM/fdnItranscriptional regulator/formate dehydrogenase nitrate‑inducible cytochrome b556(Fdn) subunit
Reads supporting (aligned to +/- strand):  ref base G (178/94);  new base T (37/88);  total (216/191)
Fisher's exact test for biased strand distribution p-value = 3.88e-11
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 419 total aligned reads displayed.

GATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATG  >  DH1/2326137‑2326437
                                                                                                                                                      |                                                                                                                                                      
gATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTg                                                                                                                                                        <  5:1231911/151‑1 (MQ=32)
 aTGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGATGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTGa                                                                                                                                                       <  8:575172/151‑1 (MQ=18)
   gcgcTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTGaaa                                                                                                                                                     <  5:1311963/151‑1 (MQ=34)
   gcgcTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTGaaa                                                                                                                                                     <  6:1650244/151‑1 (MQ=34)
         cGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTGAAACttat                                                                                                                                                >  3:1129950/1‑150 (MQ=25)
         cGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTGAAACtgatg                                                                                                                                               <  3:1402176/151‑1 (MQ=39)
          gCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGACATATGaagatcatc                                                                                                                                               >  8:283030/1‑143 (MQ=255)
            ttATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGCCGCCTTTTACAGCGCCATTTTAAACTGATGAc                                                                                                                                             >  7:1343455/1‑150 (MQ=25)
              aTCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTGAAACTGATGACaaa                                                                                                                                          <  1:551652/151‑1 (MQ=255)
                 aGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTTAAACTTATGAaaaa                                                                                                                                          >  8:1712855/1‑148 (MQ=255)
                  ggCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTGAAACTGATGACAAACGCa                                                                                                                                      <  3:193486/151‑1 (MQ=255)
                    ccTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGCCGCATTTTACAGCGCCATTTGAAACTGATGACAAAACCAAc                                                                                                                                    >  3:1745995/1‑151 (MQ=255)
                     atatattttCACTGTATTCTATTTGATTTGTATTGTTTTATTTCTGTATATGTCTTTAAAGCCTCATCCTTTATGAACAAAGCCTACTTTTTATACAATATTAAAATGACGCCTTTTACAGCGCCATTTTAAACTGATGACAAACGCAACa                                                                                                                                   <  2:938030/144‑1 (MQ=255)
                      tACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTTAAACTGATTAAAAACGCAACAt                                                                                                                                  >  1:949504/1‑151 (MQ=255)
                      tACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAACAACATAAGAAACTCATGAaaaaaaaaacat                                                                                                                                  >  8:457504/1‑143 (MQ=11)
                       acacGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGACGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTGAAACTGATGACAAACGCAACAtt                                                                                                                                 <  8:644251/151‑1 (MQ=255)
                         acGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTTAAACTTATAAAAAAAGCAACATTg                                                                                                                                >  6:530066/1‑150 (MQ=255)
                           cgtCACTGCAATCTATTGGATTTGAAGGGTTTTTTTTTTTTATTATGCGTTCACGCCGCATATGTAATGAAAAAATACTACTTTTTCTATAATATTATATTTTCGGATTTTAAAGCGCCATTTTTAACTGATGACAAACGCAACATTGCCt                                                                                                                             <  8:1821256/149‑1 (MQ=9)
                              cACTGCAATCTATTGTATTTGAAGGGTTTTCATGCCGTATAAGGCTTTAACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATTAAAATGACTCCTTTTACAGCGCCATTTTAAACTGATGACAAACGCAACATTGCCTGAt                                                                                                                          <  4:1708065/151‑1 (MQ=17)
                              cACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTTAAACTGATGACAAAAGCAACATTGCCTGAt                                                                                                                          >  5:929305/1‑151 (MQ=255)
                               actgcctggtATTGGCTTTGAAGGGTTGTCATTGCGTCTAAGCGTTTCACGCCGCATGCGGCATGAACAAAGCCTTCGTTGTCTACAATATGAAAATGACTCCTTTTACAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATg                                                                                                                         <  8:428956/142‑1 (MQ=25)
                               aCTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTTAAACTGATGAAAAACGAAACATTGACTGATg                                                                                                                         >  5:234337/1‑151 (MQ=255)
                                 tGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTTAAAGTGATGACAAACGCAACATTGCCTGATgcg                                                                                                                       >  8:251978/1‑151 (MQ=255)
                                  gttatgtATTGTATTTTAAGGGTTTTCATTCCGTTTAAGGGTTTCACGCCGCATCCGGCATGAACAAAGCCTACTTTGTCTACAATATTAAAATGTCTCCTTTTACAGCGCCATTTTAAACTGATGACAAACGCAACATTGCCTGATgcg                                                                                                                       <  2:1374496/144‑1 (MQ=17)
                                  ggattCTATTGGATTTGAAGGGTTTTCATGGGGGTTATGGCTTTCACGCCGCATGGGGGATGTACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTTAAACTGATGACAAACGCAACATTGCCTGATgcg                                                                                                                       <  4:934432/146‑1 (MQ=17)
                                  gCAATCTATTGTATTTGAAGGGTTTTCATGCCGTATAAGTCTTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATTAAAATGACGCCTTTTACAGCGCCATTTTAAACTGATGACAAACGCAACATTGCCTGATgcgc                                                                                                                      <  1:115907/151‑1 (MQ=17)
                                         aTTGTATTTTAAGGGTTTTCTTGCGGTATATGTCGTTCAAGACGGTTCGGGCATGAACAAAGCCTACGTTGTCTTCAATATTAAAATGACTCCTTTTACAGCGCCATTTTAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCtt                                                                                                               <  2:1660327/151‑1 (MQ=18)
                                         aTTGGATTTGAAGGGTTTTCATGCCGCATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTAc                                                                                                                   <  5:1390933/147‑1 (MQ=255)
                                         aTTGGATTTGAAGGGCTTTCCAGCCGGATCTGCCTTTCACGCCTGATGAGGCATGATCTTACGCTACTTTGTATTCAATATTAAAATGTCCCCTTTTACAGCGCCATTTTTAACTGATGACAAACGCAACATTGCCTGATGCGCTACg                                                                                                                  <  7:1658521/148‑1 (MQ=12)
                                          ttGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCtt                                                                                                               >  1:1382406/1‑150 (MQ=255)
                                            ggATTTGAAGGGTTTTCTTGCCGGATAAGTCTTTCACGCCGCATACGGCATGAACAAAGCCTACGTTGTCTACAATATTAAAATGACTCCTTTTACAGCGCCATTTTAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATc                                                                                                            <  4:1627553/151‑1 (MQ=18)
                                            ggATTTGAAGGGTTTTCATGCCGGATAAGTCTTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATTAAAATGACGCCTTTTACAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTAt                                                                                                             <  1:36452/150‑1 (MQ=255)
                                              aggtcaccggTTTTCACGCCGCCTAAGCCTTTCACCCCGCATCCGGCATGAACAAAGCCTACGTTTTCAACAATATTAAAATGCCGCCTTTTACAGCGCCATTTTAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAg                                                                                                          <  7:1548501/143‑1 (MQ=18)
                                              aTTTGAATGTTTTTCATTCCGTATAAGCCGTTCACGCCGCTTCGGGCATGAACAAAGCCTACTTTGTCTACAATATGAAAATGACCCCTTTTACAGCGCCATTTTAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAg                                                                                                          <  2:1490248/151‑1 (MQ=18)
                                                 tGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGcc                                                                                                       >  7:1434707/1‑151 (MQ=255)
                                                           tCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTGAAACTGATGACACACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTc                                                                                              >  3:461235/1‑150 (MQ=255)
                                                           tCATGCCGGATAAGCCGTTCACGGCGCATACTGGATGTACAAAGCCTACGTTGTCTATAATATGAAAATGACGCCTTTTACAGCGCCATTTTAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTc                                                                                              <  7:1248896/150‑1 (MQ=25)
                                                            cATGCCGTTTAAGTCGATCACGCCGCATCCGGCATGAACAAAGCCTACGTTTTATACAATATTAAAATGACTCCTTTTACAGCGCCATTTTAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCAc                                                                                            <  2:1123462/151‑1 (MQ=21)
                                                            cATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCAc                                                                                            >  3:737871/1‑151 (MQ=255)
                                                            cATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTGAAACTGATGAAAAACGCAACATTGACTGATGCGCTACGCTTATCAGGCCTACACGTTCa                                                                                             >  6:30064/1‑150 (MQ=255)
                                                             aTGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGATGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACt                                                                                           >  3:58106/1‑151 (MQ=255)
                                                                 cGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACTTTCACTGCa                                                                                        >  4:100219/1‑150 (MQ=255)
                                                                  ggATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAAt                                                                                      >  5:831252/1‑151 (MQ=255)
                                                                     tAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTTAAACTGATGACAAAAGCAACATTAAATCATGAGATAAGCTTATCATGACTAAACGTTCACTGCAAtat                                                                                    >  8:1647134/1‑148 (MQ=17)
                                                                       aGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTTAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCATGCCTACACGTTCACTGCAATCTAtt                                                                                 >  8:1058921/1‑151 (MQ=255)
                                                                           gTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATTAAAATGACGCCTTTTACAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGAt                                                                             <  7:708501/151‑1 (MQ=255)
                                                                           gTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTTAAACTGATGACAAACGCAACATTGACTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGAt                                                                             >  2:300283/1‑151 (MQ=255)
                                                                           gTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGAt                                                                             >  4:1110342/1‑151 (MQ=255)
                                                                              cACTCCGCATCCTGCATGAACAAAGCCTACTTTGTATACAATATGAAAATGCCGCCTTATACAGCGCCATTTTAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGAttt                                                                           <  7:184753/150‑1 (MQ=18)
                                                                              cACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTTAAACTGAACACAAACACAACATTACATAATGCGCAACACTTATCAGTACTACAAGTTCACTGCAATCTATTTTATTTg                                                                          >  8:1194296/1‑151 (MQ=18)
                                                                                cgccgcATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTTAAACTGATGAAAAAAACAACATTCCATGATGAGGTACGCTTATAAGGCCTACACGTTCAATGCAATATATAGCATATGa                                                                         >  2:373981/1‑150 (MQ=18)
                                                                                cgccgcATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTTAAACTGATCAAAAAAGCAACATTTACTGATGCGCTACCATTATCAGTACTACACTTACACTGCAATCTATTGCATTAGa                                                                         >  4:1270038/1‑150 (MQ=18)
                                                                                cgccgcATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTGAAACTGATGAAAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGa                                                                         >  6:104926/1‑150 (MQ=255)
                                                                                    gcATCCTGAATGAACATATCATACGTTCTCTACAATATTAAAATGACCCCTTTTAAAGCGCCATTTTAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGt                                                                    <  1:449910/151‑1 (MQ=17)
                                                                                     tattatgaaTGAACAAATCCTACTTTGTCTACAATATGAAAATGTCTCTTTTTACAGCGCCATTTTTAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGt                                                                    <  6:1543141/142‑1 (MQ=17)
                                                                                      aTCAGGCATGAACAAAGCCTACGTTGTATACAATATCAAAATGCCCCCTTTTACAGCGCCATTTTAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGttt                                                                  <  7:1001518/151‑1 (MQ=18)
                                                                                           gCATGAACAAAGCCTACGTTGTATACAATATCAAAATGACCCCTTTTACAGCGCCATTTTAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATg                                                             <  7:324834/151‑1 (MQ=255)
                                                                                            cATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTGCAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACTTTCACTGCAATATATTGGATTTGAAGGGTTTTCATgg                                                            >  1:1272410/1‑150 (MQ=255)
                                                                                            cATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTGAAACTGATGAAAAAAGCAACATTGCCTGAATCGCTACGATTATAATGACTACACGTTCAATGCAATCTAATGGATATGAAGGGTTTTCATGc                                                            >  6:606150/1‑151 (MQ=25)
                                                                                             aTGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTGAAACTGAAGACAAAAGAAACATTGCATGATGCGCTACGCTTATCAGTACTACACGTTCACTGCAATATATTGGATATGAAGGGTTTTCATGcc                                                           >  6:130236/1‑151 (MQ=21)
                                                                                                aaCAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCgg                                                         >  3:1545036/1‑150 (MQ=255)
                                                                                                       ccTACGTTGTCTACAATATGAAAATGCCGCCTTTTACAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCg                                                 <  7:477972/151‑1 (MQ=255)
                                                                                                          aCGTTGTCTACAATATGAAAATGACTCATTTTACAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTc                                              <  8:1265569/151‑1 (MQ=255)
                                                                                                          aCGTTGTAAACAATATAAAAATGTCTCATTTTAAAGCGCCATTTTAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGtt                                               <  1:1579376/150‑1 (MQ=11)
                                                                                                             ttGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCAcg                                           >  5:1567298/1‑151 (MQ=255)
                                                                                                             ttGTATACAATATTAAAATGACTCCTTTTACAGTCTTATTTTAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCAc                                            <  4:301605/150‑1 (MQ=14)
                                                                                                             gtgtattctaaatGAACATGACGCCTTTTACTGCGCTATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCAc                                            <  2:1228691/139‑1 (MQ=14)
                                                                                                              tGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCAcgc                                          >  3:639317/1‑151 (MQ=255)
                                                                                                                tattcAATATTAAAATGTCTCTTTTTATAGCGCCATTTTTAACTGATGACAAAGGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCAcgcc                                         <  2:80949/146‑1 (MQ=11)
                                                                                                                tCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCAcgccg                                        >  1:634086/1‑151 (MQ=255)
                                                                                                                 cTACAATATGAAAATGACGCCTTTTACAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCAcgccgc                                       >  2:195858/1‑151 (MQ=255)
                                                                                                                   aCAATAAAAAAATGCCCCCTTTTACAGCGCCATTTTAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCAt                                     <  1:213143/151‑1 (MQ=255)
                                                                                                                    cAATATGAAAATGACGCCTTTTACAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCCTc                                    >  3:1701895/1‑151 (MQ=255)
                                                                                                                    cAATATGAAAATGACGCCTTTTACAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATGGGATTTGAAGGGTATTCATGCCGGATAAGGCGTTCACGCCGCagc                                    >  2:552959/1‑149 (MQ=255)
                                                                                                                      atatGAAAATGACGCCTTTTACAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCg                                  >  7:922515/1‑151 (MQ=255)
                                                                                                                           aaaaTGACGCCTTTTACAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATg                             >  5:1006469/1‑151 (MQ=255)
                                                                                                                           aaaaTGACGCCTTTTACAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATg                             >  6:391298/1‑151 (MQ=255)
                                                                                                                           aaaaTGACGCCTTTTACAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATg                             >  5:789545/1‑151 (MQ=255)
                                                                                                                            aaaTGACGCCTTTTAGATCTCCATTTTAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATg                             <  2:1750545/150‑1 (MQ=255)
                                                                                                                             aatgtctcattttACAGCTCCATTTTTAACTTATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTAGACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGaa                           <  4:353735/142‑1 (MQ=11)
                                                                                                                             aaTGACGCCTTTTACAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGa                            >  8:1535542/1‑150 (MQ=255)
                                                                                                                              aggacccctggaaCAGCCCCATCTCAAACTGCACACCACCGCACCCTTGCCCGATGCCCAACTCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAAc                          <  5:1094269/139‑1 (MQ=11)
                                                                                                                              aTGACGCCTTTTACAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAAc                          >  3:295546/1‑151 (MQ=255)
                                                                                                                              aTGACGCCTTTTACAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAAc                          >  5:1775580/1‑151 (MQ=255)
                                                                                                                               tGACGCCTTTTACAGCGACATTTGAAACTGATGACAAACGCAACATTGCCTGATTCGCTACGCTTATCAGGCCTAGACGTTCACTGAAATCTATTGGATTTGAAGGGTTTTCATGACGGATAAGGCGTTCACGCCGCATCCGGCATCAAc                          >  6:1151924/1‑150 (MQ=14)
                                                                                                                                 ccGCCTTATACAGCGCCATTTAAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACaaa                       <  5:413629/150‑1 (MQ=17)
                                                                                                                                 accccTTTTACAGCGCCATTTGAAACTGAGGACAAACGCACCATTGCCTGATGCGCCACGCCTCACAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACaaa                       <  1:1730556/148‑1 (MQ=14)
                                                                                                                                       tttACAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCt                   <  3:1271174/149‑1 (MQ=38)
                                                                                                                                         tACAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGACGCATCCGGCATGAACAAACCCTACGt               >  8:757935/1‑151 (MQ=25)
                                                                                                                                          aCAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGGTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGtt              <  3:944376/151‑1 (MQ=25)
                                                                                                                                             gcgcCATTTTAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTc           <  1:805143/151‑1 (MQ=14)
                                                                                                                                             gcgcCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTc           >  1:1005445/1‑151 (MQ=33)
                                                                                                                                                tttttttAAACTTATGACTTCCGCTATATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTa         <  4:715382/147‑1 (MQ=11)
                                                                                                                                                ccATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTAc        >  7:1143484/1‑151 (MQ=33)
                                                                                                                                                  aTTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACa       >  5:1608615/1‑150 (MQ=33)
                                                                                                                                                   tttGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAat     >  6:808194/1‑151 (MQ=33)
                                                                                                                                                    ttGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAat     <  5:371458/150‑1 (MQ=32)
                                                                                                                                                      gAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGTCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATg  >  5:536986/1‑151 (MQ=14)
                                                                                                                                                      gAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAatat   <  7:238099/150‑1 (MQ=31)
                                                                                                                                                      gAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATg  >  7:60526/1‑151 (MQ=31)
                                                                                                                                                      |                                                                                                                                                      
GATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATGAAAATGACGCCTTTTACAGCGCCATTTGAAACTGATGACAAACGCAACATTGCCTGATGCGCTACGCTTATCAGGCCTACACGTTCACTGCAATCTATTGGATTTGAAGGGTTTTCATGCCGGATAAGGCGTTCACGCCGCATCCGGCATGAACAAAGCCTACGTTGTCTACAATATG  >  DH1/2326137‑2326437

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: