Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 422,062 G→A 5.4% V359V (GTG→GTA proY → proline‑specific permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913422,0620GA 5.4% 247.2 / 3.1 74V359V (GTG→GTAproYproline‑specific permease
Reads supporting (aligned to +/- strand):  ref base G (34/36);  new base A (2/2);  total (36/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.82e-01

GTATGGCAGAGCAGGGCAGCGCGCCGAAAATTTTCAGCAAAACGTCGCGTCGCGGTATTCCGTGGGTTACGGTGCTGGTGATGACTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTA  >  NC_000913/421917‑422191
                                                                                                                                                 |                                                                                                                                 
gTATGGTAGAGCAGGGCAGCGCGCCGAAAATTTTCAGCAATCCGTCGCGTCGCGGTATTCCGTGGGTTACGGTGCTGGTGATGACTACCGCGCTGCTGTATGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGAt                                                                                                                                 <  8:239283/148‑1 (MQ=255)
gTATGGTAGAGCAGGGCAGCGCGCCGAAAATTTTCAGCAAAACCTCGCGTCGCGGTATTCCGTGGGTTACGGTGCTGGAGATGACTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGAt                                                                                                                                 <  8:245462/148‑1 (MQ=255)
        gcgCAGGGCAGCGCGCCGAAAATTTTCAGCAAAACGTCGCGTCGCGGTATTCCGTGGGTTACGGTGCTGGTGATGACTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGAtcgcttcgc                                                                                                                         >  3:192670/3‑148 (MQ=255)
          gCAGGGCAGCGCGCCGAAAATTTTCAGCAAAACGTCGCGTCGCGGTATTCCGTGGGTTACGGTGCTGGTGATGACTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTg                                                                                                                       <  7:54723/148‑1 (MQ=255)
               gCAGCGCGCCGAAAATTTTCAGCAAAACGTCGCGTCGCGGTATTCCGTGGGTTACGGTGCTGGTGATGACTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGAtcgcttcg                                                                                                                          <  4:288272/140‑1 (MQ=255)
               gCAGCGCGCCGAAAATTTTCAGCAAAACGTCGCGTCGCGGTATTCCGTGGGTTACGGTGCTGGTGATGACTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGAtcgcttcg                                                                                                                          >  3:288272/1‑140 (MQ=255)
               gCAGCGCGCCGAAAATTTTCAGCAAAACGTCGCGTCGCGGTATTCCGTGGGTTACGGTGCTGGTGATGACTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAAc                                                                                                                  >  5:249659/1‑148 (MQ=255)
                   cgcgcCGAAAATTTTCAGCAAAACGTCGCGTCGCGGTATTCCGTGGGTTACGGTGCTGGTGATGACTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTc                                                                                                              <  8:66982/148‑1 (MQ=255)
                   cgcgcCGAAAATTTTCAGCAAAACGTCGCGTCGCGGTATTCCGTGGGTTACGGTGCTGGTGATGACTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTc                                                                                                              >  2:378458/1‑148 (MQ=255)
                      gcCGAAAATTTTCAGCAAAACGTCGCGTCGCGGTATTCCGTGGGTTACGGTGCTGGTGATGACTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGcc                                                                                                           <  5:279505/148‑1 (MQ=255)
                           aaaTTTTCAGCAAAACGTCGCGTCGCGGTATTCCGTGGGTTACGGTGCTGGTGATGACTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACggt                                                                                                      <  6:372921/148‑1 (MQ=255)
                           aaaTTTTCAGCAAAACGTCGCGTCGCGGTATTCCGTGGGTTACGGTGCTGGTGATGACTACCGCGCTGCGGTTTGCGGTGTATCTGCACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACggt                                                                                                      >  3:197707/1‑148 (MQ=255)
                           aaaTTTCCAGCAAAACGTCGCGTCGCGGTATTCCGGGGGTTACGGTGCTGGTGATGACTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAACACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACggt                                                                                                      >  4:337835/1‑148 (MQ=255)
                                  cAGCAAAACGTCGCGTCGCGGTATTCCGTGGGTTACGGTGCTGGTGATGACTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACggtgtgggt                                                                                                <  6:249659/147‑1 (MQ=255)
                                     cAAAACGTCGCGTCGCGGTATTCCGTGGGTTACGGTGCTGGTGATGACTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACggt                                                                                                      >  7:324549/1‑138 (MQ=255)
                                     cAAAACGTCGCGTCGCGGTATTCCGTGGGTTACGGTGCTGGTGATGACTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACggt                                                                                                      <  8:324549/138‑1 (MQ=255)
                                                gtcgcgGTATTCCGTGGGTTACGGTGCTGGTGATGACTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCct                                                                                 >  7:103956/1‑148 (MQ=255)
                                                 tcgcgGTATTCCGTGGGTTACGGTGCTGGTGATGACTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCctg                                                                                <  1:378458/148‑1 (MQ=255)
                                                   gcgGTATTCCGTGGGTTACGGTGCTGGGGATGACTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCctgc                                                                               <  7:108090/147‑1 (MQ=255)
                                                                  ttACGGTGCTGGTGATGACTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGAt                                                                                    <  4:18990/127‑1 (MQ=255)
                                                                  ttACGGTGCTGGTGATGACTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGAt                                                                                    >  3:18990/1‑127 (MQ=255)
                                                                    aCGGTGCTGGTGATGACTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTcc                                                             >  3:21124/1‑148 (MQ=255)
                                                                            ggTGATGACTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGACACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGt                                                     >  4:181880/1‑148 (MQ=255)
                                                                                atgaCTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTg                                                  <  8:103956/147‑1 (MQ=255)
                                                                                  gaCTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTg                                                                   >  6:168965/1‑128 (MQ=255)
                                                                                  gaCTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTg                                                                   <  5:168965/128‑1 (MQ=255)
                                                                                  gaCTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTgccg                                               >  3:218247/1‑148 (MQ=255)
                                                                                  gaCTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGACTATGATCCTGCTGTCGCaaa                                                                      >  2:256860/1‑125 (MQ=255)
                                                                                  gaCTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGACTATGATCCTGCTGTCGCaaa                                                                      <  1:256860/125‑1 (MQ=255)
                                                                                   aCTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTgccgc                                              >  7:235934/1‑148 (MQ=255)
                                                                                      aCCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTgccgc                                              >  8:25506/1‑145 (MQ=255)
                                                                                      aCCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTgc                                                 <  8:344157/142‑1 (MQ=255)
                                                                                      aCCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCAGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTgc                                                 >  7:344157/1‑142 (MQ=255)
                                                                                       ccGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCag                                           >  2:25529/1‑147 (MQ=255)
                                                                                              gctgTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCtt                                                               >  8:68517/1‑120 (MQ=255)
                                                                                              gctgTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCtt                                                               <  7:68517/120‑1 (MQ=255)
                                                                                               ctgTTTGCGGTGTATCTGAACTCCATCATGCCGGAAACCGGCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTa                                  <  4:197707/148‑1 (MQ=255)
                                                                                               ctgTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTa                                  <  4:192670/148‑1 (MQ=255)
                                                                                                 gTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAg                                >  3:30868/1‑148 (MQ=255)
                                                                                                   ttGCGGTGTATCTGAAATACATCATGCCGGAAAACGTCTTCCTGTTGATCGCTTCGCTGGAAACCTTCGCCACGGTGTGGGTGTGGATTATTATc                                                                                   >  2:269723/1‑95 (MQ=255)
                                                                                                    tGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGgcg                             >  5:329863/1‑148 (MQ=255)
                                                                                                             tCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCaaa                                                                      <  7:298184/98‑1 (MQ=255)
                                                                                                             tCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCaaa                                                                      >  8:298184/1‑98 (MQ=255)
                                                                                                                 aaCTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAg                                <  1:361208/132‑1 (MQ=255)
                                                                                                                 aaCTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAg                                >  2:361208/1‑132 (MQ=255)
                                                                                                                 aaCTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTaa                  <  4:218247/146‑1 (MQ=255)
                                                                                                                 aaCTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTaa                  <  3:181880/146‑1 (MQ=255)
                                                                                                                 aaCTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGGTGATTATGATCCTGCTGGCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTa                   <  4:30868/145‑1 (MQ=255)
                                                                                                                  aCTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGt               <  4:21124/148‑1 (MQ=255)
                                                                                                                     aCATCATGCCGGAACACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGcc                                                       >  5:176590/1‑105 (MQ=255)
                                                                                                                     aCATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGcc                                                       <  6:176590/105‑1 (MQ=255)
                                                                                                                     aCATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTgc                                                 <  2:353214/111‑1 (MQ=255)
                                                                                                                     aCATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTgc                                                 >  1:353214/1‑111 (MQ=255)
                                                                                                                     aCATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTaa                  <  1:257683/142‑1 (MQ=255)
                                                                                                                     aCATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTaa                  >  2:257683/1‑142 (MQ=255)
                                                                                                                     aCATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGcc            >  7:136287/1‑148 (MQ=255)
                                                                                                                     aCATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGcc            >  3:336404/1‑148 (MQ=255)
                                                                                                                     aCATCATGCCGGAAAACGTCTTCCTGGTAATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTaa                  >  6:265305/1‑142 (MQ=255)
                                                                                                                     aCATCATGCCGGAAAACGTCTTCCTGGTAATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTaa                  >  6:265303/1‑142 (MQ=255)
                                                                                                                     aCATCATGCCGGAAAACGTCTTCCTGGTAATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTaa                  <  5:265305/142‑1 (MQ=255)
                                                                                                                     aCATCATGCCGGAAAACGTCTTCCTGGTAATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTaa                  <  5:265303/142‑1 (MQ=255)
                                                                                                                      catcatGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCg           <  8:136287/148‑1 (MQ=255)
                                                                                                                       atcatGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCgg          <  7:25506/148‑1 (MQ=255)
                                                                                                                        tcatGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCggg         >  6:367094/1‑148 (MQ=255)
                                                                                                                        tcatGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCggg         <  4:354069/148‑1 (MQ=255)
                                                                                                                        tcatGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTAAAATTTAAAGTGCCggg         >  5:326685/1‑148 (MQ=255)
                                                                                                                         catGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGaagaag                                     >  6:151951/1‑119 (MQ=255)
                                                                                                                         catGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGaagaag                                     <  5:151951/119‑1 (MQ=255)
                                                                                                                          atGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTg       <  8:205846/148‑1 (MQ=255)
                                                                                                                               ggAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTa  <  8:235934/148‑1 (MQ=255)
                                                                                                                                             tGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTa  >  2:280290/1‑134 (MQ=255)
                                                                                                                                             tGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTa  <  1:280290/134‑1 (MQ=255)
                                                                                                                                              ggTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTa                                  >  8:399051/1‑101 (MQ=255)
                                                                                                                                              ggTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTa                                  <  7:399051/101‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                 
GTATGGCAGAGCAGGGCAGCGCGCCGAAAATTTTCAGCAAAACGTCGCGTCGCGGTATTCCGTGGGTTACGGTGCTGGTGATGACTACCGCGCTGCTGTTTGCGGTGTATCTGAACTACATCATGCCGGAAAACGTCTTCCTGGTGATCGCTTCGCTGGCAACCTTCGCCACGGTGTGGGTGTGGATTATGATCCTGCTGTCGCAAATTGCCTTCCGTCGCCGTTTGCCGCCAGAAGAAGTTAAGGCGCTGAAATTTAAAGTGCCGGGTGGGGTA  >  NC_000913/421917‑422191

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