Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 790,931 C→G 9.1% W33C (TGG→TGC)  galT ← galactose‑1‑phosphate uridylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913790,9310CG 9.1% 222.2 / 12.0 66W33C (TGG→TGC) galTgalactose‑1‑phosphate uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (28/32);  new base G (4/2);  total (32/34)
Fisher's exact test for biased strand distribution p-value = 4.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.24e-01

CAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACTGCC  >  NC_000913/790784‑791073
                                                                                                                                                   |                                                                                                                                               
cAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGc                                                                                                                                                 >  2:604878/1‑148 (MQ=255)
  aaGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCa                                                                                                                                               >  1:385128/1‑148 (MQ=255)
   aGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑g                                                                                                                                             >  1:908158/1‑148 (MQ=255)
        cAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGgcgc                                                                                                                                       >  3:592167/1‑149 (MQ=255)
        cAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGgcgc                                                                                                                                       <  4:150501/149‑1 (MQ=255)
         aaaGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCTGGCAG‑GGGgcgca                                                                                                                                      <  3:203623/149‑2 (MQ=255)
         aaaGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGgcgc                                                                                                                                       <  1:781629/148‑1 (MQ=255)
          aaGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCTGGCAG‑GGGgcgc                                                                                                                                       >  1:748723/1‑147 (MQ=255)
          aaGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCG‑CCCCTGGCA‑gggggcgc                                                                                                                                      <  2:748723/147‑5 (MQ=255)
            gTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAg                                                                                                                                   <  3:737827/149‑1 (MQ=255)
             tCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAg                                                                                                                                   >  2:637643/1‑148 (MQ=255)
                ttAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGcgcg                                                                                                                               >  4:723427/1‑149 (MQ=255)
                    tGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCgcggtgc                                                                                                                           >  4:642895/1‑149 (MQ=255)
                         aCGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCG‑CCCCTGGCA‑gggggcg                                                                                                                                       <  1:121822/131‑4 (MQ=255)
                         aCGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCG‑CCCCTGGCA‑gggggcg                                                                                                                                       >  2:121822/1‑128 (MQ=255)
                              aGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAAcc                                                                                                                  <  1:604878/148‑1 (MQ=255)
                                tCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGa                                                                                                               >  3:450385/1‑149 (MQ=255)
                                  ccGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGaa                                                                                                              <  1:837807/148‑1 (MQ=255)
                                  ccGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAAt                                                                                                             >  4:797180/1‑149 (MQ=255)
                                  ccGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAAt                                                                                                             <  3:642895/149‑1 (MQ=255)
                                   cGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAAcc                                                                                                                  >  2:727451/1‑143 (MQ=255)
                                   cGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAAcc                                                                                                                  <  1:727451/143‑1 (MQ=255)
                                      tgtAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCAt                                                                                                         <  3:723427/149‑1 (MQ=255)
                                           tCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGcc                                                                                                     <  2:235657/148‑1 (MQ=255)
                                               ggTTTTTATCGCCTGTCACCCGTACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTg                                                                                                <  4:746242/149‑1 (MQ=255)
                                                  ttttATCGCCTGTCACCCGCACATTACCTGCGCATAGGAAGCAATCTGGATCGTGCTCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCT‑GGGGCGCCTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGc                                                                                             <  4:557701/149‑1 (MQ=255)
                                                       tCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCTGGCAG‑GGGgcgcag                                                                                                                                     >  3:862245/1‑102 (MQ=255)
                                                       tCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGcgcg                                                                                                                               <  4:615055/110‑1 (MQ=255)
                                                       tCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGcgcg                                                                                                                               >  3:615055/1‑110 (MQ=255)
                                                       tCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGtt                                                                                        <  3:95177/149‑1 (MQ=255)
                                                       tCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGtt                                                                                        <  4:450385/149‑1 (MQ=255)
                                                                     cacaTTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGAtgtg                                                                          <  3:171522/149‑1 (MQ=255)
                                                                        aTTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGgat                                                                       <  4:110237/149‑1 (MQ=255)
                                                                        aTTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGga                                                                        <  2:908158/148‑1 (MQ=255)
                                                                                cgcAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGAtgtg                                                                          >  3:97419/1‑138 (MQ=255)
                                                                                cgcAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGAtgtg                                                                          <  4:97419/138‑1 (MQ=255)
                                                                                   agagGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACggg                                                             >  1:712803/1‑148 (MQ=255)
                                                                                       gAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTaa                                                        >  3:876310/1‑149 (MQ=255)
                                                                                        aaGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCgg                                                                                           >  2:868755/1‑113 (MQ=255)
                                                                                        aaGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCgg                                                                                           <  1:868755/113‑1 (MQ=255)
                                                                                            aaTCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGc                                                   <  4:592167/149‑1 (MQ=255)
                                                                                                     cGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGt                                           >  2:279354/1‑148 (MQ=255)
                                                                                                     cGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGt                                           >  2:299188/1‑148 (MQ=255)
                                                                                                                  cacCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCg                              <  1:655169/148‑1 (MQ=255)
                                                                                                                   acCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCggg                            <  3:598707/149‑1 (MQ=255)
                                                                                                                   acCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCggg                            >  3:863395/1‑149 (MQ=255)
                                                                                                                    ccTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCtgg                            >  2:291975/1‑148 (MQ=255)
                                                                                                                     cTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGat                          >  3:35874/1‑149 (MQ=255)
                                                                                                                       gTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGatat                        >  3:191500/1‑149 (MQ=255)
                                                                                                                               ggcGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGtt                                       <  1:173310/126‑1 (MQ=255)
                                                                                                                               ggcGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGtt                                       >  2:173310/1‑126 (MQ=255)
                                                                                                                               ggcGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCtgtg                 >  2:573620/1‑148 (MQ=255)
                                                                                                                                gcGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGa               <  4:876310/149‑1 (MQ=255)
                                                                                                                                   tttCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATgg            >  4:218635/1‑149 (MQ=255)
                                                                                                                                    ttCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATgg            <  2:6247/148‑1 (MQ=255)
                                                                                                                                    ttCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATgg            <  1:279354/148‑1 (MQ=255)
                                                                                                                                      ccTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCg          <  2:385128/148‑1 (MQ=255)
                                                                                                                                      ccTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGt         <  4:863395/149‑1 (MQ=255)
                                                                                                                                       taagcgccccTGGCAGGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGt         >  4:813332/7‑149 (MQ=255)
                                                                                                                                        tGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGa       >  4:798579/1‑149 (MQ=255)
                                                                                                                                         gcgcCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGAc      <  3:605424/149‑1 (MQ=255)
                                                                                                                                          gcgccccTGGCAGGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGa       >  2:102901/4‑148 (MQ=255)
                                                                                                                                             ccccTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATg                                             >  4:31568/1‑106 (MQ=255)
                                                                                                                                             ccccTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATg                                             <  3:31568/106‑1 (MQ=255)
                                                                                                                                             ccccTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACtgcc  >  3:289925/1‑149 (MQ=255)
                                                                                                                                             ccccTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACtgc   <  1:299188/148‑1 (MQ=255)
                                                                                                                                              cccTGGCAGGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCg                              >  2:518456/1‑121 (MQ=255)
                                                                                                                                                   |                                                                                                                                               
CAAAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCA‑GGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACTGCC  >  NC_000913/790784‑791073

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