Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,809,982 2 bp→AG coding (158‑159/480 nt) coaD → pantetheine‑phosphate adenylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,809,9820CA95.2% 122.4 / ‑6.0 42A53E (GCA→GAG) coaDpantetheine‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (2/0);  new base A (17/23);  total (19/23)
Fisher's exact test for biased strand distribution p-value = 1.99e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.63e-01
*NC_0009133,809,9830AG97.6% 131.2 / ‑4.3 42A53E (GCA→GAG) coaDpantetheine‑phosphate adenylyltransferase
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (18/23);  total (19/23)
Fisher's exact test for biased strand distribution p-value = 4.52e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.28e-01

GTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGGTTGTTATGCAAAAACGGGCGATTTATCCGGGTACTTTCGATCCCATTACCAATGGTCATATCGATATCGTGACGCGCGCCACGCAGATGTTCGATCACGTTATTCTGGCGATTGCCGCCAGCCCCAGTAAAAAACCGATGTTTACCCTGGAAGAGCGTGTGGCA‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGTCGGGTTTAGTGATTTAATGGCGAACTTCGCCCGTAATCAACACGCTACGGTGCTGATTCGTGGCCTGCGTGCGGTGGCAGATTTTGAATATGAAATGCAGCTGGCGCATATGAATCGCCACTTAATGCCGGAACTGGAAAGTGTGTTTCTGATGCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCGTT  >  NC_000913/3809740‑3810216
                                                                                                                                                                                                                                                  ||                                                                                                                                                                                                                                          
