Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 265,418 A→G 55.9% R42R (CGT→CGC ykfF ← CP4‑6 prophage; uncharacterized protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913265,4180AG55.9% 62.6 / 150.7 132R42R (CGT→CGCykfFCP4‑6 prophage; uncharacterized protein
Reads supporting (aligned to +/- strand):  ref base A (32/26);  new base G (42/31);  total (75/57)
Fisher's exact test for biased strand distribution p-value = 8.60e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 135 total aligned reads displayed.

TCTGTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCC  >  NC_000913/265320‑265516
                                                                                                  |                                                                                                  
tCTGTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCAc                                                                                                   <  2:2345811/100‑1 (MQ=255)
 cTGTACGGATGCCTGATGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGACTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATcgcg                                                                                                  <  1:1976784/100‑1 (MQ=255)
  tGTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGa                                                                                                 >  1:211317/1‑100 (MQ=255)
   gTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGaa                                                                                                >  2:2408734/1‑100 (MQ=255)
   gTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGaa                                                                                                >  1:141699/1‑100 (MQ=255)
      cGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAacca                                                                                             >  2:1392275/1‑100 (MQ=255)
            ccTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGAc                                                                                       >  2:3306993/1‑100 (MQ=255)
            ccTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGAc                                                                                       <  2:3096214/100‑1 (MQ=255)
             cTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACg                                                                                      >  2:3119091/1‑100 (MQ=255)
             cTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACg                                                                                      >  1:3164275/1‑100 (MQ=255)
             cTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACg                                                                                      >  1:2658929/1‑100 (MQ=255)
               gAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGGCCATCCGGCCTTCAGTATCACGAACCACCAGACGGa                                                                                    >  1:203849/1‑100 (MQ=255)
               gAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGa                                                                                    >  2:800572/1‑100 (MQ=255)
                aGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGaa                                                                                   >  2:2473099/1‑100 (MQ=255)
                 ggTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAg                                                                                  >  2:2712324/1‑100 (MQ=255)
                  gTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGt                                                                                 >  2:784624/1‑100 (MQ=255)
                  gTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGt                                                                                 <  1:800569/100‑1 (MQ=255)
                   tGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTg                                                                                >  2:265425/1‑100 (MQ=255)
                   tGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTg                                                                                <  1:799996/100‑1 (MQ=255)
                    gCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGa                                                                               <  2:1740105/100‑1 (MQ=255)
                    gCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGa                                                                               >  2:3327802/1‑100 (MQ=255)
                     cGGATATAGCGGTTCAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCGGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGAc                                                                              <  1:3242812/100‑1 (MQ=255)
                     cGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGAc                                                                              <  1:2566796/100‑1 (MQ=255)
                     cGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGAc                                                                              <  2:2881373/100‑1 (MQ=255)
                       gATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTg                                                                            <  2:2770720/100‑1 (MQ=255)
                         ttccgatcttAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGcc                                                                          <  1:2217227/92‑1 (MQ=255)
                         tataGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGc                                                                           >  2:2644545/1‑99 (MQ=255)
                          ataGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGccc                                                                         <  2:2084735/100‑1 (MQ=255)
                           taGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCt                                                                        >  2:2985127/1‑100 (MQ=255)
                             gCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTgg                                                                      >  2:3027405/1‑100 (MQ=255)
                             gCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTgg                                                                      <  1:1453161/100‑1 (MQ=255)
                                gTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGgtcg                                                                   >  1:2125927/1‑100 (MQ=255)
                                gTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGgtcg                                                                   >  2:1529704/1‑100 (MQ=255)
                                gTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGgtcg                                                                   <  1:2507157/100‑1 (MQ=255)
                                  tAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGccagatcggaa                                                                 >  2:2217242/1‑91 (MQ=255)
                                   aggACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCt                                                                >  2:2301180/3‑100 (MQ=255)
                                   aaGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACGATCCGGCCTTCAGTATCGCGAACCAGCAGACGGAGGGGACTGCCCTGgtcgt                                                                  >  1:625159/1‑98 (MQ=21)
                                     gACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTc                                                              <  2:2931813/100‑1 (MQ=255)
                                     gACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTc                                                              >  2:125387/1‑100 (MQ=255)
                                      aCCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCg                                                             <  2:3023773/100‑1 (MQ=255)
                                      aCCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCg                                                             >  2:530831/1‑100 (MQ=255)
                                       ccTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGa                                                            >  1:1124338/1‑100 (MQ=255)
                                       ccTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGa                                                            >  1:75207/1‑100 (MQ=255)
                                       ccTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGa                                                            <  1:3326416/100‑1 (MQ=255)
                                        cTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAg                                                           <  2:518390/100‑1 (MQ=255)
                                        cTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAg                                                           <  2:63901/100‑1 (MQ=255)
                                        cTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAg                                                           <  1:1741520/100‑1 (MQ=255)
                                         ttCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAgt                                                          >  1:15282/1‑100 (MQ=255)
                                          tCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAgtg                                                         <  1:3192662/100‑1 (MQ=255)
                                           cACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAgtgt                                                        <  1:1480978/100‑1 (MQ=255)
                                               ggCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCCCCAGACGGCAGTGACTGCCCTGGTCGTCTTCGAGTGTGata                                                    >  1:1588711/1‑100 (MQ=255)
                                               ggCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGata                                                    >  1:389832/1‑100 (MQ=255)
                                               ggCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGata                                                    <  2:1398624/100‑1 (MQ=255)
                                               ggCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGGCTTCGAGTGTGat                                                     >  1:1605153/1‑99 (MQ=255)
                                               gacATCCGGCTCATAGTTCCATGCCCGCCCGACCATCCGGCCTTCAGTATCACGAACCACCATCCGGAAGTGTGTGCCCTGGTCGTCTTCGAGtgtg                                                       >  1:1167928/3‑97 (MQ=25)
                                                gCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGatat                                                   >  1:1872450/1‑100 (MQ=255)
                                                    ccGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGc                                                >  1:914527/1‑99 (MQ=255)
                                                     cGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTg                                              <  2:1704836/100‑1 (MQ=255)
                                                     cGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTg                                              >  2:1344108/1‑100 (MQ=255)
                                                       gCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGt                                             <  2:1311982/99‑1 (MQ=255)
                                                        cTCAAAGTCCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGCAGTGCCTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTa                                            >  1:1094624/1‑99 (MQ=37)
                                                          cAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGt                                         >  1:349659/1‑100 (MQ=255)
                                                            aaGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTgg                                       <  1:120507/100‑1 (MQ=255)
                                                             aGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGt                                      <  2:455574/100‑1 (MQ=255)
                                                             aGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCGTCGAGTGTGATATTGCTGTACGTgg                                       >  1:112479/1‑99 (MQ=255)
                                                              gTTCCATGCCCGTCAGACCAGCCGGCCTTCAGTATCGCGAACCACCAGTCGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTgg                                       <  2:1670103/98‑1 (MQ=37)
                                                                tCCATGCCCGCCAGACCAGCCGGCCTTCAGGATCGCGAACCACGAGGCGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAg                                    <  1:1165683/99‑1 (MQ=37)
                                                                 ccATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTg                                  >  2:1234350/1‑100 (MQ=255)
                                                                 ccATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTg                                  <  1:2314499/100‑1 (MQ=255)
                                                                  cATGCCCGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACGGCCCTGGTCGTCTTCGAGTGTGCTATTGCTGTACGTGGTAGTg                                  >  1:870946/1‑99 (MQ=255)
                                                                  cATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGa                                 >  2:2132792/1‑100 (MQ=255)
                                                                  aaTGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGa                                 >  2:212866/2‑100 (MQ=255)
                                                                      cccGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcgc                             >  2:2617461/1‑100 (MQ=255)
                                                                        cGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcgctt                           >  1:2508867/1‑100 (MQ=255)
                                                                        cGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcgctt                           >  2:3266828/1‑100 (MQ=255)
                                                                        cGCCAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcgct                            >  2:3149974/1‑99 (MQ=255)
                                                                         gCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcgcttg                          >  2:2561749/1‑100 (MQ=255)
                                                                          ccAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcgcttgc                         <  2:2939134/100‑1 (MQ=255)
                                                                           cAGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcgcttgcg                        <  2:2222379/100‑1 (MQ=255)
                                                                            aGACCATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGCTAGTGACcgcttgcgc                       <  2:2369577/100‑1 (MQ=255)
                                                                             gACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcgcttgcgct                      <  1:2737041/100‑1 (MQ=255)
                                                                               ccATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcgcttgcgcttg                    <  2:3206028/100‑1 (MQ=255)
                                                                               ccATCCGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACcgcttgcgcttg                    <  2:2486287/100‑1 (MQ=255)
                                                                                cATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGt                   >  1:3140222/1‑100 (MQ=255)
                                                                                 aTCCGGCCTTCAGTATCTCGCACCACCGGGCTGATTTGACTACCCTTTTCGTCTTCCAGTGTGCTATTGCTGTTCTTCGTACTGACcccttcagc                       >  1:374822/1‑91 (MQ=255)
                                                                                   ccGGCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGtctc                >  1:852106/1‑100 (MQ=255)
                                                                                     ggCCTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCg              <  1:885809/100‑1 (MQ=255)
                                                                                       ccTTCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCggg            >  2:1577581/1‑100 (MQ=255)
                                                                                       ccTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCggg            <  1:730179/100‑1 (MQ=255)
                                                                                        cTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGt           <  2:1523481/100‑1 (MQ=255)
                                                                                          tCAGTATCGCGATCCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCggg            <  1:2750005/97‑1 (MQ=255)
                                                                                          tCAGTATCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGa         >  2:428657/1‑100 (MQ=255)
                                                                                          tCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGa         <  2:1897095/100‑1 (MQ=255)
                                                                                               aTCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGc    >  1:543616/1‑100 (MQ=255)
                                                                                               aTCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGc    >  1:310227/1‑100 (MQ=255)
                                                                                               aTCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGc    >  2:2715580/1‑100 (MQ=255)
                                                                                               aTCGCGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGc    >  2:2650983/1‑100 (MQ=255)
                                                                                                 cACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGccc  <  1:530828/100‑1 (MQ=255)
                                                                                                 cACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGccc  >  2:1923573/1‑100 (MQ=255)
                                                                                                 cACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTACTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGccc  >  2:35069/1‑100 (MQ=255)
                                                                                                  |                                                                                                  
TCTGTACGGATGCCTGAGGTGCGGATATAGCGGTTAAGACCTTCACCGGCATCCGGCTCAAAGTTCCATGCCCGCCAGACCATCCGGCCTTCAGTATCACGAACCACCAGACGGAAGTGACTGCCCTGGTCGTCTTCGAGTGTGATATTGCTGTACGTGGTAGTGACCGCTTGCGCTTGTCTCCGGGTGAAAGGCCC  >  NC_000913/265320‑265516

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