Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 790,932 C→A 6.3% W33L (TGG→TTG) ‡ galT ← galactose‑1‑phosphate uridylyltransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913790,9320CA 6.3% 303.9 / 9.0 95W33L (TGG→TTG) ‡galTgalactose‑1‑phosphate uridylyltransferase
Reads supporting (aligned to +/- strand):  ref base C (39/50);  new base A (2/4);  total (41/54)
Fisher's exact test for biased strand distribution p-value = 6.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.25e-02

AAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACTGCCAGTGC  >  NC_000913/790786‑791078
                                                                                                                                                  |                                                                                                                                                  
aaGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAg                                                                                                                                                  >  3:360332/1‑149 (MQ=255)
 aGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGGGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAgg                                                                                                                                                 <  4:426517/149‑1 (MQ=255)
  gCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAggg                                                                                                                                                >  1:386046/1‑149 (MQ=255)
   ccGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAgggg                                                                                                                                               <  5:14510/149‑1 (MQ=255)
   ccGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAgggg                                                                                                                                               <  6:59864/149‑1 (MQ=255)
   ccGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAgggg                                                                                                                                               <  3:202567/149‑1 (MQ=255)
    cGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGgc                                                                                                                                              >  5:315921/1‑149 (MQ=255)
       aaaGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCTGGCAGGGGgcgca                                                                                                                                           <  8:165598/149‑2 (MQ=255)
        aaGGCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCTGGCAGGGGgcgcag                                                                                                                                          <  3:148665/149‑3 (MQ=255)
         aGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTa                                                                                                                                         <  6:315921/149‑1 (MQ=255)
             aTTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGCGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAgcgc                                                                                                                                     >  3:76417/1‑149 (MQ=255)
             aTTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCTGGCAGGGGgcgca                                                                                                                                           <  7:247990/143‑2 (MQ=255)
             aTTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCTGGCAGGGGgcgca                                                                                                                                           <  7:241399/143‑2 (MQ=255)
             aTTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCTGGCAGGGGgcgca                                                                                                                                           >  8:247990/1‑142 (MQ=255)
             aTTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCTGGCAGGGGgcgca                                                                                                                                           >  8:241399/1‑142 (MQ=255)
             aTTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAgcgc                                                                                                                                     >  3:13120/1‑149 (MQ=255)
                 gTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAg                                                                                                                                                  <  8:376192/132‑1 (MQ=255)
                 gTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAg                                                                                                                                                  >  7:376192/1‑132 (MQ=255)
                  tGAAAACGTAAGTCCCGGTGTATTCGGGGTTTTTATCGCCTGTCACCCTCACATTACCTGCGCTGAGGATGCAATCTGGAGCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCgcggtgg                                                                                                                                <  8:92056/149‑2 (MQ=255)
                    aaaaCGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCgcggtgcgg                                                                                                                              >  5:373637/1‑149 (MQ=255)
                     aaaCGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCgcggtgcggt                                                                                                                             <  8:294717/149‑1 (MQ=255)
                        cGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGa                                                                                                                           >  3:304847/1‑148 (MQ=255)
                         gTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGaaa                                                                                                                         <  8:328306/149‑1 (MQ=255)
                           aaGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAAc                                                                                                                        >  5:316670/1‑148 (MQ=255)
                                ccGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGTAATCTGGATCGTGCGCAAGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAAt                                                                                                                  >  6:319853/1‑149 (MQ=255)
                                ccGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGaa                                                                                                                   <  3:10612/148‑1 (MQ=255)
                                ccGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGaa                                                                                                                   <  5:170878/148‑1 (MQ=255)
                                ccGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAAt                                                                                                                  <  3:421138/149‑1 (MQ=255)
                                    tgtAAGCGCGGTGTATCTGGCCTGTCTCGCGCTCATTACCAGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCGTGCCCCCCGAGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATcctt                                                                                                              <  7:446011/149‑3 (MQ=255)
                                       aaTCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGc                                                                                                           >  8:200951/1‑149 (MQ=255)
                                            gggTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGt                                                                                                      >  3:123927/1‑149 (MQ=255)
                                               tttttATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAg                                                                                                   >  2:318856/1‑149 (MQ=255)
                                                 tttATCGCCTGTCACCCGCTCATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGc                                                                                                  <  5:425490/148‑1 (MQ=255)
                                                 tttATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGc                                                                                                  <  6:316670/148‑1 (MQ=255)
                                                  ttATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAAGACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCg                                                                                                 <  3:252441/148‑1 (MQ=255)
                                                    aTCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGGTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCggg                                                                                               <  1:236872/148‑1 (MQ=255)
                                                     tCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGtt                                                                                             <  1:318856/149‑1 (MQ=255)
                                                      cGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGCATCCATTGCCCGGTGAGCGGGTTg                                                                                            >  5:19316/1‑149 (MQ=255)
                                                                ccGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGTGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGAt                                                                                  <  5:26016/149‑1 (MQ=255)
                                                                ccGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGAt                                                                                  <  4:364202/149‑1 (MQ=255)
                                                                  gctcaTTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGAtgt                                                                                <  4:45378/146‑1 (MQ=255)
                                                                  gCACATTACCTGCGCAGAGGAAGTAATCTGGATCGTGCGCAAGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGAtgt                                                                                <  5:319853/149‑1 (MQ=255)
                                                                  gCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCGTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGGGCGCGGGTTGTAGCGGCGAggt                                                                                >  7:364290/1‑149 (MQ=255)
                                                                            tGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGga                                                                             >  8:264034/1‑142 (MQ=255)
                                                                            tGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGga                                                                             <  7:264034/142‑1 (MQ=255)
                                                                              cgcAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAAc                                                                     <  6:67633/148‑1 (MQ=255)
                                                                              cgcAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACg                                                                    <  4:13120/149‑1 (MQ=255)
                                                                                 agagGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGa                                                                 <  5:447265/149‑1 (MQ=255)
                                                                                 agagGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGa                                                                 <  7:35800/149‑1 (MQ=255)
                                                                                    ggAAGCAATCTGGATCGTGCGCAGGTACCACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGAta                                                               >  4:375556/1‑147 (MQ=255)
                                                                                     gAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTa                                                              >  4:211051/1‑148 (MQ=255)
                                                                                     gAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCCGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATc                                                                        >  1:282006/1‑138 (MQ=255)
                                                                                     gAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCCGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGACCGGGTTGTAGCGGCGATGTGGATGATc                                                                        <  2:282006/138‑1 (MQ=255)
                                                                                       aGCAATCTGGAGCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCGCTGCCAGTGGCGCTTAGCGCGGTGCGGTGTAACCAGAATCCATTGCCCGGTGAGCggg                                                                                               <  5:186030/113‑1 (MQ=255)
                                                                                        gCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGcggtgcggtg                                                                                                                            >  7:218720/1‑83 (MQ=255)
                                                                                        gCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGcggtgcggtg                                                                                                                            <  8:218720/83‑1 (MQ=255)
                                                                                        gCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGa                                                                                                    >  2:244627/1‑107 (MQ=255)
                                                                                        gCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGa                                                                                                    <  1:244627/107‑1 (MQ=255)
                                                                                          aaTCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGc                                                        <  1:176222/149‑1 (MQ=255)
                                                                                          aaTCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGc                                                        >  6:413954/1‑149 (MQ=255)
                                                                                          aaTCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTGCCTGCTCCCCCTGCCAGGGGCGCTTATCGCGGTGCGGTGTAACCATAATCGATTGCCGTGTGAGCGGGTTGTAGCGGCGTTGTTGATGATCTACGGGATTATATTGc                                                        <  8:373811/149‑1 (MQ=255)
                                                                                           aTCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCg                                                       >  5:267719/1‑149 (MQ=255)
                                                                                           aTCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCg                                                       <  5:413954/149‑1 (MQ=255)
                                                                                             cTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAgcggc                                                                                     >  4:53112/1‑117 (MQ=255)
                                                                                             cTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAgcggc                                                                                     <  3:53112/117‑1 (MQ=255)
                                                                                              tGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCa                                                    >  8:142861/1‑149 (MQ=255)
                                                                                                gATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATg                                                  <  7:142861/149‑1 (MQ=255)
                                                                                                gATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGGTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGt                                                      <  4:76417/145‑1 (MQ=255)
                                                                                                gATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGAGCGGTGAAACCAGAATCCATTGCCCGGTGAGCCGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCa                                                    >  3:222849/1‑147 (MQ=255)
                                                                                                  tcctGCGCAGGTAACACCTGATTGGCTGGCGTTTCCTGCGCCCCCAGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGt                                                >  5:297487/4‑149 (MQ=255)
                                                                                                  tCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGt                                                <  4:123927/149‑1 (MQ=255)
                                                                                                     tGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTgg                                                                              >  6:23099/1‑116 (MQ=255)
                                                                                                     tGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGt                                             <  2:386046/149‑1 (MQ=255)
                                                                                                       cgcAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGtt                                            <  6:373637/148‑1 (MQ=255)
                                                                                                       cgcAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGt                                             >  4:281864/1‑147 (MQ=255)
                                                                                                       cgcAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGt                                             <  3:281864/147‑1 (MQ=255)
                                                                                                        gcAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTcc                                          <  7:385396/149‑1 (MQ=255)
                                                                                                               acacCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAATCAATAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTGCTTAATc                                    <  5:73501/148‑1 (MQ=255)
                                                                                                                cacCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTGGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCgg                                  <  4:360332/149‑1 (MQ=255)
                                                                                                                 cccTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCgg                                  >  8:418012/2‑148 (MQ=255)
                                                                                                                 acCTGTGAGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCAGTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCACCGGGATTAAATTGCGGCATGGTCGTTCCTTAATCgg                                  >  6:380640/1‑148 (MQ=255)
                                                                                                                   cTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGat                               <  7:200951/149‑1 (MQ=255)
                                                                                                                   cTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGAGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGa                                >  8:454011/1‑148 (MQ=255)
                                                                                                                     gtgtggctggcGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCg                                                                                                        <  6:312233/71‑1 (MQ=255)
                                                                                                                            tggcGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCGACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCtgtg                      >  3:393531/1‑149 (MQ=255)
                                                                                                                            tggcGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCtgtg                      >  2:337906/1‑149 (MQ=255)
                                                                                                                            tggcGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCtgtg                      >  3:26381/1‑149 (MQ=255)
                                                                                                                              gggTTGCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCGATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGTGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTgg                     <  6:297487/146‑1 (MQ=255)
                                                                                                                                gTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGACCCCAGACTCCATTGCCCGGTGCGCGGGTTGTAGCGGCGATGTGGATGATCACCGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATg                  >  5:147461/1‑149 (MQ=255)
                                                                                                                                    ccTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGt              >  7:336090/1‑149 (MQ=255)
                                                                                                                                     cTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGt              <  3:211051/148‑1 (MQ=255)
                                                                                                                                      tGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGa            >  2:154899/1‑149 (MQ=255)
                                                                                                                                          cccccTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACtgc        <  2:101052/149‑1 (MQ=255)
                                                                                                                                            cccTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGGTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTACTCGGGATATCCCTGTGGATGGCGTGACtgcca      <  6:147461/149‑1 (MQ=255)
                                                                                                                                                gCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTAATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACtgccagtgc  >  4:160919/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
AAGCCGCAAAGTCATTAGTGAAAACGTAAGTCCCGGTGTAATCGGGGTTTTTATCGCCTGTCACCCGCACATTACCTGCGCAGAGGAAGCAATCTGGATCGTGCGCAGGTAACACCTGTTTGGCTGGCGTTTCCTGCGCCCCCTGCCAGGGGCGCTTAGCGCGGTGCGGTGAAACCAGAATCCATTGCCCGGTGAGCGGGTTGTAGCGGCGATGTGGATGATCAACGGGATTAAATTGCGTCATGGTCGTTCCTTAATCGGGATATCCCTGTGGATGGCGTGACTGCCAGTGC  >  NC_000913/790786‑791078

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: