Predicted mutation
evidence seq id position mutation annotation gene description
RA DH1 4,204,741 C→A H381N (CAT→AAT)  alsA → D‑allose import ATP‑binding protein AlsA

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*DH14,204,7410CA98.8% 2189.1 / 12.2 586H381N (CAT→AAT) alsAD‑allose import ATP‑binding protein AlsA
Reads supporting (aligned to +/- strand):  ref base C (2/5);  new base A (300/279);  total (302/284)
Fisher's exact test for biased strand distribution p-value = 2.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.36e-01

Only 100 of 586 total aligned reads displayed.

ACCCCTGGATGCCGTGAAAAAAGGGATGGCTTACATCACTGAAAGCCGCCGGGATAACGGTTTTTTCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTCATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTCCTGATCTCCAAATGGCTGTGCTGTTGCCCGGAA  >  DH1/4204593‑4204887
                                                                                                                                                    |                                                                                                                                                  
 ccccTGGATGCCGTGAAAAAAGGGATGGCTTACATCACTGAAAGCCGCCGGGATAACGGTTTTTTCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATg                                                                                                                                                 <  1:799116/151‑1 (MQ=255)
      ggATGCCGTGAAAAAAGGGATGGCTTACATCACTGAAAGCCGCCGGGATAACGGTTTTTTCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGtt                                                                                                                                            >  2:1432503/1‑151 (MQ=255)
      ggATGCCGTGAAAAAAGGGATGGCTTACATCACTGAAAGCCGCCGGGATAACGGTTTTTTCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGt                                                                                                                                             <  7:1772883/150‑1 (MQ=255)
        aTGCCGTGAAAAAAGGGATGGCTTACATCACTGAAAGCCGCCGGGATAACGGTTTTTTCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGa                                                                                                                                          <  1:1619188/151‑1 (MQ=255)
         tGCCGTGAAAAAAGGGATGGCTTACATCACTGAAAGCCGCCGGGATAACGGTTTTTTCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGAc                                                                                                                                         >  8:131976/1‑151 (MQ=255)
          gCCGTGAAAAAAGGGATGGCTTACATCACTGAAAGCCGCCGGGATAACGGTTTTTTCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGAc                                                                                                                                         <  3:937011/150‑1 (MQ=255)
            cGTGAAAAAAGGGATGTCTTACATCACTGAAAGCCGCCGGGATAACGGTTTTTTCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAg                                                                                                                                      <  5:1826178/151‑1 (MQ=255)
            cGTGAAAAAAGGGATGGCTTACATCACTGAAAGCCGCCGGGATAACGGTTTTTTCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACg                                                                                                                                        >  1:1717777/1‑149 (MQ=255)
               gAAAAAAGGGATGGCTTACATCACTGAAAGCCGCCGGGATAACGGTTTTTTCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCa                                                                                                                                    <  8:1066613/150‑1 (MQ=255)
                  aaaaGGGATGGCTTACATCACTGAAAGCCGCCGGGATAACGGTTTTTTCCCCAACTTTTCCATCCCTCAGAACATGGCGATCAGCCGCAGTCTCAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGATGACCAGCAACGt                                                                                                                                >  1:373580/1‑151 (MQ=255)
                    aaGGGATGGCTTACATCACTGAAAGCCGCCGGGATAACGGTTTTTTCCCCAACTTTTCCATCGCTCGGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTAc                                                                                                                              >  7:73698/1‑151 (MQ=255)
                    aaGGGATGGCTTACATCACTGAAAGCCGCCGGGATAACGGTTTTTTCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTAc                                                                                                                              >  4:1900276/1‑151 (MQ=255)
                    aaGGGATGGCTTACATCACTGAAAGCCGCCGGGATAACGGTTTTTTCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTAc                                                                                                                              >  2:1886938/1‑151 (MQ=255)
                        gATGGCTTACATCACTGAAAGCCGCCGGGATAACGGTTTTTTCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCt                                                                                                                          <  3:681196/151‑1 (MQ=255)
                          tGGCTTACATCACTGAAAGCCGCCGGGATAACGGTTTTTTCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTg                                                                                                                         >  1:1018943/1‑150 (MQ=255)
                          tGGCTTACATCACTGAAAGCCGCCGGGATAACGGTTTTTTCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTg                                                                                                                         >  7:1003662/1‑150 (MQ=255)
                             cTTACATCACTGAAAGCCGCCGGGATAACGGTTTTTTCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGCCGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGaaa                                                                                                                      >  6:1460283/1‑150 (MQ=255)
                                 cATCACTGAAAGCCGCCGGGATAACGGTTTTTTCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCaa                                                                                                                 <  8:540701/151‑1 (MQ=255)
                                  aTCACTGAAAGCCGCCGGGATAACGGTTTTTTCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAAc                                                                                                                <  1:1855319/151‑1 (MQ=255)
                                    cACTGAAATCCGCCGGGATAACGGTTTTTTGCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAATACTGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAAcgc                                                                                                              <  6:330365/151‑1 (MQ=255)
                                    cACTGAAAGCCTCCGGGATAACGGTTTTTTCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAAcg                                                                                                               <  3:1654053/150‑1 (MQ=255)
                                          aaGCCGCCGGGATAACGGTTTTTTCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAActg                                                                                                        <  8:143303/151‑1 (MQ=255)
                                             ccgccgGGATAACGGTTTTTTCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAActgctg                                                                                                     <  5:907802/151‑1 (MQ=255)
                                              cgccgGGATAACGGTTTTTTCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTgg                                                                                                    <  7:610792/151‑1 (MQ=255)
                                               gccgGGATAACGGTTTTTTCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTTAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTTAAAATCAACGCGAACTGCTGgc                                                                                                   <  1:1360701/151‑1 (MQ=255)
                                               gccgGGATAACGGTTTTTTCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGgc                                                                                                   >  5:1321488/1‑151 (MQ=255)
                                                ccgGGATAACGGTTTTTTCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGgcg                                                                                                  <  4:27241/151‑1 (MQ=255)
                                                ccgGGATAACGGTTTTTTCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGgcg                                                                                                  >  4:1261613/1‑151 (MQ=255)
                                                ccgGGATAACGGTTTTTTCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGgcg                                                                                                  <  2:1172942/151‑1 (MQ=255)
                                                      tAACGGTTTTTTCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGaa                                                                                             >  3:88917/1‑150 (MQ=255)
                                                       aaCGGTTTTTTCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGCACTGCTGGCGCTGAAAt                                                                                           >  8:1178147/1‑151 (MQ=255)
                                                        aCGGTTTTTTCCCGAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTAGCGCTGAAAATCAACGCGAACTGCTGGCTCTGAAATg                                                                                          <  3:345926/151‑1 (MQ=255)
                                                           gTTTTTTCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTc                                                                                        >  4:1455881/1‑150 (MQ=255)
                                                            ttttttCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCAt                                                                                      >  7:928452/1‑151 (MQ=255)
                                                              ttttCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCAtt                                                                                     <  2:167821/150‑1 (MQ=255)
                                                                  ccccAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTa                                                                                <  3:1055535/151‑1 (MQ=255)
                                                                    ccAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTaa                                                                               >  3:648503/1‑150 (MQ=255)
                                                                     cAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCACTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAAc                                                                             >  2:1517144/1‑151 (MQ=255)
                                                                        cTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAAcc                                                                            >  7:18666/1‑149 (MQ=255)
                                                                           ttCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGaa                                                                        <  4:1559202/150‑1 (MQ=255)
                                                                            tCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAat                                                                       <  2:1259628/150‑1 (MQ=255)
                                                                             ccATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAatat                                                                     >  1:509270/1‑151 (MQ=255)
                                                                             ccATCGCTCAGAACATGGCGAGCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAatat                                                                     <  6:535540/151‑1 (MQ=255)
                                                                              cATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATc                                                                    >  3:1278006/1‑151 (MQ=255)
                                                                                 cGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTCATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCAcc                                                                 >  5:1957487/1‑151 (MQ=255)
                                                                                    tCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGa                                                               >  4:1672515/1‑150 (MQ=255)
                                                                                     cAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGaa                                                              <  6:770725/150‑1 (MQ=255)
                                                                                     cAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAAc                                                             >  7:339000/1‑151 (MQ=255)
                                                                                      aGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTTATTCGGTAAACCAGAATATCACCGAAct                                                            <  1:1233787/151‑1 (MQ=255)
                                                                                      aGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAAct                                                            >  8:1909143/1‑151 (MQ=255)
                                                                                       gAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAActc                                                           >  2:1089216/1‑151 (MQ=255)
                                                                                        aaCATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAActct                                                          >  7:158124/1‑151 (MQ=255)
                                                                                          cATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTcc                                                        <  5:28872/151‑1 (MQ=255)
                                                                                           aTGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCg                                                       >  7:1962051/1‑151 (MQ=255)
                                                                                           aTGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCg                                                       >  3:1420300/1‑151 (MQ=255)
                                                                                           aTGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCg                                                       <  2:662868/151‑1 (MQ=255)
                                                                                               cGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCggggg                                                   <  8:787858/151‑1 (MQ=255)
                                                                                                gATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGa                                                  <  5:1528561/151‑1 (MQ=255)
                                                                                                   cAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAAt                                               >  3:250246/1‑151 (MQ=255)
                                                                                                      ccGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGGGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATcag                                            <  1:114399/151‑1 (MQ=255)
                                                                                                      ccGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATcag                                            <  1:668007/151‑1 (MQ=255)
                                                                                                        gCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATcaga                                           >  1:959632/1‑149 (MQ=255)
                                                                                                         cAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATcagcag                                         <  4:536907/151‑1 (MQ=255)
                                                                                                          aGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATcagc                                           >  6:1210056/1‑148 (MQ=255)
                                                                                                              tGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGt                                    >  5:1010200/1‑151 (MQ=255)
                                                                                                              tGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGt                                    >  6:1898509/1‑151 (MQ=255)
                                                                                                               gAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTc                                   <  8:867190/151‑1 (MQ=255)
                                                                                                                aaaGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTcc                                  >  2:719546/1‑151 (MQ=255)
                                                                                                                aaaGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTcc                                  <  4:917495/151‑1 (MQ=255)
                                                                                                                 aaGACGGCGGCTATAAAGGCGCGATGGGCTTCTTTAATGAAGTTGACGAGCAACGTACCCCAGAAAATCAACGCGAACTGCTGGCCCTGAAATGTCATTCAGTAAACAAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTccg                                 >  4:77682/1‑150 (MQ=255)
                                                                                                                   gACGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGATCTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTCCTGa                               <  8:363698/151‑1 (MQ=255)
                                                                                                                    aCGGCGGCTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTCCTGAt                              >  3:618644/1‑151 (MQ=255)
                                                                                                                      ggcggcTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTCCTGAtct                            <  5:614150/151‑1 (MQ=255)
                                                                                                                      ggcggcTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTCCTGAtct                            <  7:1279111/151‑1 (MQ=255)
                                                                                                                       gcggcTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTCCTGAtctc                           <  4:1752523/151‑1 (MQ=255)
                                                                                                                        cggcTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTCCTGATCTcc                          <  2:862579/151‑1 (MQ=255)
                                                                                                                        cggcTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTCCTGATCTcc                          >  2:322878/1‑151 (MQ=255)
                                                                                                                          gcTATAAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTCCTCATCTCCa                         >  5:1639983/1‑150 (MQ=255)
                                                                                                                             ataAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTCCTGATCTCCAAAt                      <  2:220688/150‑1 (MQ=255)
                                                                                                                             ataAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTCCTGATCTCCAAATg                     <  6:693072/151‑1 (MQ=255)
                                                                                                                             ataAAGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTCCTGATCTCCAAATg                     <  6:261252/151‑1 (MQ=255)
                                                                                                                               aaaGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTCCTGATCTCCAAATgg                    <  4:672937/150‑1 (MQ=255)
                                                                                                                               aaaGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTCCTGATCTCCAAATGgc                   >  8:1784005/1‑151 (MQ=255)
                                                                                                                                aaGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTCCTGATCTCCAAATGgct                  >  7:1101360/1‑151 (MQ=255)
                                                                                                                                 aGGCGCGATGGGCTTGTTTAATGAAGTTGACGAGCATCGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTCCTGATCTCCAAATGgctg                 <  2:51001/151‑1 (MQ=255)
                                                                                                                                  ggCGCGATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTCCTGATTTCCAAATGgctgt                <  4:204907/151‑1 (MQ=255)
                                                                                                                                     gcgATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTCCTGATCTCCAAATGgctgtgct             <  6:449742/151‑1 (MQ=255)
                                                                                                                                       gATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTCCTGATCTCCAAATGgctgtgctgt           >  7:1327296/1‑151 (MQ=255)
                                                                                                                                       gATGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTCCTGATCTCCAAATGgctgtgctgt           >  5:457145/1‑151 (MQ=255)
                                                                                                                                        aTGGGCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTCCTGATCTCCAAATGGCTGTGCTGtt          >  6:1747364/1‑151 (MQ=255)
                                                                                                                                          gggCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTCCTGATCTCCAAATGGCTGTGCTGTTg         <  7:1456374/150‑1 (MQ=255)
                                                                                                                                          gggCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTCCTGATCTCCAAATGGCTGTGCTGTTg         <  6:855384/150‑1 (MQ=255)
                                                                                                                                            gCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTCCTGATCTCCAAATGGCTGTGCTGTTGccc      <  3:150849/151‑1 (MQ=255)
                                                                                                                                            gCTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTCCTGATCTCCAAATGGCTGTGCTGTTGccc      >  1:21687/1‑151 (MQ=255)
                                                                                                                                             cTTGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAACCCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTCCTGATCTCCAAATGGCTGTGCTGTTGcccc     >  3:184476/1‑150 (MQ=255)
                                                                                                                                               tGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTTGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGTGGAAATCAGCAGAAAGTCCTGATCTCCAAATGGCTGTGCTGTTGCCCg     <  4:1074228/149‑1 (MQ=255)
                                                                                                                                               tGTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTCCTGATCTCCAAATGGCTGTGCTGTTGCCCg     <  5:1730331/149‑1 (MQ=255)
                                                                                                                                                gTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTCCTGATCTCCAAATGGCTGTGCTGTTGCCCGGa   >  3:175415/1‑150 (MQ=255)
                                                                                                                                                gTTTAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTCCTGATCTCCAAATGGCTGTGCTGTTGCCCGGa   <  5:75774/150‑1 (MQ=255)
                                                                                                                                                 tttAATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTCCTGATCTCCAAATGGCTGTGCTGTTGCCCGGa   >  5:1199504/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                  
ACCCCTGGATGCCGTGAAAAAAGGGATGGCTTACATCACTGAAAGCCGCCGGGATAACGGTTTTTTCCCCAACTTTTCCATCGCTCAGAACATGGCGATCAGCCGCAGTCTGAAAGACGGCGGCTATAAAGGCGCGATGGGCTTGTTTCATGAAGTTGACGAGCAACGTACCGCTGAAAATCAACGCGAACTGCTGGCGCTGAAATGTCATTCGGTAAACCAGAATATCACCGAACTCTCCGGGGGAAATCAGCAGAAAGTCCTGATCTCCAAATGGCTGTGCTGTTGCCCGGAA  >  DH1/4204593‑4204887

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