Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,886,749 G→A intergenic (+27/‑79) yidC → / → mnmE membrane protein insertase YidC/5‑carboxymethylaminomethyluridine‑tRNA synthase GTPase subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,886,7490GA97.4% 126.1 / NA 40intergenic (+27/‑79)yidC/mnmEmembrane protein insertase YidC/5‑carboxymethylaminomethyluridine‑tRNA synthase GTPase subunit
Reads supporting (aligned to +/- strand):  ref base G (0/1);  new base A (20/18);  total (21/19)
Fisher's exact test for biased strand distribution p-value = 4.87e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TCATTCAGAAGATGTCGCCGACCACAGTGACCGACCCGATGCAGCAGAAGATCATGACCTTTATGCCGGTCATCTTCACCGTGTTCTTCCTGTGGTTCCCGTCAGGTCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGAC  >  NC_000913/3886505‑3886992
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
tcattcaGAAGATGTCGCCGACCACAGTGACCGACCCGATGCAGCAGAAGATCATGACCTTTATGCCGGTCATCTTCACCGTGTTCTTCCTGTGGTTCCCGTCAGGTCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGACGGTc                                                                                                                                                                                                                                                >  1:44476/1‑250 (MQ=255)
 cattcaGAAGATGTCGCCGACCACAGTGACCGACCCGATGCAGCAGAAGATCATGACCTTTATGCCGGTCATCTTCACCGTGTTCTTCCTGTGGTTCCCGTCAGGTCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGACGGTCa                                                                                                                                                                                                                                               <  2:209635/250‑1 (MQ=255)
      agaagaTGTCGCCGACCACAGTGACCGACCCGATGCAGCAGAAGATCATGACCTTTATGCCGGTCATCTTCACCGTGTTCTTCCTGTGGTTCCCGTCAGGTCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGTTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAGAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGACGGTCAGTTGAc                                                                                                                                                                                                                                         <  1:255975/251‑1 (MQ=255)
             gTCGCCGACCACAGTGACCGACCCGATGCAGCAGAAGATCATGACCTTTATGCCGGTCATCTTCACCGTGTTCTTCCTGTGGTTCCCGTCAGGTCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTGTTCGCTAAAATAAGGACGGTCAGTTGACCGCAttt                                                                                                                                                                                                                                  <  2:289426/251‑1 (MQ=255)
                   gACCACAGTGACCGACCCGATGCAGCAGAAGATCATGACCTTTATGCCGGTCATCTTCACCGTGTTCTTCCTGTGGTTCCCGTCAGGTCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGACGGTCAGTTGACCGCCTTttttctt                                                                                                                                                                                                                            >  1:162070/1‑251 (MQ=255)
                       acagtgcgcgaaccGATGCAGCAGAAGATCATGACATATATGCCGGTCATATTCACCGTGTCCTTCCTGTGGTTCCCGTCAGGTCTAGTGCTGTACTATATCGTCACCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAACATCCTGATTCGGTGAGTTTTCGCTAAAATACGGCCGGTCAGTTGACCGCCTTTTTTCTTTTCg                                                                                                                                                                                                                        >  2:160567/13‑251 (MQ=255)
                         aGTGACCGACCCGATGCAGCAGAAGATCATGACCTTTATGCCGGCCATCTTCACCGTGTTCTTCCTGTGGTTCCCGTCAGGTCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCAGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGACGGTCAGTTGACCGCCTTTTTTCTTTTCGt                                                                                                                                                                                                                       >  2:195179/1‑250 (MQ=255)
                               cgaccCGATGCAGCAGAAGATCATGACCTTTATGCCGGTCATCTTCACCGTGTTCTTCCTGTGGTTCCCGTCAGGTCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGACGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCgg                                                                                                                                                                                                                >  2:69958/1‑251 (MQ=255)
                                             agaagaTCATGACCTTTATGCCGGTCATCTTCACCGTGTTCTTCCTGTGGTTCCCGTCAGGTCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGACGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCTCCGCGCa                                                                                                                                                                                                  <  2:39355/251‑1 (MQ=255)
                                             agaagaTCATAACCTTTATGCCGGTCAACTTCACCGTGGTCTTCCTGTGGATCCCGTAAGGTCTGGTGCTGTACTATATCGTCAGCAAACTGGTAACCATTATTCAGCCGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCGTGATTCGGTGAGTTTTCGCTAAAATAAGGACGGTAAGTTGACCGCCTTTTTTCTTTTAGTAGGGCGGATAAGCACCgcgc                                                                                                                                                                                                   >  2:247615/1‑250 (MQ=255)
                                                agaTCATGACCTTTATGCCGGTCATCTTCACCGTGTTCTTCCTGTGGTTCCCGTCAGGGCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGACGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATcc                                                                                                                                                                                               <  2:101066/251‑1 (MQ=255)
                                                             tATGCCGGTCATCTTCACCGTGTTCTTCCTGTGGTTCCCGTCAGGTCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGATTTTTCGCTAAAATAAGGACGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGc                                                                                                                                                                                   <  1:324683/250‑1 (MQ=255)
                                                              aTGCCGGTCTTCGTCACCGTGTTCTTCCTGTGGTTCCCGTCAGGTCTGGTGCTGTACTATATCGTCAGCATCCTGGTAACCATTATTCAGCAACAGCTGATTTACCGGGGTCAGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAACCCTGATTCGGTGAGTTTTCGCTAAAATAAGGATGGTCAGTTGACTGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCaa                                                                                                                                                                                 <  1:244815/251‑1 (MQ=255)
                                                              aTGCCGGTCATCTTCACCGTGTTCTTCCTGTGGTTCCCGTCAGGTCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGACGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCa                                                                                                                                                                                  >  2:358092/1‑250 (MQ=255)
                                                                 ccGGTCATCTTCACCGTGTTCTTCCTGTGGTTCCCGTCAGGTCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGACGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAAc                                                                                                                                                                                >  2:344496/1‑249 (MQ=255)
                                                                      cATCTTCACCGTGTTCTTCCTGTGGTTCCCGTCAGGTCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGACGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAAc                                                                                                                                                                          >  1:14766/1‑250 (MQ=255)
                                                                                   ttnttcCTGTGGTTCCCGTCAGGTCTGGTGCTATACTATATCGTCAGCAACCTGGTAACCATTAATCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCGGCATAGCCGCGAGAAGAAAAAATCCGGATTCGGGGAGTTTTCGCTAAAATAAGGACGGTCAGTGGACCGCCTTTTTTCTTTTCGTAGGGCGTATAAGCACCGCGTATCCGACACACAAAGCAACAGGAACATCATGAGCGATa                                                                                                                                                             <  2:380816/250‑1 (MQ=255)
                                                                                              gTTCCCGTCAGGTCTGGTGCTGTACTATATCGGCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCTCTAAAATAAGGACGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCg                                                                                                                                                 <  2:244439/251‑1 (MQ=255)
                                                                                                   cGTCAGGTCCGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGACGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGcc                                                                                                                                            >  1:381119/1‑251 (MQ=255)
                                                                                                      cAGGTCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGACGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAg                                                                                                                                         <  2:323972/251‑1 (MQ=255)
                                                                                                       aGGTCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGCCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTACAATAAGGACGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAgg                                                                                                                                        >  2:270824/1‑251 (MQ=255)
                                                                                                                gCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGACGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCt                                                                                                                                <  1:323885/250‑1 (MQ=255)
                                                                                                                   gTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGACGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGACACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCACGGCCACGCCTCCgg                                                                                                                            >  2:123346/1‑251 (MQ=255)
                                                                                                                          atCGTAAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGACTGGAAAAACGTGGCCTGAATAGCCGCGAGAAGAAAAAATCCTGATTCGGAGAGTTTTCGCTAAAATAAGGACGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCACAGGCCAAGCCTCCGGGACGTgg                                                                                                                     >  2:60554/1‑251 (MQ=255)
                                                                                                                               cAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCTGTGAGTTTTCGCTAAAATAAGGACGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGTCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCAGGGACGTGgcggcg                                                                                                                >  1:218680/1‑251 (MQ=255)
                                                                                                                                   aaCCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGACTGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGCCGGCGTTgg                                                                                                            >  1:147322/1‑251 (MQ=255)
                                                                                                                                    aCCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGACGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGc                                                                                                           >  1:333439/1‑251 (MQ=255)
                                                                                                                                      cTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGACGGTAAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGAGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCAt                                                                                                         >  1:42695/1‑251 (MQ=255)
                                                                                                                                                  attCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGACGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCAtctc                                                                                             <  2:44476/251‑1 (MQ=255)
                                                                                                                                                     cangagcagcTGATTTACGGTGGTCTGAAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTGCGCTAAAATAAGGACGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGTCGGCGTTGGCATCCTTCGCATCTCCg                                                                                           <  2:374705/246‑1 (MQ=255)
                                                                                                                                                        acgatgctaaTTTACGGTGGTCTGGAATAACGAGGCCTGCTTAGCCTGGAGAAGAAATAATCCTTATTCGGTGTGTTTTCGCTAAAATAAGGAGGGTCAGTTGACCGCTTTGTTTCTTTTCGGAGGGCGGATAAGCTCCGCGCATCCGCCGCACAAAGCAACAGGAACCTCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGTACGTGGCGGCGTTGGCTTCCTTAGCATCTCCGGCtt                                                                                       <  2:46122/242‑1 (MQ=255)
                                                                                                                                                         agcagcTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGACGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCGGCTTc                                                                                      >  2:358653/1‑251 (MQ=255)
                                                                                                                                                            agcTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGACGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCGGCTTCaa                                                                                    >  1:376204/1‑250 (MQ=255)
                                                                                                                                                                     aCCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGACGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCGGCTTCAAAGCCCGTGa                                                                           >  2:123886/1‑250 (MQ=255)
                                                                                                                                                                                              gCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGACGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTaa                                                   <  2:14766/249‑1 (MQ=255)
                                                                                                                                                                                                    ccGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGACGGTCAGTTGACCGCCTTTTTTCTTTGCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTa                                           <  2:165444/251‑1 (MQ=255)
                                                                                                                                                                                                       ggagaAGAAAAGATCCTGATTCGGTCTGTTTCGGCCAAAAGAAGGGCGGTCATGTGACCGCCTTTTTTCTGTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAGCAGGAGCATCGTGAGCGATAATGACGCTAGCGGAGCCCAGGCCACGCCGCCGGGTCGTGGCGGCGTTGGCATCCTGCGCATCTCCGGCTTCAAAGCCCGTGAAGTTGCCGCAACCGTGCTGGGTAAACTGCCTAAGc                                        <  2:101422/250‑1 (MQ=255)
                                                                                                                                                                                                       cgAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGACGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGc                                        <  1:356682/251‑1 (MQ=255)
                                                                                                                                                                                                                                         cTAAAATAAGGACGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGTACGTGGCGGCGTTGGCATCCTGCGCATCTCCGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGc      >  2:112954/1‑251 (MQ=255)
                                                                                                                                                                                                                                             aaTAAGGACGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGAc  <  1:150558/251‑1 (MQ=255)
                                                                                                                                                                                                                                                    |                                                                                                                                                                                                                                                   
TCATTCAGAAGATGTCGCCGACCACAGTGACCGACCCGATGCAGCAGAAGATCATGACCTTTATGCCGGTCATCTTCACCGTGTTCTTCCTGTGGTTCCCGTCAGGTCTGGTGCTGTACTATATCGTCAGCAACCTGGTAACCATTATTCAGCAGCAGCTGATTTACCGTGGTCTGGAAAAACGTGGCCTGCATAGCCGCGAGAAGAAAAAATCCTGATTCGGTGAGTTTTCGCTAAAATAAGGGCGGTCAGTTGACCGCCTTTTTTCTTTTCGTAGGGCGGATAAGCACCGCGCATCCGCCACACAAAGCAACAGGAACATCATGAGCGATAATGACACTATCGTAGCCCAGGCCACGCCTCCGGGACGTGGCGGCGTTGGCATCCTGCGCATCTCCGGCTTCAAAGCCCGTGAAGTTGCCGAAACCGTGCTGGGTAAACTGCCTAAGCCGCGCTACGCCGATTATCTTCCGTTTAAAGACGCCGAC  >  NC_000913/3886505‑3886992

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