Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,701,108 A→T 12.9% N51I (AAT→ATT)  malY → PLP‑dependent beta‑cystathionase and maltose regulon regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,701,1080AT12.9% 372.5 / 39.8 124N51I (AAT→ATT) malYPLP‑dependent beta‑cystathionase and maltose regulon regulator
Reads supporting (aligned to +/- strand):  ref base A (55/53);  new base T (10/6);  total (65/59)
Fisher's exact test for biased strand distribution p-value = 4.32e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.94e-01

Only 100 of 130 total aligned reads displayed.

TTCGATTTTTCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTT  >  NC_000913/1700960‑1701256
                                                                                                                                                    |                                                                                                                                                    
ttCGATTTTTCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGaa                                                                                                                                                      <  1:212290/149‑1 (MQ=255)
ttCGATTTTTCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGaa                                                                                                                                                      >  4:521208/1‑149 (MQ=255)
ttCGATTTTTCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGaa                                                                                                                                                      >  4:261574/1‑149 (MQ=255)
ttCGATTTTTCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGaa                                                                                                                                                      >  4:243687/1‑149 (MQ=255)
ttCGATTTTTCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGaa                                                                                                                                                      >  1:767848/1‑149 (MQ=255)
 tCGATTTTTCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAAt                                                                                                                                                     <  4:368125/149‑1 (MQ=255)
  cGATTTTTCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATc                                                                                                                                                    >  4:872598/1‑149 (MQ=255)
     tttttCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAgc                                                                                                                                                 >  2:436244/1‑149 (MQ=255)
     tttttCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGAGGGCTTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAgc                                                                                                                                                 >  3:210839/1‑149 (MQ=255)
        ttCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGATTCAGCGcc                                                                                                                                              <  4:174492/149‑1 (MQ=255)
           aaaGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGATTCAGCGCCtc                                                                                                                                            >  3:1163/1‑147 (MQ=255)
           aaaGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGATTCAGCGCCtc                                                                                                                                            <  4:1163/148‑2 (MQ=255)
            aaGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGAt                                                                                                                                          <  2:26738/149‑1 (MQ=255)
                     gATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGt                                                                                                                                 >  1:808035/1‑149 (MQ=255)
                      aTCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTa                                                                                                                                <  4:468319/149‑1 (MQ=255)
                       tCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTAt                                                                                                                               >  2:853490/1‑149 (MQ=255)
                       tCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTAt                                                                                                                               <  2:162040/149‑1 (MQ=255)
                          tCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTg                                                                                                                            >  4:772042/1‑149 (MQ=255)
                              ggCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTTACCTGTTACCGTTCACGATTCCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTa                                                                                                                        >  2:800775/1‑149 (MQ=255)
                              ggCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTa                                                                                                                        >  1:354827/1‑149 (MQ=255)
                                    tGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCg                                                                                                                  >  1:893606/1‑149 (MQ=255)
                                      gtgtACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCt                                                                                                                <  2:348687/149‑1 (MQ=255)
                                          acacAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGaa                                                                                                            >  2:223482/1‑149 (MQ=255)
                                             caGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGaaaaa                                                                                                         <  1:821225/149‑1 (MQ=255)
                                              aGTGGGATTATGTCGCTGACCGTTTCGGCACCGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAAc                                                                                                        <  1:439992/149‑1 (MQ=255)
                                                tGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAAc                                                                                                        >  1:809/1‑147 (MQ=255)
                                                tGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAAc                                                                                                        <  2:809/147‑1 (MQ=255)
                                                tGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGa                                                                                                      <  3:872598/149‑1 (MQ=255)
                                                tGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGa                                                                                                      <  1:367733/149‑1 (MQ=255)
                                                 gggATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGAt                                                                                                     <  1:555940/149‑1 (MQ=255)
                                                  ggATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATg                                                                                                    <  1:853490/149‑1 (MQ=255)
                                                  ggATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATg                                                                                                    >  4:183568/1‑149 (MQ=255)
                                                   gATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGa                                                                                                   <  2:893606/149‑1 (MQ=255)
                                                   gATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGa                                                                                                   <  4:495308/149‑1 (MQ=255)
                                                     ttATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTg                                                                                                                                            <  1:634808/106‑1 (MQ=255)
                                                     ttATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTg                                                                                                                                            >  2:634808/1‑106 (MQ=255)
                                                     ttATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGt                                                                                                 <  4:494747/149‑1 (MQ=255)
                                                     ttATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGt                                                                                                 >  2:21870/1‑149 (MQ=255)
                                                              tGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCgg                                                                                        <  2:651520/149‑1 (MQ=255)
                                                                   gTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTAtt                                                                                   >  2:744995/1‑149 (MQ=255)
                                                                    tttCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTg                                                                                  <  3:879140/149‑1 (MQ=255)
                                                                      tCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGcc                                                                                <  4:210839/149‑1 (MQ=255)
                                                                           aCTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTg                                                                           <  3:521208/149‑1 (MQ=255)
                                                                               ctgACCTGTTACCGTTCACGCTTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAAACAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGttt                                                                       <  3:159541/149‑1 (MQ=255)
                                                                                    cTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCAc                                                                  >  4:348684/1‑149 (MQ=255)
                                                                                       ttACCGTTCACGATTTCCGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATACACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCa                                                               <  1:707708/149‑1 (MQ=255)
                                                                                        tACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAg                                                              <  2:831297/149‑1 (MQ=255)
                                                                                           cGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCAt                                                           <  2:593585/149‑1 (MQ=255)
                                                                                           cGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCAt                                                           >  2:239202/1‑149 (MQ=255)
                                                                                            gTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCAtt                                                          >  3:486173/1‑149 (MQ=255)
                                                                                             ttCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTa                                                         >  2:484562/1‑149 (MQ=255)
                                                                                             ttCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTa                                                         >  3:471982/1‑149 (MQ=255)
                                                                                                  gATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCg                                                    >  4:711796/1‑149 (MQ=255)
                                                                                                  gATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCg                                                    >  2:909001/1‑149 (MQ=255)
                                                                                                   aTTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGc                                                   >  4:766241/1‑149 (MQ=255)
                                                                                                    tttCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGcc                                                  >  3:119391/1‑149 (MQ=255)
                                                                                                    tttCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGcc                                                  >  3:899119/1‑149 (MQ=255)
                                                                                                    tttCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGcc                                                  >  2:549774/1‑149 (MQ=255)
                                                                                                    tttCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGcc                                                  >  1:519778/1‑149 (MQ=255)
                                                                                                              ggATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATtctc                                        >  2:428188/1‑149 (MQ=255)
                                                                                                                aTTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAg                                      >  3:193018/1‑149 (MQ=255)
                                                                                                                aTTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAg                                      >  1:629699/1‑149 (MQ=255)
                                                                                                                 ttttGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGttt                                                                       <  4:538319/115‑1 (MQ=255)
                                                                                                                 ttttGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGttt                                                                       >  3:538319/1‑115 (MQ=255)
                                                                                                                 ttttGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGa                                     <  2:629699/149‑1 (MQ=255)
                                                                                                                  tttGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGAc                                    <  2:767848/149‑1 (MQ=255)
                                                                                                                     gCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACggt                                 <  2:124730/149‑1 (MQ=255)
                                                                                                                      ccACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACggtg                                <  3:183568/149‑1 (MQ=255)
                                                                                                                       cACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACggtgg                               >  4:191837/1‑149 (MQ=255)
                                                                                                                             ccccTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATg                         <  1:428188/149‑1 (MQ=255)
                                                                                                                             ccccTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATg                         <  3:711796/149‑1 (MQ=255)
                                                                                                                             ccccTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCACGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATg                         >  4:548437/1‑149 (MQ=255)
                                                                                                                              cccTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATgg                        >  1:900829/1‑149 (MQ=255)
                                                                                                                              cccTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATg                         >  2:163658/1‑148 (MQ=255)
                                                                                                                              cccTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATg                         <  1:163658/148‑1 (MQ=255)
                                                                                                                               ccTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGc                       >  1:697316/1‑149 (MQ=255)
                                                                                                                               ccTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGc                       >  2:656404/1‑149 (MQ=255)
                                                                                                                                cTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGcc                      >  4:61555/1‑149 (MQ=255)
                                                                                                                                cTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGcc                      <  4:287605/149‑1 (MQ=255)
                                                                                                                                   cATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCtt                   >  4:849701/1‑149 (MQ=255)
                                                                                                                                    aTTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATTAGTTTCTCGCGGCTATTGCCCACTGGTGTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCTTTc                  >  1:925305/1‑149 (MQ=255)
                                                                                                                                    aTTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTc                  <  4:133943/149‑1 (MQ=255)
                                                                                                                                    aTTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATGTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTc                  >  4:461775/1‑149 (MQ=255)
                                                                                                                                     ttATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCt                 <  3:471433/149‑1 (MQ=255)
                                                                                                                                     ttATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCt                 <  2:637878/149‑1 (MQ=255)
                                                                                                                                     tgcatcaGGCGCTGATTCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCt                 >  3:288462/7‑149 (MQ=255)
                                                                                                                                       aTCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGt               <  3:653984/149‑1 (MQ=255)
                                                                                                                                        atcaGGCGCTGATTCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTc              >  1:475402/4‑149 (MQ=255)
                                                                                                                                        atcaGGCGCTGATTCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTc              >  3:53655/4‑149 (MQ=255)
                                                                                                                                         tcaGGCGCTGATTCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCa             >  3:834186/3‑149 (MQ=255)
                                                                                                                                         tcaGGCGCTGATTCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCa             >  4:551107/3‑149 (MQ=255)
                                                                                                                                          caGGCGCTGATTCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACggtg                                >  1:360316/2‑129 (MQ=255)
                                                                                                                                          caGGCGCTGATTCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACggtg                                <  2:360316/128‑1 (MQ=255)
                                                                                                                                           aGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATc           <  3:243356/149‑1 (MQ=255)
                                                                                                                                            ggCGCTGATTCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCt          >  3:582960/1‑149 (MQ=255)
                                                                                                                                            ggCGCTGATTCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCt          >  2:262558/1‑149 (MQ=255)
                                                                                                                                               gcTGATTCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCtata       <  1:262558/149‑1 (MQ=255)
                                                                                                                                                cTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTatat      <  1:21870/149‑1 (MQ=255)
                                                                                                                                                cTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTatat      >  4:804819/1‑149 (MQ=255)
                                                                                                                                                    aTCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGtt  >  3:39404/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TTCGATTTTTCAAAGGTCGTGGATCGTCATGGCACATGGTGTACACAGTGGGATTATGTCGCTGACCGTTTCGGCACTGCTGACCTGTTACCGTTCACGATTTCAGACATGGATTTTGCCACTGCCCCCTGCATTATCGAGGCGCTGAATCAGCGCCTGATGCACGGCGTATTTGGCTACAGCCGCTGGAAAAACGATGAGTTTCTCGCGGCTATTGCCCACTGGTTTTCCACCCAGCATTACACCGCCATCGATTCTCAGACGGTGGTGTATGGCCCTTCTGTCATCTATATGGTT  >  NC_000913/1700960‑1701256

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: