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,9690GA87.3% 140.7 / 12.3 55L196L (CTG→TTG) wcaJUDP‑glucose:undecaprenyl‑phosphate glucose‑1‑phosphate transferase
Reads supporting (aligned to +/- strand):  ref base G (4/3);  new base A (26/22);  total (30/25)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.51e-01

GCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAGCTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGCCATGATTACGCAGCCAGCCCGCCCCAATGCGAATACACGAACGGCAAACCACCAGTCCGATGCTGGTCAGCGCATACCACG  >  NC_000913/2120728‑2121214
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      
gCACGTTTGAGCAGGCGGTTAACCACGGAAAGCTGCGTGTCGTAAAGCGTCACCACCGGTACGCCGTTCATCTCTTCGATGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGCTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCAGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGGTTTCGCGTCCTCGACCAACTGTTGCAg                                                                                                                                                                                                                                               <  2:14500/251‑1 (MQ=255)
      ttGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTAc                                                                                                                                                                                                                                          <  1:372831/250‑1 (MQ=255)
         aGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑Gc                                                                                                                                                                                                                                     >  1:239518/1‑251 (MQ=255)
               cGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTAACC‑GCCCAGTc                                                                                                                                                                                                                               >  2:371492/1‑251 (MQ=255)
                ggTTACCCCCGGAAAGCGGCGTGTCGTAAAGCGGCGCCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCGTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCg                                                                                                                                                                                                                              <  1:390174/251‑1 (MQ=255)
                      ccccGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGTCCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGaa                                                                                                                                                                                                                        <  1:383493/251‑1 (MQ=255)
                        ccGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAGCTGTTGCAGGTTACCC‑GCCCAGTCGTTAGnacc                                                                                                                                                                                                                      >  2:444090/1‑249 (MQ=255)
                           gAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAAcgcc                                                                                                                                                                                                                   >  1:241629/1‑251 (MQ=255)
                                         gTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCTCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGt                                                                                                                                                                                                     <  1:286137/251‑1 (MQ=255)
                                         gTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGt                                                                                                                                                                                                     >  1:102253/1‑251 (MQ=255)
                                           aaaGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGAACCGAACAGGTGGTGTCCGGCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCTTCCTCGACCAGCTGTTGCAGGTTTCCC‑TCCCAGTCGGTAGAAACGCCGCCCGTTTTCGGGTCg                                                                                                                                                                                                   >  1:342407/1‑251 (MQ=255)
                                                 gCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTa                                                                                                                                                                                              >  1:187627/1‑250 (MQ=255)
                                                       ccGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAGCTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGccc                                                                                                                                                                                       >  2:100755/1‑251 (MQ=255)
                                                         ggTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGCATGTTGACGGTCAACACCACGGGGATCAGCACCCCCGAACAGGTGGGGTCCGCCAGTTGTTAGCCCCGTTTTTTCACTCGCGGGCCGACGCACATCTGCCTCGCGATATAGCCGTTATGAATCTTGCACGCTTTCGCGTCCTCGACGAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCCCCCGGCTTCGGGCCGTGGTAAACGccccc                                                                                                                                                                                     >  2:65874/1‑249 (MQ=255)
                                                              gCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACtt                                                                                                                                                                                >  2:394762/1‑251 (MQ=255)
                                                                          ttCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAAcc                                                                                                                                                                    >  1:483150/1‑251 (MQ=255)
                                                                              aggcgacaaTGAAGAATGGTGAAGGCGAAGACGTCGGGGATCAGCCGCACCGAACCGGAAGGGTCCGCCAGTTGATCGACCAGCTTTTTAACTCGCGCGACGTCGCACATTTGCATCGCGATATAGCAGCAATGAATCTTGCACGCCCTCGCGTCATCGCCCAACTGTTGCCGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCTGGTTTAGGGTCGTGGTCAACGCCAACTACTTCACACCCTAACCNCgg                                                                                                                                                                >  2:441259/8‑251 (MQ=255)
                                                                                gCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAGCTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCt                                                                                                                                                              <  1:471344/251‑1 (MQ=255)
                                                                                  gTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAGCTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTgg                                                                                                                                                            >  2:244947/1‑251 (MQ=255)
                                                                                   tGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGt                                                                                                                                                           >  1:370138/1‑251 (MQ=255)
                                                                                         gagaATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGa                                                                                                                                                     <  1:248300/251‑1 (MQ=255)
                                                                                           gaATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAg                                                                                                                                                   <  1:290560/251‑1 (MQ=255)
                                                                                            aaTGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGc                                                                                                                                                  >  2:383230/1‑251 (MQ=255)
                                                                                                 tGAAGGTAAATACGACGGGGATAAGCAGCACCGAACAGGTGGTGTACGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGGCGCACATTTGCATCGAGATATAAACGTTATGAATCATGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCAGGTTTCGGGTCGTGGTAAACACCCAACACTTCAAACCCTAACCACGGCTGGTACCGGAAGCTCTcc                                                                                                                                             >  2:74541/1‑251 (MQ=255)
                                                                                                     ggTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATc                                                                                                                                          >  1:84751/1‑250 (MQ=255)
                                                                                                           gACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAGCTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCAttt                                                                                                                                   <  1:444090/251‑1 (MQ=255)
                                                                                                                  gggATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGCTACGGAAGCTCTCCATCCGCATTTACCcggc                                                                                                                            >  1:256603/1‑251 (MQ=255)
                                                                                                                       cagcagCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCt                                                                                                                        >  1:47094/1‑250 (MQ=255)
                                                                                                                           agcaCCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGTCGGCTAAATc                                                                                                                   >  1:397118/1‑251 (MQ=255)
                                                                                                                              cccGAACAGGTGGTGTCCGCCAGTTGATGTACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATcccc                                                                                                                <  1:388088/250‑1 (MQ=255)
                                                                                                                               ccGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGc                                                                                                                                 >  1:377204/1‑233 (MQ=255)
                                                                                                                               ccGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGc                                                                                                                                 <  2:377204/233‑1 (MQ=255)
                                                                                                                                 gAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGc                                                                                                             <  2:241629/251‑1 (MQ=255)
                                                                                                                                  aaCAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGGCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCCGGCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGcc                                                                                                            <  1:109331/251‑1 (MQ=255)
                                                                                                                                     aGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCAcc                                                                                                         >  2:22085/1‑251 (MQ=255)
                                                                                                                                      ggtggtGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCAccc                                                                                                        >  1:136066/1‑250 (MQ=255)
                                                                                                                                                  cAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCt                                                                                            <  2:102253/251‑1 (MQ=255)
                                                                                                                                                    gTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTg                                                                                          <  1:111937/251‑1 (MQ=255)
                                                                                                                                                             ccAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAAc                                                                                                                                                                           >  2:113391/1‑161 (MQ=255)
                                                                                                                                                             ccAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACGTCAAAc                                                                                                                                                                           <  1:113391/161‑1 (MQ=255)
                                                                                                                                                                       cACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAGCTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGCCATGATTACGc                                                                        <  1:244947/250‑1 (MQ=255)
                                                                                                                                                                            gcgcgcCGGCGTACATTTGAATCGCGATATAGACGTTATGAATCTTGCCCTCTTTCGCGTCCCCGACCAACTGTTGCAGGTTCCCC‑GCCCAGTCGTTAGAAAAGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCTATTCAAACCCTACCCACGGCTGGTTACTGAAGCTCTCCATCAGTATGTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGCCATGATTACGcagccag                                                                  <  2:426498/251‑1 (MQ=255)
                                                                                                                                                                               cgcCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCTTCCTCTACCAACTGTTGCAGGTTACCC‑GCGCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTTCCCGGCGGCGAAGTCCCCCGTCACCGCGCCCATGCGCTTGTTATAGCCATGATTACGCAGCCAgccc                                                               <  1:65874/251‑1 (MQ=255)
                                                                                                                                                                                 ccGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGCCATGATTACGCAGCCAgcccgc                                                             >  1:412111/1‑251 (MQ=255)
                                                                                                                                                                                                aTCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGCCATGATTACGCAGCCAGCCCGCCCCACTGCGAATac                                               >  1:352894/1‑250 (MQ=255)
                                                                                                                                                                                                  cgcgATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGCCATGATTACGCAGCCAGCCCGCCCCAATGCGAATacac                                             <  1:305820/250‑1 (MQ=255)
                                                                                                                                                                                                    tgATATAGACGTTATTAATCTTGCCCGCTTTCGCGTCCTCTACCAACTGTTGC‑CGTTACCCCGCCCAGTCGTTAGAAACTCCGCCCGGTTTCGGGTCGGGTTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGCCATGATTACGCAGCCAGCCCGCCCCAATGCGAATACACg                                            <  2:370138/248‑1 (MQ=255)
                                                                                                                                                                                                     gATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGCCATGATTACGCAGCCAGCCCGCCCCAATGCGAATACACGAAc                                         <  1:371492/251‑1 (MQ=255)
                                                                                                                                                                                                        ataGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGCCATGATTACGCAGCCAGCCCGCCCCAATGCGAATACACGAACgg                                       <  1:383230/250‑1 (MQ=255)
                                                                                                                                                                                                             cGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGCCATGATTACGCAGCCAGCCCGCCCCAATGCGAATACACGAACGGCAAacc                                 >  2:89602/1‑251 (MQ=255)
                                                                                                                                                                                                                    aaTCTTGCCCGCTTTCGCGTCCGCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTATAACCGCCGCCCGGTTTCGGGTCGTGGTAAACGACCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGCCATGATTACGCAGCCAGCCCGCCCCAATGCGAATACACGAACGGCAAACCACCAGTc                          <  1:25910/251‑1 (MQ=255)
                                                                                                                                                                                                                       cTTGCCCGCTTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGCCATGATTACGCAGCCAGCCCGCCCCAATGCGAATACCCGAACGGCAAACCACCAGTCCg                        >  2:81638/1‑250 (MQ=255)
                                                                                                                                                                                                                          gCCCGATTTCGCGTCCTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGCCATGATTACGCAGCCAGCCCGCCCCAATGCGAATACACGAACGGCAAACCACCAGTCCGNTGc                    >  2:436311/1‑251 (MQ=255)
                                                                                                                                                                                                                                        ccTCGACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCGAACCACGGCTGGTTACGGAAGCTCTCCATCATCATTTGCCCGGCGGCTAAATCCCCCGCCGCCGCGACCATGCGCTTGTTATAGCCATGATTACGCAGCCAGCCCGCCCCAATGCGAATACACGAACGGCAAACCACCAGTCCGATGCTGGTCAGCGCATa       <  1:481266/250‑1 (MQ=255)
                                                                                                                                                                                                                                            gACCAACTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCCTCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGCCATGATTACGCAGCCAGCCCGCCCCAATGCGAATACACGAACTGCAAACCACCAGTCCCATGCTGGTCAGCGCATACCACg  >  1:466024/1‑251 (MQ=255)
                                                                                                                                                                                                                                                 |                                                                                                                                                                                                                                                      
GCACGTTTGAGCAGGCGGTTAACCCCGGAAAGCGGCGTGTCGTAAAGCGGCACCACCGGTACGCCGTTCATCTCTTCGAGGCGTGAATGGAGAATGTTGAAGGTAAAGACGTCGGGGATCAGCAGCACCGAACAGGTGGTGTCCGCCAGTTGATGGACCAGTTTTTTCACTCGCGCGCCGTCGCACATTTGCATCGCGATATAGACGTTATGAATCTTGCCCGCTTTCGCGTCCTCGACCAGCTGTTGCAGGTTACCC‑GCCCAGTCGTTAGAAACGCCGCCCGGTTTCGGGTCGTGGTAAACGCCCACCACTTCAAACCCTAACCACGGCTGGTTACGGAAGCTCTCCATCAGCATTTGCCCGGCGGCTAAATCCCCCGCCACCGCGACCATGCGCTTGTTATAGCCATGATTACGCAGCCAGCCCGCCCCAATGCGAATACACGAACGGCAAACCACCAGTCCGATGCTGGTCAGCGCATACCACG  >  NC_000913/2120728‑2121214

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: ‑