Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 4,026,356 A→G N457D (AAT→GAT)  ubiD → 3‑octaprenyl‑4‑hydroxybenzoate decarboxylase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,026,3560AG89.2% 190.8 / 10.5 74N457D (AAT→GAT) ubiD3‑octaprenyl‑4‑hydroxybenzoate decarboxylase
Reads supporting (aligned to +/- strand):  ref base A (2/6);  new base G (31/35);  total (33/41)
Fisher's exact test for biased strand distribution p-value = 2.86e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.09e-01

CCGGACACGCGAAGCGCGTCATGATGGGCGTCTGGTCGTTCTTACGCCAGTTTATGTACACTAAATTTGTGATCGTTTGCGATGATGACGTTAACGCACGCGACTGGAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCA  >  NC_000913/4026111‑4026601
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
ccGGACACGCGAAGCGCGTCATGATGGGCGTCTGGTCGTTCTTACGCCAGTTTATGTACACTAAATTTGTGATCGTTTGCGATGATGACGTTAACGCACGCGACTGGAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTCTTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATaa                                                                                                                                                                                                                                                   >  1:331753/1‑250 (MQ=255)
    acacGCGAAGCGCGTCATGATGGGCGTCTGGTCGTTCTTACGCCAGTTTATGTACACTAAATTTGTGATCGTTTGCGATGATGACGTTAACGCACGCGACTGGAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAACGGc                                                                                                                                                                                                                                              <  2:177799/251‑1 (MQ=255)
          gaatcgcgTCATGATGGGCGTCTGGTCGTTCTTACGCCAGTTTATGTACACTAAATTTGTGATCGTTTGCGATGATGACGTTAACGCACGCGACTGGAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCggggg                                                                                                                                                                                                                                        <  2:8853/247‑1 (MQ=255)
               gcgTCATGATGGGCGTCTGGTCGTTCTTACGCCAGTTTATGTACACTAAATTTGTGATCGTTTGCGATGATGACGTTAACGCACGCGACTGGAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAAcc                                                                                                                                                                                                                                   <  1:469423/251‑1 (MQ=255)
                cgTCATGATGGGCGTCTGGTCGTTCTTACGCCAGTTTATGTACACTAAATTTGTGATCGTTTGCGATGATGACGTTAACGCACGCGACTGGAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAAccc                                                                                                                                                                                                                                  <  1:128082/251‑1 (MQ=255)
                                 gggcGTTCTGACGCCAGTTTATGTACACTAAATTTGTGATCGTTTGCGATGATGACGTTAACGCACGCGACTGGAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGAc                                                                                                                                                                                                                   <  2:80958/246‑1 (MQ=255)
                                 ggTCGTTCTTACGCCAGTTTATGTGCACTAAAGTTGTGATCGTTTGCGATGATGACGTTAACGCACGCGACTGGAACGAGGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGa                                                                                                                                                                                                                    <  2:341444/248‑1 (MQ=255)
                                 ggTCGTTCTTACGCCAGTTTATGTACACTAAATGTGTGATCGTTTGCGATGATGACGTTAACGCACGCGACTGGAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGAGTATCTGGAGTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGt                                                                                                                                                                                                                 <  2:263340/251‑1 (MQ=255)
                                    cGTTCTTACGCCAGTTTATGTACACTAAATTTGTGATCGTTTGCGATGATGACGTTAACGCACGCGACTGGAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTccc                                                                                                                                                                                                              >  2:502828/1‑251 (MQ=255)
                                    cGTTCTTACGCCAGTTTATGTACACTAAATTTGTGATCGTTTGCGATGATGACGTTAACGCACGCGACTGGAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTcc                                                                                                                                                                                                               <  1:116560/250‑1 (MQ=255)
                                      ttCTTACGCCAGTTTATGTACACTAAATTTGTGATCGTTTGCGATGATGACGTTAACGCACGCGACTGGAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCAt                                                                                                                                                                                                            >  1:348544/1‑251 (MQ=255)
                                        cTTACGCCAGTTTATGTACACTAAATTTGTGATCGTTTGCGATGCTAACGTTAACGCACGCGACTGGAAAGATGAGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATCCTGTTTTGGTCGACAATAAGCCTATTGATTATCTGGTTTTTGCCTCGCATGTCTCCAGGCTGGGTTCACAAATGGGGCTGGATGCCACGGCTAAATGGGCGGGGGACACCCAGCGTGAATGGGGACGTCCCATc                                                                                                                                                                                                           >  1:293212/1‑250 (MQ=255)
                                          tACGCCAGTTTATGTACACTAAATTTGTGATCGTTTGCGATGATGACGTTAACGCACGCGACTGGAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATTCCACGGATAAATGGCCGTGGGAAACCCAGCGTTAATGGGGACGTCCCATc                                                                                                                                                                                                           <  2:301882/248‑1 (MQ=255)
                                               caggttATGTACACTAAATTTGTGATCGGTTGCGATGATGACGTTAACGCACGCGACTGGAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTTATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGa                                                                                                                                                                                                   <  1:398567/247‑1 (MQ=255)
                                                aGTTTATGTACACTAAATTTGTGATCGTTTGCGATGATGACGTTAACGCACGCGACTGGAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGCTGCCACGGATAACTGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCaaaaa                                                                                                                                                                                                      >  1:25366/1‑247 (MQ=255)
                                                     aTGTACACTAAATTTGTGATCGTTTGCGATGATGACGTTAACGCACGCGACTGGAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAg                                                                                                                                                                                              <  1:231281/250‑1 (MQ=255)
                                                       gTACACTAAATTTGTGATCGTTTGCGATGATGACGTTAACGCACGCGACTGGAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGAt                                                                                                                                                                                            <  1:513000/250‑1 (MQ=255)
                                                          cacTAAATTTGTGATCGTTTGCGATGATGACGTTAACGCACGCGACTGGAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGAtgttg                                                                                                                                                                                        <  1:410529/251‑1 (MQ=255)
                                                           acTAAATTTGTGATCGTTTGCGATGATGACGTTAACGCACGCGACTGGAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGAtgttgt                                                                                                                                                                                       >  1:367370/1‑251 (MQ=255)
                                                             tAAATTTGTGATCGTTTGCGATGATGACGTTAACGCACGCGACTGGAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCACAAAAGATCCAGATGTTGTcg                                                                                                                                                                                     >  1:300621/1‑251 (MQ=255)
                                                              aaaTGTGTGATCGTTTGCGATGATGACGTTAACGCACGCGACTGGAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTcgc                                                                                                                                                                                    <  1:597762/251‑1 (MQ=255)
                                                               aaTTTGTGATCGTTTGCGATGATGACGTTAACGCACGCGACTGGAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTcgcg                                                                                                                                                                                   <  2:244042/251‑1 (MQ=255)
                                                               aaTTTGTGATCCATTGCGATGCTGACGTTAACGCCCGCGACTGGAACGATGTTCTTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTTGTAAATCCGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGTTAAATGGCCGGGGGAAACCCAACGAGAATGGTGACGTCCAAGCAAAAAAGATCCAGAAGTAGTcgcg                                                                                                                                                                                   >  2:500746/1‑251 (MQ=255)
                                                                         cGTTTGCGATGATGACGTTAACGCACGCGACTGGAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGAc                                                                                                                                                                          <  2:432868/250‑1 (MQ=255)
                                                                         cGTTTGCGATGATGACGTTAACGCACGCGACTGGAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACg                                                                                                                                                                         >  2:223261/1‑251 (MQ=255)
                                                                          gTTTGCGATGATGACGTTAACGCACGCGACTGGAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGc                                                                                                                                                                        >  2:527468/1‑251 (MQ=255)
                                                                           tttGCGATGATGACGTTAACGCACGCGACTGGAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAACAGATCCAGATGTTGTCGCGCATATTGACGcc                                                                                                                                                                       >  2:363478/1‑251 (MQ=255)
                                                                            ttGCGATGATGACGTTAACGCACGCGACTGGAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGcc                                                                                                                                                                       >  1:580138/1‑250 (MQ=255)
                                                                                       aCGTTAACGCACGCGACTGGAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGa                                                                                                                                                            >  1:292593/1‑250 (MQ=255)
                                                                                        cGTTAACGCACGCGACTGGAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAAc                                                                                                                                                          >  2:66546/1‑251 (MQ=255)
                                                                                         gTTAACGCCCGGGACTGGAACGATGGGCTTTGGGCGATTAACACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCACAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGACACCCAGCGTGAATGGGGACGTCCCATCAAAAAAAATCCAGATGTTGTCGCGCATATTGAAGCCCTCTGGGATGAACt                                                                                                                                                         >  2:380004/1‑251 (MQ=255)
                                                                                            aaCGCACGCGACTGGAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGc                                                                                                                                                      <  1:267459/251‑1 (MQ=255)
                                                                                                      aCTGGAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTaac                                                                                                                                            <  1:63109/251‑1 (MQ=255)
                                                                                                          gAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGTCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGTACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACg                                                                                                                                        <  1:383513/251‑1 (MQ=255)
                                                                                                           aaCGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACgg                                                                                                                                       >  1:172782/1‑251 (MQ=255)
                                                                                                           aaCGATGGGATGTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTGCTGGTAGAAAATACGCCTATTGATTATCTGGATGTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGGTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACgg                                                                                                                                       <  2:58178/251‑1 (MQ=255)
                                                                                                             cGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGCTTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAACTGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGt                                                                                                                                      >  1:394297/1‑250 (MQ=255)
                                                                                                               aTGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTaaa                                                                                                                                   <  1:527468/251‑1 (MQ=255)
                                                                                                                tgtgATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTaaa                                                                                                                                   >  1:95268/1‑250 (MQ=255)
                                                                                                                               tACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAAGTCGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGTCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCAGATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGTTAATAGCGCCTGAGGcgcg                                                                                                                    <  2:308571/250‑1 (MQ=255)
                                                                                                                                accaccCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCGCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCGAGATGTTGTCGCGCATATTGACGCCAGGTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCgtt                                                                                                                  <  1:269869/251‑1 (MQ=255)
                                                                                                                                accaccCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTAGTGATGATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGAGGTTGGCGCGCATATTGACGCCAGCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCgtt                                                                                                                  <  1:563368/251‑1 (MQ=255)
                                                                                                                                  caccCGTATGGTCCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTAGCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCgtttg                                                                                                                <  2:230267/251‑1 (MQ=255)
                                                                                                                                         aTGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGcc                                                                                                         <  1:179491/251‑1 (MQ=255)
                                                                                                                                           ggACCCGGCGCGGGATACTGTTCTGGTAGACAATACGCCTATTGATTATCTGGCTTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAACATGGGGCTGTATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAACGATCCAACTGTTGTCGCGCATATTGACGCCATCTGGGATGTACTGTCTATTTTTCACAACGGTAAAAGCGCCTGCTGCTCGTTTGTTTTTCCCt                                                                                                       >  2:26747/1‑251 (MQ=255)
                                                                                                                                            gACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTa                                                                                                      <  1:66546/251‑1 (MQ=255)
                                                                                                                                                 cgcgcGGGATACTGTTCTGGTAGAACCTACGCCTATTGATTATCTGGGTTTTCCCTCGACAGTCTCCAGGCTGGGATCAAAAATGGGGCTGGCTGCCACGGATAACTGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCACAAACGATCCAGATGTTGCCGCGAATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAACCGCCTGATGCGCGTTTGTTTTGCCCTATTTAt                                                                                                 >  2:58456/2‑251 (MQ=255)
                                                                                                                                                      aggATACTGTTCTGGTAGAAAATATGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACAGGCTATTTTTAACAACGGTCAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGAt                                                                                             >  2:443952/2‑250 (MQ=255)
                                                                                                                                                       ggATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATcc                                                                                           <  1:223947/251‑1 (MQ=255)
                                                                                                                                                        gATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCg                                                                                          >  1:326678/1‑251 (MQ=255)
                                                                                                                                                         aTACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCg                                                                                          >  1:237692/1‑250 (MQ=255)
                                                                                                                                                          tACTGTTCTGGTCGAAAATACGCCTATTGCTTATCTGGATTTTGCCTCGACTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGCTGCCACGGCTAACTGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCCTATTgccgcc                                                                                                                                                                       >  1:375477/1‑172 (MQ=255)
                                                                                                                                                              gTTCTGGTAGAATATATGCCTAGTGATTATCTGGATTTTGTCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGTCTGTAGGCCACGGAAAAATGGCCGGTGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGAGTCAGATGTTGTCGCGCAGATTGACGTCAGCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTATCCATATTTATCGATCCGACaga                                                                                     <  2:468679/250‑1 (MQ=255)
                                                                                                                                                                      aGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAgcgc                                                                             <  1:397562/250‑1 (MQ=255)
                                                                                                                                                                               gCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGAGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCt                                                                   >  2:386494/1‑251 (MQ=255)
                                                                                                                                                                                   aTTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTa                                                                 >  1:200450/1‑249 (MQ=255)
                                                                                                                                                                                   aTTGATTAGCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTTTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAg                                                               <  1:569774/251‑1 (MQ=255)
                                                                                                                                                                                   aTTGATTAGCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGGTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTa                                                                 <  1:394453/249‑1 (MQ=255)
                                                                                                                                                                                   aTTGAGTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATTTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAGAGCGCGTGATGCGCGGTTGTTTTGCCCTATTTAGCGATCCGACAGAGTAAGCGCATGACAACCTTaa                                                                <  2:310423/250‑1 (MQ=255)
                                                                                                                                                                                              ggATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTg                                                     >  1:209647/1‑250 (MQ=255)
                                                                                                                                                                                              ggATTTTGCATCGCCTGTCTCCGGGCTGGGTTCAAAAAGGGGGGTGGATGCTACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTaaa                                                        <  1:46183/247‑1 (MQ=255)
                                                                                                                                                                                                               tctcCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTGTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGAATGACAACCTTATGCTGTAAAGTGACAACGGTAGAAGCTAtt                                   >  2:513929/1‑250 (MQ=255)
                                                                                                                                                                                                                  ccGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACg                                >  2:614834/1‑251 (MQ=255)
                                                                                                                                                                                                                   cGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACgg                               <  1:580310/251‑1 (MQ=255)
                                                                                                                                                                                                                     ggCTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGGTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGAt                             <  1:599440/251‑1 (MQ=255)
                                                                                                                                                                                                                       cTGGGTTCAAAAATGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATAc                           >  1:451178/1‑251 (MQ=255)
                                                                                                                                                                                                                         gggTTCAAAAATGGGGCTGGATGGCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATAcc                          <  2:134330/250‑1 (MQ=255)
                                                                                                                                                                                                                              cAAAAATGGGGCTGGAAGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACACCTGTAAAAGCGCCTGATGCGCGTTTGTTGTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACTGATACCGtata                     <  1:323500/250‑1 (MQ=255)
                                                                                                                                                                                                                                  aaTGGGGCTGGATGCCACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGCCGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGACAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCgtg                >  2:227545/1‑251 (MQ=255)
                                                                                                                                                                                                                                               gCCACGGATAAATGGCCGGGGGAACCCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGGCGCCATCTGGGATGAACTGGCTATTTTTATCAACGGTAAAAGCGTCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGTCAACGTTAAGCTGTAAATTGATCTCGGTAGAAGCTATCACGGATACCGTATAGCGTGTCCGCATCGTGcc   <  1:58456/251‑1 (MQ=255)
                                                                                                                                                                                                                                                ccACGGATAAATGGCCGGGGGAAACTCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGGTTGTTTTGTCCTATTTAGCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTTACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCTTGCCa  <  2:33029/251‑1 (MQ=255)
                                                                                                                                                                                                                                                ccACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCa  <  1:99373/251‑1 (MQ=255)
                                                                                                                                                                                                                                                 cACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGcc   >  1:12830/1‑249 (MQ=255)
                                                                                                                                                                                                                                                 cACGGATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGCCGCGCATATTGCCGCCATCTGGGATGAACTGGCTATTTGTAAAAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATAGCTCCGACAGAGAAAGCGCATGACAACCTTACGCTGCAAGGTGACCTCGGAAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCATGcc   >  1:611317/1‑249 (MQ=255)
                                                                                                                                                                                                                                                     |                                                                                                                                                                                                                                                     
CCGGACACGCGAAGCGCGTCATGATGGGCGTCTGGTCGTTCTTACGCCAGTTTATGTACACTAAATTTGTGATCGTTTGCGATGATGACGTTAACGCACGCGACTGGAACGATGTGATTTGGGCGATTACCACCCGTATGGACCCGGCGCGGGATACTGTTCTGGTAGAAAATACGCCTATTGATTATCTGGATTTTGCCTCGCCTGTCTCCGGGCTGGGTTCAAAAATGGGGCTGGATGCCACGAATAAATGGCCGGGGGAAACCCAGCGTGAATGGGGACGTCCCATCAAAAAAGATCCAGATGTTGTCGCGCATATTGACGCCATCTGGGATGAACTGGCTATTTTTAACAACGGTAAAAGCGCCTGATGCGCGTTTGTTTTGCCCTATTTATCGATCCGACAGAGAAAGCGCATGACAACCTTAAGCTGTAAAGTGACCTCGGTAGAAGCTATCACGGATACCGTATATCGTGTCCGCATCGTGCCA  >  NC_000913/4026111‑4026601

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