Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 338,761 G→A V405I (GTC→ATC)  yahG → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913338,7610GA85.2% 724.5 / 82.1 265V405I (GTC→ATC) yahGhypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (19/20);  new base A (128/97);  total (147/118)
Fisher's exact test for biased strand distribution p-value = 3.85e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.42e-01

Only 100 of 265 total aligned reads displayed.

GCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACA  >  NC_000913/338671‑338851
                                                                                          |                                                                                          
gCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATc                                                                                          >  1:6849248/1‑93 (MQ=255)
gCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATc                                                                                          >  1:5916739/1‑93 (MQ=255)
gCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATc                                                                                          >  2:5884877/1‑93 (MQ=255)
  tATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCAc                                                                                        <  1:4884214/93‑1 (MQ=255)
  tATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCAc                                                                                        <  1:2329982/93‑1 (MQ=255)
  tATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCAc                                                                                        >  2:6174414/1‑93 (MQ=255)
   aTCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCAcc                                                                                       >  2:3136622/1‑93 (MQ=255)
      gtcgCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCAtt                                                                                    >  2:1108441/1‑93 (MQ=255)
       tcgCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTc                                                                                   <  1:3593604/93‑1 (MQ=255)
       tcgCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTc                                                                                   <  2:2597770/93‑1 (MQ=255)
       tcgCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTc                                                                                   <  1:406750/93‑1 (MQ=255)
       tcgCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTc                                                                                   <  1:6455175/93‑1 (MQ=255)
        cgcgCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTcc                                                                                  <  2:7222293/93‑1 (MQ=255)
            cTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCgctg                                                                              <  1:3858694/93‑1 (MQ=255)
            cTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCgctg                                                                              >  1:3976266/1‑93 (MQ=255)
             tggtggGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCgctgc                                                                             >  1:2937022/1‑93 (MQ=255)
               gtggGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTc                                                                           >  2:7127458/1‑93 (MQ=255)
               gtggGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTc                                                                           >  1:2070459/1‑93 (MQ=255)
                tggGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCg                                                                          >  1:3550162/1‑93 (MQ=255)
                tggGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCg                                                                          <  2:4025020/93‑1 (MQ=255)
                 gggCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCgg                                                                         <  2:6085720/93‑1 (MQ=255)
                  ggcggcACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGt                                                                        >  2:3542656/1‑93 (MQ=255)
                  ggcggcACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGt                                                                        >  2:1521844/1‑93 (MQ=255)
                     ggcACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGtttt                                                                     >  2:5124306/1‑93 (MQ=255)
                     ggcACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGtttt                                                                     <  1:7317846/93‑1 (MQ=255)
                      gcACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTa                                                                    <  2:2213716/93‑1 (MQ=255)
                      gcACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTa                                                                    <  1:5582878/93‑1 (MQ=255)
                      gcACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGCGAACACCCCAACATCACCATTCCGCTGCTCGGTTTTa                                                                    >  2:1954072/1‑93 (MQ=255)
                        aCGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATg                                                                  <  1:3187323/93‑1 (MQ=255)
                        aCGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATg                                                                  <  2:644018/93‑1 (MQ=255)
                            tGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCCCCATTCCGCTGCTCGGTTTTATGGGCg                                                              <  1:4339235/93‑1 (MQ=255)
                             ggAAGAAGCTATTGATTTCTCACGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGt                                                             <  1:2689393/93‑1 (MQ=255)
                               aagaagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGc                                                           >  1:3304596/1‑93 (MQ=255)
                                agaagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGcc                                                          >  2:3338639/1‑93 (MQ=255)
                                 gaagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCg                                                         >  1:4560421/1‑93 (MQ=255)
                                 gaagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCg                                                         >  2:6413967/1‑93 (MQ=255)
                                  aagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGt                                                        <  2:459343/93‑1 (MQ=255)
                                  aagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGt                                                        <  1:5120692/93‑1 (MQ=255)
                                  aagCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGt                                                        >  2:5378294/1‑93 (MQ=255)
                                   agCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTc                                                       >  1:3929550/1‑93 (MQ=255)
                                   agCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTc                                                       <  1:7116933/93‑1 (MQ=255)
                                    gCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTAAGGGCGTGCCGTCg                                                      >  2:4945873/1‑93 (MQ=255)
                                    gCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTGTATGGGCGTGCCGTCg                                                      >  1:1009862/1‑93 (MQ=255)
                                      tATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGc                                                    >  2:5404912/1‑93 (MQ=255)
                                       aTTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCa                                                   >  2:6452409/1‑93 (MQ=255)
                                       aTTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCa                                                   >  1:144324/1‑93 (MQ=255)
                                        ttGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCaa                                                  >  2:6632340/1‑93 (MQ=255)
                                        ttGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCaa                                                  <  2:890275/93‑1 (MQ=255)
                                        ttGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCaa                                                  <  1:4620026/93‑1 (MQ=255)
                                          gATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATc                                                >  2:1386605/1‑93 (MQ=255)
                                          gATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATc                                                <  1:5746293/93‑1 (MQ=255)
                                          gATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATc                                                >  2:3860176/1‑93 (MQ=255)
                                            tttCTCCCGTCAGATGCGCGAAATCCCCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGa                                              <  2:5523392/93‑1 (MQ=255)
                                             ttCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGAc                                             >  2:6210759/1‑93 (MQ=255)
                                              tctcCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACa                                            <  1:1583533/93‑1 (MQ=255)
                                               ctcCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACAt                                           >  2:6489257/1‑93 (MQ=255)
                                               ctcCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACAt                                           >  2:1205241/1‑93 (MQ=255)
                                                  ccGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCAc                                        >  2:2307682/1‑93 (MQ=255)
                                                  ccGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCAc                                        >  2:3250589/1‑93 (MQ=255)
                                                  ccGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCAc                                        >  2:5814556/1‑93 (MQ=255)
                                                   cGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCAcc                                       <  2:6962886/93‑1 (MQ=255)
                                                    gTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCAccc                                      <  1:4938254/93‑1 (MQ=255)
                                                    gTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCAccc                                      >  2:3121323/1‑93 (MQ=255)
                                                    gTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCAccc                                      <  2:7392079/93‑1 (MQ=255)
                                                     tCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCcg                                     >  1:95686/1‑93 (MQ=255)
                                                     tCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCcg                                     >  2:597573/1‑93 (MQ=255)
                                                         aTGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTg                                 <  1:3008374/93‑1 (MQ=255)
                                                         aTGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTg                                 <  2:290977/93‑1 (MQ=255)
                                                          tGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTgg                                >  2:3822288/1‑93 (MQ=255)
                                                          tGCGCGAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTgg                                <  1:7072407/93‑1 (MQ=255)
                                                             gcgAAATCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTa                             >  1:52327/1‑93 (MQ=255)
                                                                aaaTCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTagca                          >  1:5428234/1‑93 (MQ=255)
                                                                 aaTCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTagcag                         >  2:7083927/1‑93 (MQ=255)
                                                                 aaTCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCAGCAATCGACATCACCCGCGTGGGTagcag                         <  1:1490494/93‑1 (MQ=255)
                                                                  aTCCCCCTCGGTGAAAACCCCAACATCCCCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCCCCCGCGTGGGTagcagc                        <  1:241489/93‑1 (MQ=255)
                                                                  aTCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTagcagc                        >  2:3258364/1‑93 (MQ=255)
                                                                  aTCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTagcagc                        >  2:2693677/1‑93 (MQ=255)
                                                                  aTCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTagcagc                        >  1:1132427/1‑93 (MQ=255)
                                                                  aTCACCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGCGCCGTCGGCAATCGACATCACCCGCGTGGGTagcagc                        <  2:6832693/93‑1 (MQ=255)
                                                                     aCCCTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGc                     <  2:2994845/93‑1 (MQ=255)
                                                                      tcgtCGGTGAAAACCCCAACATCGCCATTCCGCTGCTCGTTTTTATGGGCGTGCCGTCGTCAATCGGCATCACCCGCGTGGGTAGCAGCGGCa                    <  2:1858818/90‑1 (MQ=255)
                                                                      cccTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCa                    >  1:822026/1‑93 (MQ=255)
                                                                      cccTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCa                    <  1:613372/93‑1 (MQ=255)
                                                                       ccTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCAt                   <  2:4339096/93‑1 (MQ=255)
                                                                       ccTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCAt                   <  2:6335490/93‑1 (MQ=255)
                                                                       ccTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCAt                   >  1:6583560/1‑93 (MQ=255)
                                                                        cTCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCAtt                  <  1:2813430/93‑1 (MQ=255)
                                                                         tCGGTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTc                 <  1:4241843/93‑1 (MQ=255)
                                                                            gTGAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGc              <  2:6268070/93‑1 (MQ=255)
                                                                              gAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCg            >  1:5261481/1‑93 (MQ=255)
                                                                              gAAAACCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCg            >  2:3641436/1‑93 (MQ=255)
                                                                               aaaaCCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATAACCCGCGTGGGGAGCAGCGGCATTCTGCCgg           >  1:6987523/1‑93 (MQ=255)
                                                                                aaaCCCCAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGt          >  2:4828851/1‑93 (MQ=255)
                                                                                   ccccAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGAt       <  2:4662455/93‑1 (MQ=255)
                                                                                    cccAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATc      <  1:4917797/93‑1 (MQ=255)
                                                                                    cccAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATc      >  1:6697755/1‑93 (MQ=255)
                                                                                      cAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCaa    >  2:6740509/1‑93 (MQ=255)
                                                                                      cAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCaa    >  2:2275983/1‑93 (MQ=255)
                                                                                      cAACATCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCaa    >  1:7193592/1‑93 (MQ=255)
                                                                                        aCGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAaca  >  1:6208049/1‑93 (MQ=255)
                                                                                          |                                                                                          
GCTATCGTCGCGCTGGTGGGCGGCACGGTGGAAGAAGCTATTGATTTCTCCCGTCAGATGCGCGAAATCACCCTCGGTGAAAACCCCAACGTCACCATTCCGCTGCTCGGTTTTATGGGCGTGCCGTCGGCAATCGACATCACCCGCGTGGGTAGCAGCGGCATTCTGCCGGTGATCAACA  >  NC_000913/338671‑338851

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