Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,120,969 G→A L196L (CTG→TTG)  wcaJ ← UDP‑glucose:undecaprenyl‑phosphate glucose‑1‑phosphate transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,120,9690GA96.4% 180.6 / NA 56L196L (CTG→TTG) wcaJUDP‑glucose:undecaprenyl‑phosphate glucose‑1‑phosphate transferase
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (26/28);  minor base C (2/0);  total (28/28)
Fisher's exact test for biased strand distribution p-value = 4.91e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

GTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAGCTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGCCATGATTACGCAGCCAGCCCGCCCCAATGCGAATACACGAACGGCAAACCACCAGTCCGATGCTGGTCAGCGCATACC  >  NC_000913/2120732‑2121211
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   
gTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGtt                                                                                                                                                                                                                                        >  1:400422/1‑250 (MQ=255)
        cggCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCATCACCGGTACGCCGTTCATCGCTTCGAGGCGTGAATTGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAGCAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGCTTACCCGCCCa                                                                                                                                                                                                                               <  2:343360/250‑1 (MQ=255)
         ggCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGCATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGccc                                                                                                                                                                                                                                >  1:528650/1‑249 (MQ=255)
          gCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACACTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGt                                                                                                                                                                                                                             >  1:543774/1‑251 (MQ=255)
             gTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTGCGCGTCCTCGACCAACTGTTGCAGGTTCCCCGCCCAGTCGt                                                                                                                                                                                                                          <  2:403780/251‑1 (MQ=255)
                aaCCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGGTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTa                                                                                                                                                                                                                        <  2:473062/250‑1 (MQ=255)
                 aCCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGCATGGAGAATGTTGAAGGTAAAGCCGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGtt                                                                                                                                                                                                                         >  2:277644/1‑248 (MQ=255)
                   cccGGAAAGCGGCTTGTCGTAAAGCTGCACCACCGGCACGCCGTTCATCTCTTCGGCGCGTGAATGGAGAATGTTGAAGGTAACGCCGTCGGGGATCAGCAGCACCGAACAGGAGGTGTCCGCCAGTTGATGGCCCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCACCTGTTGCAGGTTACCCGCCCAGTCGTTAGaac                                                                                                                                                                                                                    >  2:177465/1‑250 (MQ=255)
                      ggAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGAGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAAcgc                                                                                                                                                                                                                 >  1:348078/1‑251 (MQ=255)
                           gcggcgTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGccc                                                                                                                                                                                                            <  1:483491/251‑1 (MQ=255)
                            cggcgTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCg                                                                                                                                                                                                           <  2:528650/251‑1 (MQ=255)
                                  gtCGGAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCTTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAAGCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTc                                                                                                                                                                                                      <  1:516053/250‑1 (MQ=255)
                                       aaaGCGGCACCACCGTTACGCCGTTCAGCTCTTCGAGGCGTGAATTGAGAATGTTGAAGGGAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTGTTCACTCGCGCGCCGGCGCTCATTTGCATCGCGATATAGACGTTAGGAATCTTGCCCGCGTTCGCGTGCTCGACCAACTGTTGCAGTGTTCCCGCCCAGGCGTTAGAAACGCCTCCCCGTTTCggggc                                                                                                                                                                                                 <  2:515194/250‑3 (MQ=255)
                                        aaGCGGCACCACCGGTAAGACGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGCCGTCGGGGCTCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCCGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGt                                                                                                                                                                                               >  2:365093/1‑251 (MQ=255)
                                           cGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGt                                                                                                                                                                                            <  1:296868/251‑1 (MQ=255)
                                                                        cGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCa                                                                                                                                                                <  2:140418/250‑1 (MQ=255)
                                                                                  aTGGAGAATGTTGAAGGTAAAGCCGGCGGGGATCAGCAGCCCCGAACAGGTGGGGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGCCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCGCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTAc                                                                                                                                                     <  2:310088/251‑1 (MQ=255)
                                                                                  aTGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTAc                                                                                                                                                     <  2:340664/251‑1 (MQ=255)
                                                                                     gagaATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGa                                                                                                                                                  <  1:539914/251‑1 (MQ=255)
                                                                                       gaATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGa                                                                                                                                                  >  1:122681/1‑249 (MQ=255)
                                                                                                     aggACGTCGGTGATCAGCAGCACCGAACAGTGGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCTCGCCGTCGCACAGTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAATCCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCAt                                                                                                                                  <  2:301037/249‑1 (MQ=255)
                                                                                                            cGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTCGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCcg                                                                                                                           <  1:289725/251‑1 (MQ=255)
                                                                                                               ggATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGGTGGTTACGGAAGCTCTCCATCAGCATTTGCCcggcg                                                                                                                        >  2:49095/1‑251 (MQ=255)
                                                                                                                   cagcagCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCt                                                                                                                     <  1:295198/250‑1 (MQ=255)
                                                                                                                        gcaCCGAACCGGTGGTGACCGCCCGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGACCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATcc                                                                                                               >  2:76726/1‑251 (MQ=255)
                                                                                                                        gcaCAGAACAGGGGGTGTCCGCGAGTTGATGGACCAGTTTTTGCACGCGCGCGCCGGCGCACATTTGCATCGCGATATAGACGTTAGGAAGCATGCCCGCATTCGCGTCATCGGCCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGACGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCGCTCCATCAGCATTTGTCCGGCGGCTAAATcc                                                                                                               <  1:177465/251‑1 (MQ=255)
                                                                                                                           ccGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGGTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATccccc                                                                                                            <  1:49095/251‑1 (MQ=255)
                                                                                                                                   gtggtgTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACcgc                                                                                                    >  2:305940/1‑251 (MQ=255)
                                                                                                                                          ccgccatgtgagtgACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTGGCCCGCTTTCGCGTCCTCTACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACTCCGCCCGGTTTCGGGTCGTGGTAATCTCCCACCACTGCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATg                                                                                             <  1:337061/238‑1 (MQ=255)
                                                                                                                                          ccGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATg                                                                                             <  1:109029/251‑1 (MQ=255)
                                                                                                                                               aGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCtt                                                                                        >  1:201732/1‑251 (MQ=255)
                                                                                                                                                  tGATGGACCCGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGCTATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACAGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAACTCCCCCCCCACCACGACCATGCGCTTGtt                                                                                     >  1:192614/1‑251 (MQ=255)
                                                                                                                                                   gATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCCCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCTCTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTta                                                                                    <  1:178202/251‑1 (MQ=255)
                                                                                                                                                      ggACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAg                                                                                 <  2:119396/251‑1 (MQ=255)
                                                                                                                                                        aCCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGcc                                                                               <  1:424264/251‑1 (MQ=255)
                                                                                                                                                        aCCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGcc                                                                               >  1:146060/1‑251 (MQ=255)
                                                                                                                                                          cAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGGCGTGGTAAACGCCCACCACTTCAAACCCTAAACACGGCTGGTTACGGAAGCTCTCAATCAGCATTTACCAGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGCCAt                                                                             >  1:319038/1‑251 (MQ=255)
                                                                                                                                                               ttttCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGCCATGAtt                                                                         <  1:305940/250‑1 (MQ=255)
                                                                                                                                                                 ttCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGCCATGATTACg                                                                      >  1:183971/1‑251 (MQ=255)
                                                                                                                                                                        gcgcgcCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCa                                                                                                                                         >  2:430853/1‑177 (MQ=255)
                                                                                                                                                                        gcgcgcCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCa                                                                                                                                         <  1:430853/177‑1 (MQ=255)
                                                                                                                                                                        gcgcgcCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGATTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCCAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGCCATGATTACGcagccac                                                               >  1:359375/1‑250 (MQ=255)
                                                                                                                                                                          gcgcCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGCCATGATTACGCagccagcc                                                             <  2:400422/251‑1 (MQ=255)
                                                                                                                                                                            gcCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTATCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGCCATGATTACGCAGCCAgcccg                                                           >  1:262732/1‑251 (MQ=255)
                                                                                                                                                                                                 gATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGCCATGATTACGCAGCCAGCCCGCCCCAATGCGAATACACGAAc                                      >  2:272005/1‑251 (MQ=255)
                                                                                                                                                                                                  atatAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATgc                                                                                            >  1:501356/1‑196 (MQ=255)
                                                                                                                                                                                                  atatAGACGTTATGAAGCTTGCCCGCTTTCGCGGCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATgc                                                                                            <  2:501356/196‑1 (MQ=255)
                                                                                                                                                                                                     taGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGCCATGATTACGCAGCCAGCCCGCCCCAATGCGAATACACGAACGGc                                   >  1:93295/1‑250 (MQ=255)
                                                                                                                                                                                                           ttatgcctCTTGCCCGATTTCGCGCACTCCACCACCTGTTGCAGGTTACCCGCCCAGTCGTTAGACACGCCGCCCGGTTTCGGGACGTGGTAAACGCCCACCCCTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGAATTTGACCGGCGGCTAACTACCCCGCCACCGCGACCATGCGCTTGTTATAGCCATCATTCCGCAGCCAGCCCGCCCCAATGCGAATACACAAAAGGCAAaccac                            >  2:261722/8‑251 (MQ=255)
                                                                                                                                                                                                           ttATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTAACCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCAACACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCCTGCGCTTGTTATAGCCAGGAGTACGCAGCCAGCCCGCCCCAATGCGAATACACGAACGGCAAaacac                            >  2:399646/1‑251 (MQ=255)
                                                                                                                                                                                                                   cTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGCCCATGCGCTTGTTATAGCCATGATTACGCAGCCAGCCCGCCCCAATGCGAATACACGAACGGAAAACCACCAGTCCg                     >  1:462645/1‑250 (MQ=255)
                                                                                                                                                                                                                      gCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGCCATGATTACGCAGCCAGCCCGCCCCAATGCGAATACACGAACGGCAAACCACCAGTCCGATGc                 >  1:89440/1‑251 (MQ=255)
                                                                                                                                                                                                                            gttCGCGTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTTTTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGCCATGATTACGCAGCCAGCCCGCCCCATGGCTAATACACGAACGGCAAACGACCAGTCCTATGCTGGTCa           <  2:359375/250‑1 (MQ=255)
                                                                                                                                                                                                                                gcgTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGCCATGATTACGCAGCCAGCCCGCCCCAATGCGAATACACGAACGGCAAACCACCAGTCCGATGCTGGTCAgcgc       >  2:103426/1‑251 (MQ=255)
                                                                                                                                                                                                                                 cgTCCTCGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACAACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGCCATGATTACGCAGCCAGCCCGCCCCAATGCGAATACACGAACGGCAAACAACCAGTCCGATGCTGGTCAgcgc       >  1:472641/1‑250 (MQ=255)
                                                                                                                                                                                                                                     cccGACCAACTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTGCGGGTCGTGGTAACCGCCCACCACTTCAAACCCTAACCCCGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGCCATGATTACGCAGCCAGCCCGCCCCAATGCGAATACACGAACGGCAAACCACCAGTCCGATGCTGGTCAGCGCATAcc  <  2:183971/249‑1 (MQ=255)
                                                                                                                                                                                                                                             |                                                                                                                                                                                                                                                   
GTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAGCTGTTGCAGGTTACCCGCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGCCATGATTACGCAGCCAGCCCGCCCCAATGCGAATACACGAACGGCAAACCACCAGTCCGATGCTGGTCAGCGCATACC  >  NC_000913/2120732‑2121211

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