gTACTGTCTCAAAACCAGGGCGCGCTAAAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGGTTGTTATGCAAAAACGGGCGATTTATCCGGGTACTTTCGATCCCATTACAAATGGTCATATCGATATCGTGACGCGCGCCACGCAAATGTTCGATCACGTTATTCTGGCGCTTGCCGCCAGCCCCAGTAAAAAACCGATGTTTACACTGGAAGAGCGTGtggat‑ctggcac                                                                                                                                                                                                                                    >  2:191202/1‑251 (MQ=255)
               cAGGGCGCGCTACAGCGTATGCTTCAACTGCTGGAACCTTCCCTGCCACCGACAACGCATTGAGGTTGTTATGCAAAAACGGGCGATTTATCCGGGTACTTTCGATCCCATTACCAATGGTCATATCGATATCGTGACGCGCGCCACGCAAATGTTCGATCACGTTATTCTGCCGATTGCCGCCAGCCCCAGTAAAAAACCGATGTTTACCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCgcgc                                                                                                                                                                                                                     >  2:83334/1‑251 (MQ=255)
                          aCAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGGTTGTTATGCAAAAACGGGCGATTTATCCGGGTACTTTCGAGCCCATTACCAATGGTCATATCGATATCGTGACGCGCGCCACGCAAATGTGCGATCACGTTATTCTGGCGAGTGCCGCCAGCCCCAGGAAAAAACCGATTTTTGCCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAGCCGCGCATCTGGGGaa                                                                                                                                                                                                           <  2:437273/250‑1 (MQ=255)
                                   gCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGGTTGTTATGCAAAAACGGGCGATTTATCCGGGTACTTTCGATCCCATTACCAATGGTCATATCGATATCGTGACGCGCGCCACGCAAATGTTCGATCACGTTATTCTGGCGATTGCCGCCAGCCCCAGTAAAAAACCGATGTTTACCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTg                                                                                                                                                                                                 <  2:372524/251‑1 (MQ=255)
                                   gCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGGTTGTTATGCAAAAACGGGCGATTTAGCCGGGTACTTTCGATCCCATTACCAATGGTCATATCGATATCGTGACGCGCGCCACGCAAATGTTCGATCACGTTATTCTGGCGATTGCCGCCAGCCCCAGTAAAAAACCGAGGTGTACCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGaa                                                                                                                                                                                                    <  2:259727/248‑1 (MQ=255)
                                        aaCTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGGTTGTTATGCAAAAACGGGCGATTTATCCGGGTACTTTCGATCCCATTACCAATGGTCATATCGATATCGTGACGCGCGCCACGCAAATGTTCGATCACGTTATTCTGGCGATTGCCGCCAGCCCCAGTAAAAAACCGATGTTTACCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGTCgg                                                                                                                                                                                            >  1:67947/1‑251 (MQ=255)
                                               ggAACCTTACCTGCCACCGAAAACGCATTGAGGTTGTTATGCAAAAACGGGCGATTTATCCGGGTACTTTCGATCCCATTACCAATGGTCATATCGATATCGTGACGCGCGCCACGCAAATGTTCGATCACGTTATTCTGGCGATTGCCGCCAGCCCCAGTAAAAAACCGATGTTTACCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGTCGGGTTTAGt                                                                                                                                                                                     <  2:150687/251‑1 (MQ=255)
                                                  aCCTTACCTGCCACCGAAAACGCATTGAGGTTGTTATGCAAAAACGGGCGATTTATCCGGGTACTTTCGCTCCCATTACCAATGGTCATATCGATATCGTGACGCGCGCCACGCAAATGTTCGATCACGTTATTCTGGCGATTGCCGCCAGCCCCAGTAAAAAACCGATGTTTACCCTGGCAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGGGGTCGGGTTTAGTGAt                                                                                                                                                                                  >  1:260703/1‑251 (MQ=255)
                                                       aCCTGCCACCGAAAACGCATTGAGGTTGTTATGCAAAAACGGGCGATTTATCCGGGTACTTTCGATCCCATTACCAATGGTCATATCGATATCGTGACGCGCGCCACGCAAATGTTCGATCACGTTATTCTGGCGATTGCCGCCAGCCCCAGTAAAAAACCGATGTTTACCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGTCTGGTTTAGTGAtttgat                                                                                                                                                                             <  1:305750/251‑4 (MQ=255)
                                                           gCCACCGAAAACGCATTGAGGTTGTTATGCAAAAACGGGCGATTTATCCGGGTACTTTCGATCCCATTACCAATGGTCATATCGATATCGTGACGCGCGCCACGCAAATGTTCGATCACGTTATTCTGGCGATTGCCGTCAGCCCCAGTAAAAAACCGATGTTTACCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGTCGGGTTTAGTGATTTGATGGc                                                                                                                                                                          <  1:393255/250‑1 (MQ=255)
                                                              aCCGAAAACGCATTGAGGTTGTTATGCAAAAACGGGCGATTTATCCGGGTACTTTCGCTCCCATTACCAATGGTCATATCGATATAGTGACGCACGCCACGAAAATGTTCGATCACGTTATTCTGGCGATTGCCGCCAGCCGCAGTAAACAACCGATGTTTACCCTGGAGGAGCGTGTG‑CAGCTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGGGGTCGGGTTTatgggtttgatagcgaac                                                                                                                                                                      >  1:71956/1‑234 (MQ=255)
                                                              aCCGAAAACGCATTGAGGTTGTTATGCAAAAACGGGCGATTTATCCGGGTACTTTCGATCCCATTACCAATGGTCATATCGATATCGTGACGCGCGCCACGCAAATGTTCGATCACGTTATTCTGGCGATTGCCGCCAGCCCCAGTAAAAAACCGATGTTTACCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGTCGGGTTTAGTGATTTGATGGCGAAc                                                                                                                                                                      >  1:15166/1‑251 (MQ=255)
                                                                 gAAAACGCATTGAGGTTGTTATGCAAAAACGGGCGATTTATCCGGGTACTTTCGATCCCATTACCAATGGTCATATCGATATCGTGACGCGCGCCACGCAAATGTTCGATCACGTTATTCTGGCGATTGCCGCCAGCCCCAGTAAAAAACCGATGTTTACCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGTCGGGTTTAGTGATTTGATGGCGAACTTc                                                                                                                                                                   <  2:303161/251‑1 (MQ=255)
                                                                    aaCGCATTGAGGTTGTTATGCAAAAACGGGCGATTTATCCGGGTACTTTCGATCCCATTACCAATGGTCATATCGATATCGTGACGCGCGCCACGCAAATGTTCGATCACGTTATTCTGGCGATTGCCGCCAGCCCCAGTAAAAAACCGATGTTTACCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGTCGGGTTTAGTGATTTGATGGCGAACTTCGcc                                                                                                                                                                >  1:290297/1‑251 (MQ=255)
                                                                     aCGCATTGAGGTTGTTATGCAAAAACGGGCGATTTATCCGGGTACTTTCGATCCCATTACCAATGGTCATATCGATATCGTGACGCGCGCCACGCAAATGTTCGATCACGTTATTCTGGCGATTGCCGCCAGCCCCAGTAAAAAACCGATGTTTACCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGTCGGGTTTAGTGATTTGATGGCGAACTTCGccc                                                                                                                                                               <  1:159404/251‑1 (MQ=255)
                                                                            gAGGTTGTTATGCAAAAACGGGCGATTTATCCGGGTACTTTCGATCCCATTACCAATGGTCATATCGATATCGTGACGCGCGCCACGCAAATGTTCGATCACGTTATTCTGGCGATTGCCGCCAGCCCCAGTAAAAAACCGATGTTTACCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGTCGGGTTTAGTGATTTGATGGCGAACTTCGCCCGTAATc                                                                                                                                                         <  1:204809/250‑1 (MQ=255)
                                                                              ggTTGTTATGCAAAAACGGGCGATTTATCCGGGTACTTTCGATCCCATTACCAATGGTCATATCGATATCGTGACGCGCGCCACGCAAATGTTCGATCACGTTATTCTGGCGATTGCCGCCAGCCCCAGTAAAAAACCGATGTTTACCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGTCGGGTTTAGTGATTTGATGGCGAACTTCGCCCGTAATCAac                                                                                                                                                      >  2:416457/1‑251 (MQ=255)
                                                                                                          ccGGGTACTTTCGATCCCATTACCAATGGTCATATCGATATCGTGACGCGCGCCACGCAAATGTTCGATCACGTTATTCTGGCGCTTGCCGCCAGCCCCAGTAAAAAACCGATGTTTACCCTGGAAGAGCATGTGGAG‑CTGGCACAGCAGGCAACCGCGCCTCTGGGGAACGTGGAAGTGGTCGGGTTTAGTGATTTGATGGCGAACTTCGCCCGTAATCAACACGCTACGGTGCTGATTCGTGGCCtgcg                                                                                                                          >  1:432616/1‑251 (MQ=255)
                                                                                                              gTACTTTCGATCCCATTACCACTGGTCATATCGATATCGTGACGCGCGCCACGAAAATGTTCGATCACGTTATTCTGGCGATTGCCGACAGCCCCAGTAAAAAACCGATGTTTCCCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGTCGGGTTTAGTGATTTGATGGCGAACTTCGCCCGTAATCAACACGCTACGGTGCTGATTCGTGGCCtgcgggcg                                                                                                                      >  2:428388/1‑251 (MQ=255)
                                                                                                                         cccATTACCAATGGTCATATCGATATCGTGACGCGCGCCACGCAAATGTTCGATCACGTTATTCTGGCGATTGCCGCCAGCCCCAGTAAAAAACCGATGTTTACCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGTCGGGTTTAGTGATTTGATGGCGAACTTCGCCCGTAATCAACACGCTACGGTGCTGATTCGTGGCCTGCGTGCGGTGGCGGAttt                                                                                                           <  1:31728/251‑1 (MQ=255)
                                                                                                                         cccATTACCAATGGTCATATCGATATCGTGACGCGCGCCACGCAAATGTTCGATCACGTTATTCTGGCGATTGCCGCCAGCCCCAGTAAAAAACCGATGTTGACCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGTGTAACGTGGAAGTGGTCGGGTTTAGTGATTTGATGGCGAACTTCGCCCGTAATCAACACGCTACGGTGCTGATTCGTGGCCTGCGTGCGGTGGCGGAtt                                                                                                            <  1:248234/250‑1 (MQ=255)
                                                                                                                            aTTACCAATGGTCATATCGATATCGTGACGCGCGCCACGCAAATGTTCGATCACGTTATTCTGGCGATTGCCGCCAGCCCCAGTAAAAAACCGATGTTTACCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGTCGGGTTTAGTGATTTGATGGCGAACTTCGCCCGTAATCAACACGCTACGGTGCTGATTCg                                                                                                                                   >  2:214643/1‑224 (MQ=255)
                                                                                                                            aTTACCAATGGTCATATCGATATCGTGACGCGCGCCACGCAAATGTTCGATCACGTTATTCTGGCGATTGCCGCCAGCCCCAGTAAAAAACCGATGTTTACCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGTCGGGTTTAGTGATTTGATGGCGAACTTCGCCCGTAATCAACACGCTACGGTGCTGATTCg                                                                                                                                   <  1:214643/224‑1 (MQ=255)
                                                                                                                                   aTGGTCATATCGATATCGTGACGCGCGCCACGCAAATGTTCGATCACGTTATTCTGGCGATTGCCGCCAGCCCCAGTAAAAAACCGATGTTTACCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGTCGGGTTTAGTGATTTGATGGCGAACTTCGCCCGTAATCAACACGCTACGGTGCTGATTCGTGGCCTGCGTGCGGTGGCGGATTTTGAATATGa                                                                                                  >  2:253830/1‑250 (MQ=255)
                                                                                                                                       tCATATCGATATCGTGACGCGCGCCACGCAAATGTTCGATCACGTTATTCTGGCGATTGCCGCCAGCCCCAGTAAAAAACCGATGTTTACCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGTCGGGTTTAGTGATTTGATGGCGAACTTCGCCCGTAATCAACACGCTACGGTGCTGATTCGTGGCCTGGGTGCGGTGGCGGATTTTGAATATGAAATg                                                                                              >  2:93403/1‑250 (MQ=255)
                                                                                                                                                       aCGCGCGCCACGCAAATGTTCGATCACGTTATTCTGTCGATTGCCGCCAGCCCCAGTAAAAAACCGATGTTTTCCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGTCGGGTTTAGTGAGTTGATGGCGAACTTCGCCCGTAATCAACACGCTACGGTGCTGATTCGTGGCCTGCGTGCGGTGGCGGATTTTGAATATGAAATGCAGCTGGCGCATATGaa                                                                             <  1:321214/251‑1 (MQ=255)
                                                                                                                                                           gcgcCACGCAAATGTTCGATCACGTTATTCTGGCGATTGCCGCCAGCCCCAGTAAAAAACCGATGTTTACCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGTCGGGTTTAGTGATTTGATTGCGAACTTCGCCCGTAATCAACACGCTACGGTGCTGATTCGTGGCCTGCGTGCTGTGGCGGATTTTGAATATGAAATGCAGCGGGCGCATATGAATCGc                                                                         <  1:358669/251‑1 (MQ=255)
                                                                                                                                                                  gcaaaTGTTCGATCACGTTATTCTGGCGATTGCCGCCAGCCCCAGTAAAAAACCGATGTTTACCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCCACCGCGCATCTGGGGAACGTGGAAGTGGTCGGGTTTAGCGATTTGATGGCGAACTTCGCCCGTAATCAACACGCTACGGTGCTGATTCGTGGCCTGCGTGCGGTGGCGGATTTTGAATATGACATGCAGCTGGCGCATATGAATCGCCACTTa                                                                   >  1:155049/5‑250 (MQ=255)
                                                                                                                                                                    aaaTGTTCGATCACGTTATTCTGGCGATTGCCGCCAGCCCCAGTAAAAAACCGATGTTTACCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGTCGGGTTTAGTGATTTGATGGCGAACTTCGCCCGTAATCAACACGCTACGGTGCTGATTCGTGGCCTGCGTGCGGTGGCGGATTTTGAATATGAAATGCAGCTGGCGCATATGAATCGCCACTTAATg                                                                <  2:6976/249‑1 (MQ=255)
                                                                                                                                                                       tGTTCGATCACGTTATTCTGGCGATTGCCGCCAGCCCCAGTAAAAAACCGATGTTTACCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGTCGGGTTTAGTGATTTGATGGCGAACTTCGCCCGTAATCAACACGCTACGGTGCTGATTCGTGGCCTGCGTGCGGTGGCGGATTTTGAATATGAAATGCAGCTGGCGCATATGAATCGCCACTTAATGCCg                                                             >  1:280170/1‑251 (MQ=255)
                                                                                                                                                                           cGATCACGTTATTCTGGCGATTGCCGCCAGCCCCAGTAAAAAACCGATGTTTACCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGCAGTGGTCGGGTTTAGTGATTTGATGGCGAACTTCGCCCGTAATCAACACGCTACGGTGCTGATTCGTGGCCTGCGTGCGGTGGCGGATTTTGAATATGAAATGCAGCTGGCGCATATGAATCGCCACTTAATGCCGGAAc                                                         >  2:275341/1‑251 (MQ=255)
                                                                                                                                                                             aTCACGTTATTCTGGCGATTGCCGCCAGCCCCAGTAAAAAACCGATGTTTACCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGTCGGGTTTAGTGATTTGATGGCGAACTTCGCCCGTAATCAACACGCTACGGTGCTGATTCGTGGCCTGCGTGCGGTGGCGGATTTTGAATATGAAATGCAGCTGGCGCATATGAATCGCCACTTAATGCCGGAACTg                                                       <  1:185606/251‑1 (MQ=255)
                                                                                                                                                                                 cggtatgcTGGGGATTGCCGCCAGCCCCAGTAAAAAACCGATGTTTACCCTGGAAGGGGGTGTGGAG‑ATGGCACAGCATGCAACCGCGCATCTGGGGAACGTGGAAGTGGCCGGGTTTAGTGATTTGATGGCGAACTTCGCACGTAATCAACACGCTACGGTGCTGATTCGTGGCATGCGTGCGGTGGCGTATTTTGAATATTAAATGCAGCTGGCGCATATGAATCGCCACTTAATGCCGGAACTGGaaa                                                   <  1:236263/244‑1 (MQ=255)
                                                                                                                                                                                         tGGCGATTGCCGCCAGCCCCAGTAAAAAACCGATGTTTACCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGTCGGGTTTAGTGATTTGATGGCGAACTTCGCCCGTAATCAACACGCTACGGTGCTGATTCGTGGCCTGCGTGCGGTGGCGGATTTTGAATATGAAATGCAGCTGGCGCATATGAATCGCCACTTAATGCAGGAACTGGAAAgtgtgatt                                           >  1:414933/1‑248 (MQ=255)
                                                                                                                                                                                         tGGCGATTGCCACCAGCCCCAGTACAACCCCGATGTTTACCCTGGAAGAGCGTGTGGCG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGTCGGGTTTAGTGATTTGATGGCGCACTTCGCCCGTAATCAACACGCTACGGTGCTGATTCGTGGCCTGCGTGCGGTGGCGGATTTTGAATATGAAATGCAGCTGGCGCATATGAATCGCCACTTAATGCCGGAACTGGAAAGTGTGttt                                           >  2:221957/1‑251 (MQ=255)
                                                                                                                                                                                                                 aaaaaCCGATGTTTACCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGTCGGGTTTAGTGATTTGATGGCGAACTTCGCCCGTAATCAACACGCTACGGTGCTGATTCGTGGCCTGCGTGCGGTGGCGGATTTTGAATATGAAATGCAGCTGGCGCATATGAATCGCCACTTAATGCCGGAACTGGAAAGTGTGTTTCTGATGCCGTCGAAAGAGTGGTCg                   <  1:191202/251‑1 (MQ=255)
                                                                                                                                                                                                                    aaCCGATGTTTACCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGTCTGGTTTAGTGATTTGATGGCGAACTTCGCCCGTAATCAACACGCTACGGTGCTGATTCGTGGCCTGCGTGCGGTGGCGGATTTTGAATATGAAATGCAGCTGGCGCATATGAATCGCCACTGAATGCCGGAACTGGAAAGTGTGTTTCTGATGCCGTCGAAAGAGTGGTCGttt                <  1:423202/251‑1 (MQ=255)
                                                                                                                                                                                                                            tttACCCTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGCCGGGTTTAGTGATTTGATGGCGAACTTCGCCCGTAATCAACACGCTACGGTGCTGCTTCGTGGCCTGCGTGCGGTGGCGGATTTTGAATATGAAATGCAGCTGGCGCATATGAATCGCCACTTAATGCCGGAACTGGAAAGTGTGTTTCTGATGCCGTCGAAAGAGTGGTCGTTTATCTCTTc        >  1:31194/1‑251 (MQ=255)
                                                                                                                                                                                                                                cccTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGTCGGGTTTAGTGATTTGATGGCGAACTGCGCCCGTAATCAACACGCTACGGTGCTGATTCGTGGCCTGCGTGCGGTGGCGGATGTTGAATATGAAATGCAGCTGGCGCATATGAATCGCCACTTAATGCCGGAACTGGAAAGTGTGTTTCTGATGCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCg    <  2:87617/251‑1 (MQ=255)
                                                                                                                                                                                                                                cccTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGGCGGGTTTAGTGATTTGATGGCGAACTTCGCCCGTAATCAACACGCTACGGTGCTGATTCGTGGCCTGCGTGCGGTGGCGGATTTTGAATATGAAATGCAGCTGGCGCATATGAATCGCCACTTAATGCCGGAACTGGAAAGTGTGTTTCTGATGCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCg    <  2:290297/251‑1 (MQ=255)
                                                                                                                                                                                                                                 ccTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGTCGGGTTTAGTGATTTGATGGCGAACTTCGCCCGTAATCAACACGCTACGGTGCTGATTCGTGGCCTGCGTGCGGTGGCGGATTTTGAATATGAAATGCAGCTGGCGCATATGAATCGCCACTTAATGCCGGAACTGGAAAGTGTGTTTCTGATGCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCGt   <  1:221957/251‑1 (MQ=255)
                                                                                                                                                                                                                                  cTGGAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGTCGGGTTTAGTGATTTGATGGCGAACTTCGCCCGTAATCAACACGCTACGGTGCTGATTCGTGGCCTGCGTGCGGTGGCGGATTTTGAATATGAAATGCAGCTGGCGCATATGAATCGCCACTTAATGCCGGAACTGGAAAGTGTGTTTCTGATGCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCGtt  >  1:221566/1‑251 (MQ=255)
                                                                                                                                                                                                                                    ggAAGAGCGTGTGGAG‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGTCGGGTTTAGTGATTTGATGGCGAACTTCGCCCGTAATCAACACGCTACGGTGCTGATTCGTGGCCTGCGTGCGGTGGCGGATTTTGAATATGAAATGCAGCTGGCGCATATGAATCGCCACTTAATGCCGGAACTGGAAAGTGTGTTTCTGATGCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCGtt  <  1:5670/249‑1 (MQ=255)
                                                                                                                                                                                                                                                  ||                                                                                                                                                                                                                                          
GTACTGTATCAAAACCAGGGCGCGCTACAGCGTCTGCTTCAACTGCTGGAACCTTACCTGCCACCGAAAACGCATTGAGGTTGTTATGCAAAAACGGGCGATTTATCCGGGTACTTTCGATCCCATTACCAATGGTCATATCGATATCGTGACGCGCGCCACGCAGATGTTCGATCACGTTATTCTGGCGATTGCCGCCAGCCCCAGTAAAAAACCGATGTTTACCCTGGAAGAGCGTGTGGCA‑CTGGCACAGCAGGCAACCGCGCATCTGGGGAACGTGGAAGTGGTCGGGTTTAGTGATTTAATGGCGAACTTCGCCCGTAATCAACACGCTACGGTGCTGATTCGTGGCCTGCGTGCGGTGGCAGATTTTGAATATGAAATGCAGCTGGCGCATATGAATCGCCACTTAATGCCGGAACTGGAAAGTGTGTTTCTGATGCCGTCGAAAGAGTGGTCGTTTATCTCTTCATCGTT  >  NC_000913/3809740‑3810216

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