Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,309,833 A→G V383A (GTT→GCT)  alsA ← D‑allose ABC transporter ATP binding subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,309,8330AG97.8% 328.9 / ‑2.9 92V383A (GTT→GCT) alsAD‑allose ABC transporter ATP binding subunit
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (47/43);  total (48/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.92e-01

CCGGGCAACAGCACAGCCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTT  >  NC_000913/4309696‑4309968
                                                                                                                                         |                                                                                                                                       
ccGGGCAACAGCACAGCCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCt                                                                                                                                       >  2:985142/1‑140 (MQ=255)
ccGGGCAACAGCACAGCCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCt                                                                                                                                       >  2:159965/1‑140 (MQ=255)
   ggCAACAGCACAGCCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCa                                                                                                                                    <  1:1639296/140‑1 (MQ=255)
    gCAACAGCACAGCCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCAt                                                                                                                                   <  1:1871381/140‑1 (MQ=255)
        cagcacagcCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCa                                                                                                                                    <  1:1691179/135‑1 (MQ=255)
        cagcacagcCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCa                                                                                                                                    >  2:1691179/1‑135 (MQ=255)
        cagcacagcCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGaaa                                                                                                                               >  1:1699517/1‑140 (MQ=255)
        cagcacagcCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGaaa                                                                                                                               >  1:1638872/1‑140 (MQ=255)
        cagcacagcCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGaaa                                                                                                                               >  2:1763461/1‑140 (MQ=255)
        cagcacagcCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGaaa                                                                                                                               >  2:973436/1‑140 (MQ=255)
            acagcCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAg                                                                                                                           <  1:1763461/140‑1 (MQ=255)
                   tttGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATcg                                                                                                                    <  1:1456937/140‑1 (MQ=255)
                   tttGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATcg                                                                                                                    <  2:1699517/140‑1 (MQ=255)
                   tttGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCCGTACGTTACTCCTCAGCTTCACGAAACAAGCCCATcg                                                                                                                    >  1:838647/1‑140 (MQ=255)
                    ttGGAGATCAGGACTTTATGCTTATTTCCCTCGTAAAGTCCGGTGATAATCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATcgc                                                                                                                   <  2:932461/140‑1 (MQ=255)
                     tGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAg                                                                                                                           <  2:1585756/131‑1 (MQ=255)
                     tGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAg                                                                                                                           >  1:1585756/1‑131 (MQ=255)
                     tGGAGATCAGGACTTTCTGCTGATATCCCCCAGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATcgcg                                                                                                                  <  1:1075686/140‑1 (MQ=255)
                           tCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTta                                                                                                            >  2:1322985/1‑140 (MQ=255)
                            cAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTtat                                                                                                           >  2:238612/1‑140 (MQ=255)
                                aCTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAgcc                                                                                                       >  2:1383624/1‑140 (MQ=255)
                                  tttCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAgccgc                                                                                                     >  1:1466973/1‑140 (MQ=255)
                                   ttCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAgccgcc                                                                                                    <  2:277423/140‑1 (MQ=255)
                                   ttCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAgccgcc                                                                                                    <  1:1383624/140‑1 (MQ=255)
                                   ttCTGCTGATTTCCCCCGGAGAGTTAGGTGATATTCTCGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAgccgcc                                                                                                    <  1:817059/140‑1 (MQ=255)
                                   ttCTGCTGATTTCCCACGGAGAATTAGGTAATATTCTCGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAgccgcc                                                                                                    <  1:928681/140‑1 (MQ=255)
                                    tCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAgccgcc                                                                                                    <  1:1118899/139‑1 (MQ=255)
                                    tCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAgccgcc                                                                                                    >  2:1118899/1‑139 (MQ=255)
                                    tCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGATTCATGAGACAGGCACATCGCACCTTTATAGccgccg                                                                                                   >  1:485697/1‑140 (MQ=255)
                                      tgctgATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTc                                                                                                 >  2:1758898/1‑140 (MQ=255)
                                         tgATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCttt                                                                                              <  2:1638872/140‑1 (MQ=255)
                                         tgATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCttt                                                                                              >  2:1459510/1‑140 (MQ=255)
                                           aTTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAgccgcc                                                                                                    <  1:1366021/132‑1 (MQ=255)
                                           aTTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAgccgcc                                                                                                    >  2:1366021/1‑132 (MQ=255)
                                                ccccGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTg                                                                                       >  2:1559317/1‑140 (MQ=255)
                                                   cGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCgg                                                                                    <  1:1241780/140‑1 (MQ=255)
                                                   cGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCgg                                                                                    <  1:230720/140‑1 (MQ=255)
                                                   cGGAGAGTTCGGTGAAATTCTGTTTGACCTAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCgg                                                                                    <  1:710069/140‑1 (MQ=255)
                                                    ggAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGc                                                                                   <  1:973436/140‑1 (MQ=255)
                                                     gagaGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATTAAACAAGCGCATCGCGCCTTTATAGGCGCCATCTTTCAGACTGCGGCt                                                                                  >  2:645637/1‑140 (MQ=255)
                                                     gagaGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCt                                                                                  >  2:8035/1‑140 (MQ=255)
                                                     gagaGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCt                                                                                  <  1:1044852/140‑1 (MQ=255)
                                                       gagTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGa                                                                                <  1:380241/140‑1 (MQ=255)
                                                       gagTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGa                                                                                <  1:238612/140‑1 (MQ=255)
                                                          ttCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCg                                                                             >  1:308176/1‑140 (MQ=255)
                                                          ttCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCg                                                                             >  1:275322/1‑140 (MQ=255)
                                                          ttCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCg                                                                             >  2:1491753/1‑140 (MQ=255)
                                                          ttCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCg                                                                             >  2:406591/1‑140 (MQ=255)
                                                            cGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGcc                                                                           <  1:1758898/140‑1 (MQ=255)
                                                            cGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGcc                                                                           <  2:1899494/140‑1 (MQ=255)
                                                            cGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGcc                                                                           <  1:985142/140‑1 (MQ=255)
                                                            cGGTGATATTCTGGTTCACCGAATGACATTTCAACGCCAGCAGTTCGCTTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGcc                                                                           <  2:838647/140‑1 (MQ=255)
                                                            aggTGATATTCTGGTGTACCGACTGACATGTCAGCGCTAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGcc                                                                           <  1:904975/139‑1 (MQ=255)
                                                               tGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCGATCGCGCCGTTAAAGCTGCCGTCTTTCAGAGTGCGGCTGATCGCCATg                                                                        >  2:556542/1‑140 (MQ=255)
                                                                    ttCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTCTCAGACTGCGGCTGATCGCCATGTTCTg                                                                   >  2:244707/1‑140 (MQ=255)
                                                                        ggTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCg                                                               >  1:464165/1‑140 (MQ=255)
                                                                         gTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGGCGTATCTGTGACTGCGGCTGAGCGCCATGTACTGAGcct                                                              >  2:904975/1‑138 (MQ=255)
                                                                         gTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGa                                                              >  1:1597500/1‑140 (MQ=255)
                                                                            tACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAg                                                                 >  1:955967/1‑134 (MQ=255)
                                                                            tACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAg                                                                 <  2:955967/134‑1 (MQ=255)
                                                                              ccGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGACTGATCGCCATGTTCTGAGCGATGGaa                                                         <  1:1877444/140‑1 (MQ=255)
                                                                              ccGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGACTGATCGCCATGTTCTGAGCGATGGaa                                                         >  2:1877444/1‑140 (MQ=255)
                                                                              ccGAATGACATTTCAGCGCCAGAAGTTCGCGCTGATTTTTAGCGGTACTTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGaa                                                         <  2:692729/140‑1 (MQ=255)
                                                                                   tGACATTTCAGCGCCCGCAGTTCGCCTTGCTTTTCAGTGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGtt                                                    <  1:556542/140‑1 (MQ=255)
                                                                                      cATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTggg                                                 >  1:1731188/1‑140 (MQ=255)
                                                                                      cATTTCAGCGCCAGCAGCTCGCTTTGATTTTAAGCTGTACGTTGCTCGCCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTggg                                                 <  2:532658/140‑1 (MQ=255)
                                                                                         ttCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGATGATCGCCATGTTCTGAGCGAGGGAAAAGTTGGGGaa                                              >  1:692729/1‑140 (MQ=255)
                                                                                           cAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGaaaa                                            <  1:1491753/140‑1 (MQ=255)
                                                                                           cAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGaaaa                                            <  2:312134/140‑1 (MQ=255)
                                                                                              cgcCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATAGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAAc                                         >  1:758015/1‑140 (MQ=255)
                                                                                                ccAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCg                                       <  1:1559317/140‑1 (MQ=255)
                                                                                                 cagcagTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGt                                      >  1:951851/1‑140 (MQ=255)
                                                                                                  agcagTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGtt                                     >  2:1698279/1‑140 (MQ=255)
                                                                                                   gcagTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTAGGGAAAAAACCGTTa                                    >  1:700708/1‑140 (MQ=255)
                                                                                                        tCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCcg                               <  2:275322/140‑1 (MQ=255)
                                                                                                        tCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCcg                               <  2:1466973/140‑1 (MQ=255)
                                                                                                         cgcgTTGATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCcgg                              >  2:183894/1‑140 (MQ=255)
                                                                                                           cgTTGATTTTCAGCGGTACGTTCCTCGTCAGATTCCTGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCcggcg                            <  2:758015/140‑1 (MQ=255)
                                                                                                               gATTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCtt                        >  2:1108798/1‑140 (MQ=255)
                                                                                                                aTTTTCAGCGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCttt                       <  2:308176/140‑1 (MQ=255)
                                                                                                                aTTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTCTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCttt                       <  1:244707/140‑1 (MQ=255)
                                                                                                                  tttCAGCGGTACGTTGCTCGTCAGCTTCATGTAACTAGCTCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCa                     <  2:485697/140‑1 (MQ=255)
                                                                                                                        cGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATg               >  2:298417/1‑140 (MQ=255)
                                                                                                                        cGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATg               >  1:1228661/1‑140 (MQ=255)
                                                                                                                        cGGTACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCAATGGAAAAGTTGAGGAAAAAACCGTGATCCCGGCGGCTTTCAGTGATg               >  2:723227/1‑140 (MQ=255)
                                                                                                                           tACGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTaa            <  2:1597500/140‑1 (MQ=255)
                                                                                                                             cGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAg                    <  1:1887167/130‑1 (MQ=255)
                                                                                                                             cGTTGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAg                    >  2:1887167/1‑130 (MQ=255)
                                                                                                                               ttGCTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCa        <  1:412909/140‑1 (MQ=255)
                                                                                                                                  cTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAGCCGTTATCCCGCCGTCTTTCAGTGATGTAAGCCATcc     >  1:625703/1‑140 (MQ=255)
                                                                                                                                  cTCGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATcc     >  2:1577174/1‑140 (MQ=255)
                                                                                                                                    cGTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCCTTCAGTGATGTAAGCCATCCCt   >  1:1239854/1‑140 (MQ=255)
                                                                                                                                     gTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCGGACTGCGTCTGATCGCGATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGGCATCCCtt  >  2:1886748/1‑140 (MQ=255)
                                                                                                                                     gTCAGCTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAAGCGTTATCGCGGCGGCTTTCAGTGATGTAAGTCATCCCtt  >  2:892341/1‑140 (MQ=255)
                                                                                                                                         |                                                                                                                                       
CCGGGCAACAGCACAGCCATTTGGAGATCAGGACTTTCTGCTGATTTCCCCCGGAGAGTTCGGTGATATTCTGGTTTACCGAATGACATTTCAGCGCCAGCAGTTCGCGTTGATTTTCAGCGGTACGTTGCTCGTCAACTTCATGAAACAAGCCCATCGCGCCTTTATAGCCGCCGTCTTTCAGACTGCGGCTGATCGCCATGTTCTGAGCGATGGAAAAGTTGGGGAAAAAACCGTTATCCCGGCGGCTTTCAGTGATGTAAGCCATCCCTT  >  NC_000913/4309696‑4309968

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