Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,524,463 Δ1 bp 12.5% intergenic (+95/+18) yncE → / ← ansP PQQ‑like domain‑containing protein YncE/L‑asparagine transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,524,4630A.12.5% 247.0 / 28.5 80intergenic (+95/+18)yncE/ansPPQQ‑like domain‑containing protein YncE/L‑asparagine transporter
Reads supporting (aligned to +/- strand):  ref base A (36/34);  new base . (3/7);  total (39/41)
Fisher's exact test for biased strand distribution p-value = 3.13e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.96e-01

AAACAGCAGGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCAAAATACCGA  >  NC_000913/1524318‑1524611
                                                                                                                                                  |                                                                                                                                                    
aaaCAGCAGGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCg                                                                                                                                                  >  2:266171/1‑150 (MQ=255)
  aCAGCAGGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGTGCCggg                                                                                                                                                <  2:209236/150‑1 (MQ=255)
    agcagGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGTGCCGGGCa                                                                                                                                              <  2:1974965/150‑1 (MQ=255)
    agcagGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCa                                                                                                                                              >  1:324489/1‑150 (MQ=255)
    agcagGAAGCTACCCAGCCAGACCATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGTGCCGGGCa                                                                                                                                              <  2:1392257/150‑1 (MQ=255)
         gAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGTGCCGGGCAAattc                                                                                                                                         <  1:1665033/150‑2 (MQ=255)
          aaGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCTACGTGCCGGGCAAattca                                                                                                                                        <  2:950690/150‑3 (MQ=255)
             cTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCAGGCACGTAGCCGGGCAATTTTGCaa                                                                                                                                     >  1:857571/1‑150 (MQ=255)
               aCCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAAtt                                                                                                                                   <  1:299553/150‑1 (MQ=255)
                cccAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTa                                                                                                                                  >  2:1366029/1‑150 (MQ=255)
                cccAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTa                                                                                                                                  <  2:202645/150‑1 (MQ=255)
                   agccagACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACt                                                                                                                               <  1:1083527/150‑1 (MQ=255)
                          cGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGtt                                                                                                                        >  1:1082386/1‑150 (MQ=255)
                               tgATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAgg                                                                                                                   >  2:769128/1‑150 (MQ=255)
                                gATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGc                                                                                                                  >  2:118552/1‑150 (MQ=255)
                                  ttCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCtt                                                                                                                >  2:447009/1‑150 (MQ=255)
                                    cGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTaa                                                                                                              <  2:1991989/150‑1 (MQ=255)
                                     gTATTGCGCTCTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTaaa                                                                                                             >  2:1115046/1‑150 (MQ=255)
                                       aTTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAaca                                                                                                           >  2:1054902/1‑150 (MQ=255)
                                           cgcTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAAt                                                                                                       >  1:308707/1‑150 (MQ=255)
                                             cTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAAttt                                                                                                     <  1:769128/150‑1 (MQ=255)
                                               gTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTcc                                                                                                   <  2:362145/150‑1 (MQ=255)
                                               gTAATAGATATGGCGAGGGTGAAAACGGGGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGGAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTcc                                                                                                   <  2:936547/150‑1 (MQ=255)
                                                 aaTAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTg                                                                                                 >  2:234904/1‑150 (MQ=255)
                                                  aTAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGt                                                                                                <  2:1613906/150‑1 (MQ=255)
                                                     gATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGtttt                                                                                             <  1:1198725/150‑1 (MQ=255)
                                                        atGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTttct                                                                                          <  2:449219/150‑1 (MQ=255)
                                                          ggCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCATTTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTttcttc                                                                                        <  1:1115046/150‑1 (MQ=255)
                                                           gCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCa                                                                                       <  2:128832/150‑1 (MQ=255)
                                                           gCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCa                                                                                       <  2:1082386/150‑1 (MQ=255)
                                                           gCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGACGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCa                                                                                       <  1:964494/150‑1 (MQ=255)
                                                             gAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCtt                                                                                                                >  1:1147566/1‑123 (MQ=255)
                                                             gAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCtt                                                                                                                <  2:1147566/123‑1 (MQ=255)
                                                                  tGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCt                                                                                >  2:71114/1‑150 (MQ=255)
                                                                     aaaCGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTcga                                                                             <  2:308707/150‑1 (MQ=255)
                                                                     aaaCGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTcga                                                                             <  1:1751338/150‑1 (MQ=255)
                                                                        cGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTcgacga                                                                          <  1:374089/150‑1 (MQ=255)
                                                                        cGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTcgacga                                                                          >  1:175503/1‑150 (MQ=255)
                                                                           tGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTg                                                                       >  2:532016/1‑150 (MQ=255)
                                                                              gaCTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGgcg                                                                    <  2:30473/150‑1 (MQ=255)
                                                                               aCTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGgcgc                                                                   <  1:118552/150‑1 (MQ=255)
                                                                                cTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCa                                                                  >  1:1800042/1‑150 (MQ=255)
                                                                                    cctgcctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGAGTGGCGCAGTGc                                                              >  2:665779/1‑150 (MQ=255)
                                                                                       gcctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCtgt                                                           <  1:1446652/150‑1 (MQ=255)
                                                                                        cctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGGCGGGCAATTTTGCATTTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCtgtg                                                          <  2:1939952/150‑1 (MQ=255)
                                                                                        cctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCtgtg                                                          >  2:1474080/1‑150 (MQ=255)
                                                                                        cctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCtgtg                                                          <  1:118662/150‑1 (MQ=255)
                                                                                         ctTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGa                                                         <  1:266171/150‑1 (MQ=255)
                                                                                          ttttATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGaa                                                        <  2:936123/150‑1 (MQ=255)
                                                                                             tATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAAttt                                                     >  2:612083/1‑150 (MQ=255)
                                                                                                cACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAg                                                  >  2:1418263/1‑150 (MQ=255)
                                                                                                cACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAg                                                  >  2:857286/1‑150 (MQ=255)
                                                                                                 aCGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGc                                                 <  1:1214509/150‑1 (MQ=255)
                                                                                                     cGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAAcg                                             >  2:661066/1‑150 (MQ=255)
                                                                                                        aTCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGt                                          >  2:363938/1‑150 (MQ=255)
                                                                                                                 tttCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGcaca                                 <  2:979046/150‑1 (MQ=255)
                                                                                                                   tCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACAcc                               <  1:1474080/150‑1 (MQ=255)
                                                                                                                    cGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCa                              >  1:1818252/1‑150 (MQ=255)
                                                                                                                       cAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAAc                           >  1:580976/1‑150 (MQ=255)
                                                                                                                       cAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAAc                           >  1:814923/1‑150 (MQ=255)
                                                                                                                        aGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAAcc                          >  1:1635736/1‑150 (MQ=255)
                                                                                                                            aTCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTTAATTTCATCAACGCGTTTGCGCACACCAAACCAGca                      <  2:324489/150‑2 (MQ=255)
                                                                                                                            aTCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCa                         >  2:1255962/1‑147 (MQ=255)
                                                                                                                            aTCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCa                         <  1:1255962/147‑1 (MQ=255)
                                                                                                                             tCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCt                     <  1:647305/150‑1 (MQ=255)
                                                                                                                              cTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCta                    <  1:1366029/150‑1 (MQ=255)
                                                                                                                                 aaTTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATaa                 >  1:211314/1‑150 (MQ=255)
                                                                                                                                  aattGCCCGGCTACGTGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAAc                >  2:1737278/3‑150 (MQ=255)
                                                                                                                                   tttGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAAcc               >  1:1104063/1‑150 (MQ=255)
                                                                                                                                    ttGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTTGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTCTGCGCACACCAAACCAGCCTATAACCa              <  1:717339/150‑1 (MQ=255)
                                                                                                                                     tGCCCGGCTACGTGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAg                        >  2:1821277/1‑139 (MQ=255)
                                                                                                                                     tGCCCGGCTACGTGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAg                        <  1:1821277/139‑1 (MQ=255)
                                                                                                                                     tGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAg             >  2:1074538/1‑150 (MQ=255)
                                                                                                                                      gCCCGGCTACGTGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGc            <  1:1737278/150‑1 (MQ=255)
                                                                                                                                      gCCCGGCAACGTGCCCGGCAATTTTGCAATTAACTCGCCGATTCAGGCTTAAACACAAATTCCTCTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGc            >  2:716892/1‑150 (MQ=255)
                                                                                                                                      gCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGc            >  1:486293/1‑150 (MQ=255)
                                                                                                                                          ggCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCaaaa        >  1:1694621/1‑150 (MQ=255)
                                                                                                                                          ggCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCaaaa        >  1:124902/1‑150 (MQ=255)
                                                                                                                                                 tAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCAAAATACCGa  >  2:601199/1‑150 (MQ=255)
                                                                                                                                                 tAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCAAAATACCGa  <  1:612083/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
AAACAGCAGGAAGCTACCCAGCCAGACGATGTGATTCGTATTGCGCTGTAATAGATATGGCGAGGGTGAAAACGGTGAGACTTGCCTGCCTTTTATCACGGCGTATCCGGCACTTTCGTCAGTAATCTGAATTTGCCCGGCACGTAGCCGGGCAATTTTGCAATTAACTCGCCGTTTCAGGCTTAAACACAATTTCCTGTTTTTCTTCATCTTCCTCGACGACTGGCGCAGTGCTGTGAATTTCAGCAACGCGTTTGCGCACACCAAACCAGCCTATAACCAGCAAAATACCGA  >  NC_000913/1524318‑1524611

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