Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,827,618 A→G R47R (CGT→CGC b1746 ← putative aldehyde dehydrogenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,827,6180AG96.8% 283.4 / ‑0.7 95R47R (CGT→CGCb1746putative aldehyde dehydrogenase
Reads supporting (aligned to +/- strand):  ref base A (2/1);  new base G (47/45);  total (49/46)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.11e-01

Only 100 of 101 total aligned reads displayed.

CGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTC  >  NC_000913/1827527‑1827705
                                                                                           |                                                                                       
cGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCAc                                                                                        <  2:3846023/93‑1 (MQ=255)
 ggCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCcgcgcg                                                                                       <  2:3337519/93‑1 (MQ=255)
  gCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCgg                                                                                      >  1:3126444/1‑93 (MQ=255)
  gCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCgg                                                                                      >  2:1000885/1‑93 (MQ=255)
  gCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGAGGAAACGCCGCGCgg                                                                                      >  1:4181021/1‑93 (MQ=255)
    tttATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGc                                                                                    >  1:1229202/1‑93 (MQ=255)
      tATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTg                                                                                  >  1:120739/1‑93 (MQ=255)
      tATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTg                                                                                  <  2:817024/93‑1 (MQ=255)
        ttGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGcc                                                                                <  2:2928366/93‑1 (MQ=255)
         tGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGccc                                                                               >  1:880285/1‑93 (MQ=255)
         tGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGccc                                                                               >  2:135517/1‑93 (MQ=255)
         tGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGccc                                                                               >  2:1448229/1‑93 (MQ=255)
         tGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGccc                                                                               >  2:1015523/1‑93 (MQ=255)
         tGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGccc                                                                               >  2:4638950/1‑93 (MQ=255)
          gCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCg                                                                              >  1:1290325/1‑93 (MQ=255)
           cTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGa                                                                             >  2:262689/1‑93 (MQ=255)
            tttCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGAc                                                                            >  2:2338058/1‑93 (MQ=255)
              tCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACaa                                                                          <  2:2461266/93‑1 (MQ=255)
               ccAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAg                                                                         >  1:2568559/1‑93 (MQ=255)
                 agcagTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGcc                                                                       >  1:3040799/1‑93 (MQ=255)
                  gcagTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCt                                                                      <  1:2639230/93‑1 (MQ=255)
                  gcagTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCt                                                                      <  2:2376891/93‑1 (MQ=255)
                    agTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGc                                                                    >  1:102494/1‑93 (MQ=255)
                    agTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGc                                                                    >  1:4431845/1‑93 (MQ=255)
                     gTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCt                                                                   >  1:3967929/1‑93 (MQ=255)
                     gTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCt                                                                   >  1:3141825/1‑93 (MQ=255)
                       gCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCg                                                                 <  1:3577971/93‑1 (MQ=255)
                       gCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCg                                                                 <  2:3781416/93‑1 (MQ=255)
                        cGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGa                                                                <  1:3331145/93‑1 (MQ=255)
                          gCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGAcc                                                              <  2:4042507/93‑1 (MQ=255)
                           cAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCt                                                             <  1:508931/93‑1 (MQ=255)
                           cAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCt                                                             <  1:508930/93‑1 (MQ=255)
                           cAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCt                                                             <  1:3819879/93‑1 (MQ=255)
                            aaaGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTg                                                            >  2:2744921/1‑93 (MQ=255)
                            aaaGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTg                                                            <  1:3560877/93‑1 (MQ=255)
                            aaaGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTg                                                            <  1:3560976/93‑1 (MQ=255)
                             aaGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGa                                                           <  2:1546801/93‑1 (MQ=255)
                             aaGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGa                                                           >  1:2114607/1‑93 (MQ=255)
                             aaGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGa                                                           <  2:786671/93‑1 (MQ=255)
                             aaGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGa                                                           <  1:1176366/93‑1 (MQ=255)
                             aaGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCCGCTCGACCTGa                                                           >  1:3874657/1‑93 (MQ=255)
                               gCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGc                                                         <  2:4726975/93‑1 (MQ=255)
                                   tccacacCGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggca                                                     <  1:1716018/90‑1 (MQ=255)
                                   tCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggca                                                     <  1:133556/93‑1 (MQ=255)
                                    cGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggcat                                                    >  2:2023234/1‑93 (MQ=255)
                                    cGACAACGGCATGACGTGCAGCAAATGAAAGCCGCGCCCAGCGCGGAACCGCCGCGCGGGCTGACCGACAAGCCTGGTCGACCTGAGcggcat                                                    <  2:357287/93‑1 (MQ=255)
                                      aCAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggcatcg                                                  >  2:1366721/1‑93 (MQ=255)
                                       cAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGCAACGCCGCGCGGGCTGCCCGACAAGACTGCTCGACCCGAGcgccctcgg                                                 >  2:1310286/1‑93 (MQ=255)
                                        aaCGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggcatcggc                                                <  1:2127626/93‑1 (MQ=255)
                                          cGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggcatcggcat                                              <  1:3023776/93‑1 (MQ=255)
                                          cGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGcggactcggcat                                              >  1:4484713/1‑93 (MQ=255)
                                              aTGACGTTCAGCAAATGAAGGCCGCGCCCGGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAtt                                          <  2:4190111/93‑1 (MQ=255)
                                              aTGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAtt                                          <  2:3433225/93‑1 (MQ=255)
                                              aTGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAtt                                          >  1:509962/1‑93 (MQ=255)
                                              aTGACGTTCAGCAAATGAAAGCCGAGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGACTGCTCGACCTGAGCGGCATCGGcatcatc                                          >  1:2491858/1‑92 (MQ=255)
                                               tGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAttg                                         >  2:1748470/1‑93 (MQ=255)
                                               tGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAttg                                         >  2:637379/1‑93 (MQ=255)
                                               tGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAttg                                         >  1:1946719/1‑93 (MQ=255)
                                                gACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAttgc                                        >  2:1893403/1‑93 (MQ=255)
                                                   gTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAttgcctt                                     >  1:4114/1‑93 (MQ=255)
                                                   gTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAttgcctt                                     >  1:2498273/1‑93 (MQ=255)
                                                   gTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAttgcctt                                     >  2:2719787/1‑93 (MQ=255)
                                                   gTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAttgcctt                                     <  2:3943085/93‑1 (MQ=255)
                                                   gTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGAACTGAGCGGCATCGGCATCAttgcctt                                     <  2:2567317/93‑1 (MQ=255)
                                                    ttCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAttgccttg                                    >  2:2738305/1‑93 (MQ=255)
                                                    ttCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCAttgccttg                                    <  2:2336289/93‑1 (MQ=255)
                                                       aGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCa                                 >  2:1527038/1‑93 (MQ=255)
                                                       aGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCa                                 >  1:1240849/1‑93 (MQ=255)
                                                        gCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCAt                                <  1:2969334/93‑1 (MQ=255)
                                                         cAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATa                               <  1:4653297/93‑1 (MQ=255)
                                                           aaTGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAac                             >  2:4581942/1‑93 (MQ=255)
                                                           aaTGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAac                             >  1:4700306/1‑93 (MQ=255)
                                                           aaTGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAac                             <  1:1461609/93‑1 (MQ=255)
                                                           aaTGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAac                             <  2:4499138/93‑1 (MQ=255)
                                                            aTGAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAaca                            >  1:1752641/1‑93 (MQ=255)
                                                              gAAAGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACAcc                          <  1:32750/93‑1 (MQ=255)
                                                               aaaGCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCt                         <  2:2483476/93‑1 (MQ=255)
                                                                  gCCGCGCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGc                      <  2:3256219/93‑1 (MQ=255)
                                                                    cgcgCCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGccc                    <  1:33496/93‑1 (MQ=255)
                                                                    cgcaCCCAGCGTGCAACGCCCGCGCGGGCTGCCCGACACGCCTGCGCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGccc                    <  1:33495/93‑1 (MQ=255)
                                                                     gcgcCCAGCGCGGAAACGCCGCGCGGTCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCg                   <  2:2090289/93‑1 (MQ=255)
                                                                      cgcCCAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGa                  <  1:4190341/93‑1 (MQ=255)
                                                                       ccccAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGAt                 <  2:3922184/92‑1 (MQ=255)
                                                                        cccAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATa                >  2:3142356/1‑93 (MQ=255)
                                                                         ccAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATAc               >  2:4051505/1‑93 (MQ=255)
                                                                         ccAGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATAc               >  1:1153174/1‑93 (MQ=255)
                                                                           aGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCg             <  2:767601/93‑1 (MQ=255)
                                                                           aGCGCGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCg             >  1:2458307/1‑93 (MQ=255)
                                                                            gcgcGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCgg            >  1:2322870/1‑93 (MQ=255)
                                                                            gcgcGGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCgg            >  1:1975568/1‑93 (MQ=255)
                                                                             cgcgGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGa           >  1:768409/1‑93 (MQ=255)
                                                                              gcgcAAACGCCGCGCGGGCTCCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGAt          <  1:259127/93‑1 (MQ=255)
                                                                              gcgGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGAt          >  2:869072/1‑93 (MQ=255)
                                                                               cgGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGAtt         <  2:1874183/93‑1 (MQ=255)
                                                                               cgGAAACGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGAtt         <  1:3244957/93‑1 (MQ=255)
                                                                                   aaCGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGc     >  2:3115576/1‑93 (MQ=255)
                                                                                    aCGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCt    <  2:1269028/93‑1 (MQ=255)
                                                                                    aCGCCGCGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCt    <  1:2252009/93‑1 (MQ=255)
                                                                                      gccgcGCGGGCTGCCCGCCAAGCCTGCTCGCCCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTc  <  1:370968/93‑1 (MQ=255)
                                                                                      gccgcGCGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTc  >  2:2001334/1‑93 (MQ=255)
                                                                                           |                                                                                       
CGGCTTTATTGCTTTCCAGCAGTGCGGCAAAGCGTTCGACAACGGCATGACGTTCAGCAAATGAAAGCCGCGCCCAGCGCGGAAACGCCGCACGGGCTGCCCGACAAGCCTGCTCGACCTGAGCGGCATCGGCATCATTGCCTTGCCATAACACCTCGCCCGATACCGGATTACGCTTC  >  NC_000913/1827527‑1827705

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