Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,280,531 A→G Y482C (TAC→TGC)  narG → nitrate reductase 1 alpha subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,280,5310AG99.5% 661.3 / ‑4.8 207Y482C (TAC→TGC) narGnitrate reductase 1 alpha subunit
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (95/110);  total (96/111)
Fisher's exact test for biased strand distribution p-value = 4.66e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.67e-01

Only 100 of 213 total aligned reads displayed.

TGCTGCACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATACCCCGGCCTGGGCCGAGCAGATTA  >  NC_000913/1280441‑1280620
                                                                                          |                                                                                         
tgctggACAACCTGCCGGTGAAACTCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCg                                                                                         <  2:1842845/93‑1 (MQ=255)
tgctgcACAAACTGCCGGTGAACCGCCTGCAACTGGATGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCg                                                                                         <  2:4744770/93‑1 (MQ=255)
tgctgcACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCg                                                                                         >  2:4876180/1‑93 (MQ=255)
tgctgcACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCg                                                                                         >  2:1343771/1‑93 (MQ=255)
tgctgcACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCg                                                                                         >  1:2748776/1‑93 (MQ=255)
  ctgcACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGt                                                                                       >  1:784327/1‑93 (MQ=255)
   tgcACAAACTGCCGGTGAAACGCCTGCATCTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTc                                                                                      <  2:3484668/93‑1 (MQ=255)
   tgcACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTc                                                                                      <  1:1150932/93‑1 (MQ=255)
    gcACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCt                                                                                     >  1:3940336/1‑93 (MQ=255)
     cacaAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTg                                                                                    <  1:232148/93‑1 (MQ=255)
      acaAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTgg                                                                                   >  1:3292265/1‑93 (MQ=255)
       caAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGa                                                                                  <  1:5522248/93‑1 (MQ=255)
        aaaCTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGaa                                                                                 >  2:2314671/1‑93 (MQ=255)
        aaaCTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGaa                                                                                 <  1:2982953/93‑1 (MQ=255)
        aaaCTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGaa                                                                                 >  1:6214247/1‑93 (MQ=255)
         aaCTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAAc                                                                                <  1:3253587/93‑1 (MQ=255)
          aCTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACg                                                                               >  2:1114410/1‑93 (MQ=255)
            tGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTg                                                                             >  1:2355865/1‑93 (MQ=255)
             gCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACGGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGACCGTgg                                                                            <  2:1632427/93‑1 (MQ=255)
              ccGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGc                                                                           >  1:583307/1‑93 (MQ=255)
              ccGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGc                                                                           >  2:854786/1‑93 (MQ=255)
               cGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGcc                                                                          <  2:2840140/93‑1 (MQ=255)
               cGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGcc                                                                          >  2:4858229/1‑93 (MQ=255)
                ggTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTCTGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCt                                                                         <  1:4552570/93‑1 (MQ=255)
                ggTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCt                                                                         <  1:5114973/93‑1 (MQ=255)
                ggTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCt                                                                         >  2:6250827/1‑93 (MQ=255)
                  tGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGa                                                                       <  2:4464656/93‑1 (MQ=255)
                  tGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGa                                                                       <  1:2240965/93‑1 (MQ=255)
                  tGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGa                                                                       <  2:1968808/93‑1 (MQ=255)
                   gAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGaa                                                                      <  2:4948558/93‑1 (MQ=255)
                    aaaCGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAac                                                                     <  2:1487493/93‑1 (MQ=255)
                    aaaCGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAac                                                                     <  1:1927664/93‑1 (MQ=255)
                       cGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAacgac                                                                  <  2:324717/93‑1 (MQ=255)
                        gCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAacgacg                                                                 >  2:2247743/1‑93 (MQ=255)
                         ccTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGt                                                                >  1:1022099/1‑93 (MQ=255)
                          cTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGtt                                                               <  2:3859741/93‑1 (MQ=255)
                          cTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGtt                                                               >  1:4875027/1‑93 (MQ=255)
                           tGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTa                                                              >  2:1315029/1‑93 (MQ=255)
                            gCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTaa                                                             >  1:4331479/1‑93 (MQ=255)
                              aaCTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACt                                                           <  1:4991302/93‑1 (MQ=255)
                               aCTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACtg                                                          >  2:200608/1‑93 (MQ=255)
                                  ggCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGc                                                       <  2:1804098/93‑1 (MQ=255)
                                  ggCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGc                                                       >  1:3743006/1‑93 (MQ=255)
                                   gCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGCCGTTAACTGTGCa                                                      >  1:5660141/1‑93 (MQ=255)
                                      gATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCCAcc                                                   >  2:5942309/1‑93 (MQ=255)
                                      gATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAAcc                                                   >  1:656639/1‑93 (MQ=255)
                                        tGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGGTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAg                                                 >  2:5133349/1‑93 (MQ=255)
                                        tGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAg                                                 >  1:2902562/1‑93 (MQ=255)
                                        tGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAg                                                 >  2:741497/1‑93 (MQ=255)
                                         ggCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGc                                                <  1:376364/93‑1 (MQ=255)
                                          gcagcaCCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCt                                               <  1:418907/93‑1 (MQ=255)
                                           cagcaCCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTa                                              <  2:4327843/93‑1 (MQ=255)
                                           cagcaCCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTa                                              <  1:1722174/93‑1 (MQ=255)
                                             gcaCCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATg                                            >  2:3003313/1‑93 (MQ=255)
                                             gcaCCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATg                                            >  2:3431833/1‑93 (MQ=255)
                                              caCCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGa                                           <  1:997726/93‑1 (MQ=255)
                                                ccGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACg                                         >  2:1698314/1‑93 (MQ=255)
                                                   cccTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGAtg                                      <  1:3322875/93‑1 (MQ=255)
                                                    ccgggTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGACTGGAACGTGGCGTGAACGACGTTAACTGTGCAACCAGCTATGACGAtgt                                     <  2:619106/90‑1 (MQ=255)
                                                    ccTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGAtgt                                     <  1:2486979/93‑1 (MQ=255)
                                                    ccTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGAtgt                                     <  2:3549160/93‑1 (MQ=255)
                                                    ccTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGAtgt                                     >  2:6031998/1‑93 (MQ=255)
                                                     cTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGAtgtg                                    >  2:3335482/1‑93 (MQ=255)
                                                      tGGTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGa                                   >  1:97470/1‑93 (MQ=255)
                                                       ggTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGaa                                  >  1:3640582/1‑93 (MQ=255)
                                                        gTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGaaa                                 <  1:3597239/93‑1 (MQ=255)
                                                        gTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGaaa                                 <  1:1359673/93‑1 (MQ=255)
                                                        gTGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGaaa                                 <  1:282223/93‑1 (MQ=255)
                                                         tGACCACCGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAg                                <  1:4378810/93‑1 (MQ=255)
                                                            ccaccGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCtt                             >  2:2500263/1‑93 (MQ=255)
                                                            ccaccGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCtt                             <  1:3352597/93‑1 (MQ=255)
                                                            ccaccGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCtt                             >  1:2681515/1‑93 (MQ=255)
                                                            ccaccGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCtt                             >  1:723892/1‑93 (MQ=255)
                                                             caccGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTta                            <  1:83571/93‑1 (MQ=255)
                                                             caccGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTta                            >  1:1283819/1‑93 (MQ=255)
                                                             caccGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTta                            <  1:2834093/93‑1 (MQ=255)
                                                              accGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTtat                           >  1:305156/1‑93 (MQ=255)
                                                                cGTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATAc                         <  1:6252570/93‑1 (MQ=255)
                                                                 gTTTATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATAcc                        >  2:2901987/1‑93 (MQ=255)
                                                                  tttATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATAccc                       <  2:5618186/93‑1 (MQ=255)
                                                                  tttATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATAccc                       <  2:367204/93‑1 (MQ=255)
                                                                   ttATGATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATAcccc                      <  2:3799945/93‑1 (MQ=255)
                                                                     aTGATCTGACGCTGGCAAACTGCGGTCTGCAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATACCCCgg                    <  2:2728911/93‑1 (MQ=255)
                                                                       gATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATACCCCGGcc                  >  1:6009611/1‑93 (MQ=255)
                                                                       gATCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATACCCCGGcc                  <  1:4753306/93‑1 (MQ=255)
                                                                        aTCTGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATACCCCGGCCt                 >  2:108884/1‑93 (MQ=255)
                                                                         tCTGACGCTGGCAAACTGCGGTCTGAAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATACCCCGGCCTg                <  1:3027521/93‑1 (MQ=255)
                                                                           tGACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATACCCCGGCCTggg              >  1:1916803/1‑93 (MQ=255)
                                                                            gACGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATACCCCGGCCTGGGc             <  2:1267268/93‑1 (MQ=255)
                                                                             aCGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATACCCCGGCCTGGGcc            >  1:3886423/1‑93 (MQ=255)
                                                                              cGCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATACCCCGGCCTGGGCCg           <  2:4151534/93‑1 (MQ=255)
                                                                               gCTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATACCCCGGCCTGGGCCGa          >  2:5384623/1‑93 (MQ=255)
                                                                                cTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATACCCCGGCCTGGGCCgtg         >  2:5066613/1‑91 (MQ=255)
                                                                                cTGGCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATACCCCGGCCTGGGCCGAg         <  2:792217/93‑1 (MQ=255)
                                                                                  ggCAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATACCCCGGCCTGGGCCGAGCa       <  2:3923288/93‑1 (MQ=255)
                                                                                  ggCAAACTACGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATACCCCGGCCTGGGCCGAGCa       >  2:6167064/1‑93 (MQ=255)
                                                                                    cAAACTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATACCCCGGCCTGGGCCGAGCAGa     <  1:1614331/93‑1 (MQ=255)
                                                                                     aaaCTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATACCCCGGCCTGGGCCGAGCAGAt    <  2:440512/93‑1 (MQ=255)
                                                                                     aaaCTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATACCCCGGCCTGGGCCGAGCAGAt    <  1:4662711/93‑1 (MQ=255)
                                                                                      aaCTGCGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATACCCCGGCCTGGGCCGAGCAGAtt   >  1:5355514/1‑93 (MQ=255)
                                                                                          |                                                                                         
TGCTGCACAAACTGCCGGTGAAACGCCTGCAACTGGCTGATGGCAGCACCGCCCTGGTGACCACCGTTTATGATCTGACGCTGGCAAACTACGGTCTGGAACGTGGCCTGAACGACGTTAACTGTGCAACCAGCTATGACGATGTGAAAGCTTATACCCCGGCCTGGGCCGAGCAGATTA  >  NC_000913/1280441‑1280620

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