Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 1,535,627 A→G R417R (CGT→CGC narY ← cryptic nitrate reductase 2 beta subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,535,6270AG99.5% 659.8 / ‑5.0 204R417R (CGT→CGCnarYcryptic nitrate reductase 2 beta subunit
Reads supporting (aligned to +/- strand):  ref base A (1/0);  new base G (106/97);  total (107/97)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.81e-01

Only 100 of 214 total aligned reads displayed.

GAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATATTGGCGAGATATTG  >  NC_000913/1535537‑1535716
                                                                                          |                                                                                         
gagaTAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCg                                                                                         >  1:3480774/1‑93 (MQ=255)
 agaTAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAgcgc                                                                                        >  2:4645829/1‑93 (MQ=255)
 agaTAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAgcgc                                                                                        <  2:616232/93‑1 (MQ=255)
  gaTAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCa                                                                                       >  1:3840587/1‑93 (MQ=255)
  gaTAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCa                                                                                       <  1:3471476/93‑1 (MQ=255)
    tAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCata                                                                                     <  1:1264449/93‑1 (MQ=255)
    tAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGACTACTTCGTCGATGGCACGAGTCTGAGTAACGCCTTCCACGGTTTCTGAGCGCata                                                                                     >  1:4031140/1‑93 (MQ=255)
    tAGCGATACATCTCTGCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCata                                                                                     <  1:141926/93‑1 (MQ=255)
     aGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCatat                                                                                    <  2:1675706/93‑1 (MQ=255)
      gCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCAtata                                                                                   <  2:2007510/93‑1 (MQ=255)
      gCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCAtata                                                                                   <  2:3675817/93‑1 (MQ=255)
       cGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAg                                                                                  <  2:5301365/93‑1 (MQ=255)
        gATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGt                                                                                 <  1:5596206/93‑1 (MQ=255)
        gATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGt                                                                                 <  2:2129549/93‑1 (MQ=255)
         aTACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTg                                                                                >  2:6546581/1‑93 (MQ=255)
         aTACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTg                                                                                <  2:1634770/93‑1 (MQ=255)
              tctcTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCa                                                                           >  2:4001929/1‑93 (MQ=255)
              tctcTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCa                                                                           <  2:1846491/93‑1 (MQ=255)
              tctcTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCa                                                                           <  2:1475088/93‑1 (MQ=255)
                 ctTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCg                                                                        >  1:4690413/1‑93 (MQ=255)
                  ttCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGc                                                                       <  2:1021022/93‑1 (MQ=255)
                   tCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGcc                                                                      <  1:421538/93‑1 (MQ=255)
                   tCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGcc                                                                      <  1:2441490/93‑1 (MQ=255)
                    cGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCca                                                                     <  1:278205/93‑1 (MQ=255)
                    cGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCca                                                                     >  1:2732457/1‑93 (MQ=255)
                    cGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCca                                                                     >  1:2590111/1‑93 (MQ=255)
                     gACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCcat                                                                    <  2:291223/93‑1 (MQ=255)
                      aCCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCcatc                                                                   >  2:56723/1‑93 (MQ=255)
                      aCCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCcatc                                                                   <  2:2704633/93‑1 (MQ=255)
                      aCCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCcatc                                                                   <  2:2863769/93‑1 (MQ=255)
                          gggCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATcc                                                               >  1:941248/1‑93 (MQ=255)
                          gggCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATcc                                                               <  2:2185133/93‑1 (MQ=255)
                          gggCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATcc                                                               >  2:3811815/1‑93 (MQ=255)
                           ggCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCg                                                              <  2:2567311/93‑1 (MQ=255)
                            gCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGt                                                             >  2:2214165/1‑93 (MQ=255)
                            gCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGt                                                             <  2:3318101/93‑1 (MQ=255)
                              gACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGttt                                                           <  2:4477466/93‑1 (MQ=255)
                               aCGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTc                                                          <  2:6039773/93‑1 (MQ=255)
                               aCGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTc                                                          >  1:5807601/1‑93 (MQ=255)
                                 gCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAg                                                        <  1:5108184/93‑1 (MQ=255)
                                 gCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAg                                                        <  1:3432117/93‑1 (MQ=255)
                                  cTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAgc                                                       <  2:6239980/93‑1 (MQ=255)
                                    cAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCTGTTTCAgcgc                                                     >  1:4233616/1‑93 (MQ=255)
                                      ggCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGccc                                                   <  1:360460/93‑1 (MQ=255)
                                      ggCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGccc                                                   <  1:2872776/93‑1 (MQ=255)
                                         cTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCa                                                >  2:1444805/1‑93 (MQ=255)
                                         cTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCa                                                >  1:3019202/1‑93 (MQ=255)
                                          tACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAg                                               >  2:5089000/1‑93 (MQ=255)
                                          tACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAg                                               >  2:205990/1‑93 (MQ=255)
                                          tACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAg                                               <  1:1736593/93‑1 (MQ=255)
                                           aCTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGt                                              <  1:1792557/93‑1 (MQ=255)
                                            cTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTa                                             >  1:4913954/1‑93 (MQ=255)
                                            cTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTa                                             >  2:4863081/1‑93 (MQ=255)
                                            cTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTa                                             >  1:5189883/1‑93 (MQ=255)
                                            cTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGCTTCAGCGCCCGCAGTa                                             >  2:1175450/1‑93 (MQ=255)
                                             ttCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTac                                            >  1:1349713/1‑93 (MQ=255)
                                               cgtcgATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTaccg                                          >  1:6432226/1‑93 (MQ=255)
                                               cgtcgATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTaccg                                          <  1:3648448/93‑1 (MQ=255)
                                                 tcgATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTaccgga                                        >  2:5555315/1‑93 (MQ=255)
                                                 tcgATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTaccgga                                        >  1:6011587/1‑93 (MQ=255)
                                                  cgATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTaccggac                                       >  2:5861402/1‑93 (MQ=255)
                                                  cgATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTaccggac                                       >  2:5000460/1‑93 (MQ=255)
                                                   gATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTaccggacc                                      >  1:4947630/1‑93 (MQ=255)
                                                   gATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTaccggacc                                      >  2:4246400/1‑93 (MQ=255)
                                                   gATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTaccggacc                                      <  2:6328787/93‑1 (MQ=255)
                                                     tGGCACGAGTATCAGTAACGCCTTCCACGGTTTCTGAGCGCATATAGTGGCGCATCGCCATCATCCGCTTCAGCGCCCGCAGTaccggaccgg                                    >  2:2239679/1‑93 (MQ=255)
                                                      ggCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGt                                   >  1:807903/1‑93 (MQ=255)
                                                      ggCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGt                                   >  1:3707202/1‑93 (MQ=255)
                                                       gCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTa                                  <  2:5227679/93‑1 (MQ=255)
                                                       gCACGAGTATCAGTAACCCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTa                                  <  1:270563/93‑1 (MQ=255)
                                                        cACGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTAt                                 >  1:6716053/1‑93 (MQ=255)
                                                         aCGAGTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATc                                <  1:1905397/93‑1 (MQ=255)
                                                             gTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCg                            >  2:4767289/1‑93 (MQ=255)
                                                             gTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCg                            >  2:5707124/1‑93 (MQ=255)
                                                             gTATCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCg                            >  1:4999691/1‑93 (MQ=255)
                                                               aTCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGc                          <  1:5745684/93‑1 (MQ=255)
                                                                tCAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGAGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGACGGCa                         >  1:91870/1‑93 (MQ=255)
                                                                 cAGTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCAc                        >  2:800588/1‑93 (MQ=255)
                                                                   gTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTc                      >  1:2920419/1‑93 (MQ=255)
                                                                   gTAACGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTc                      >  1:3028122/1‑93 (MQ=255)
                                                                       cGCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACa                  <  1:3117603/93‑1 (MQ=255)
                                                                        gCCTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACat                 >  2:376987/1‑93 (MQ=255)
                                                                         ccTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACata                >  1:4449971/1‑93 (MQ=255)
                                                                         ccTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACata                <  1:6191805/93‑1 (MQ=255)
                                                                          cTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACatat               >  2:352115/1‑93 (MQ=255)
                                                                          cTTCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACatat               >  1:2102094/1‑93 (MQ=255)
                                                                           ttCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATAtt              <  1:11504/93‑1 (MQ=255)
                                                                           ttCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATAtt              <  2:5797050/93‑1 (MQ=255)
                                                                            tCCACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATATTg             <  1:2296622/93‑1 (MQ=255)
                                                                             ccACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATATTgg            <  1:3569019/93‑1 (MQ=255)
                                                                              cACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATATTGGc           <  1:693714/93‑1 (MQ=255)
                                                                              cACGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATATTGGc           >  1:3330867/1‑93 (MQ=255)
                                                                               aCGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATATTGGCg          >  2:343951/1‑93 (MQ=255)
                                                                                cGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATATTGGCga         <  2:4711180/93‑1 (MQ=255)
                                                                                cGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATATTGGCga         <  2:2645540/93‑1 (MQ=255)
                                                                                cGGTTTGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATATTGGCga         <  1:1838594/93‑1 (MQ=255)
                                                                                    ttGTGAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATATTGGCGAGata     <  2:713491/93‑1 (MQ=255)
                                                                                     tgtgAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATATTGGCGAGatat    <  1:1275872/93‑1 (MQ=255)
                                                                                      gtgAGCTCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATATTGGCGAGATAtt   <  1:212963/93‑1 (MQ=255)
                                                                                      gtgAGCGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATATTGGCGAGATAtt   <  2:242202/93‑1 (MQ=255)
                                                                                          |                                                                                         
GAGATAGCGATACATCTCTTCGACCTGGGCGACGCTCAGGCCTACTTCGTCGATGGCACGAGTATCAGTAACGCCTTCCACGGTTTGTGAACGCATATAGTGGCGCATCGCCATCATCCGTTTCAGCGCCCGCAGTACCGGACCGGTATCGCCGGCACTCAACATATTGGCGAGATATTG  >  NC_000913/1535537‑1535716

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: