Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,356,616 T→C A1450A (GCT→GCC)  gltB → glutamate synthase large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,356,6160TC97.8% 425.4 / ‑0.5 137A1450A (GCT→GCC) gltBglutamate synthase large subunit
Reads supporting (aligned to +/- strand):  ref base T (2/1);  new base C (65/69);  total (67/70)
Fisher's exact test for biased strand distribution p-value = 6.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.33e-01

Only 100 of 144 total aligned reads displayed.

GACCGGCGGTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCTCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATACCGGCTCTCAGCGCGGTGAAGAGATTCTGGCGAA  >  NC_000913/3356526‑3356705
                                                                                          |                                                                                         
taCCGGCGGTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGccct                                                                                         <  1:2630896/92‑1 (MQ=255)
 aCCGGCGGTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTg                                                                                        <  1:348031/93‑1 (MQ=255)
 aCCGGCGGTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTg                                                                                        >  1:2765061/1‑93 (MQ=255)
      cggTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGAt                                                                                   <  1:2608738/93‑1 (MQ=255)
       ggTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATc                                                                                  >  2:4308189/1‑93 (MQ=255)
        gTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATcc                                                                                 >  1:3614402/1‑93 (MQ=255)
        gTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATcc                                                                                 >  2:2702261/1‑93 (MQ=255)
        gTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATcc                                                                                 >  1:2033638/1‑93 (MQ=255)
        gTTTCGCTTACGTTCTCGAGGAAAGCGGCGATTTCCGCAAACGCGTTAACCTGGAACTGGTCGAGGTCTGAAGCGTTGATGCCCTGGCGATcc                                                                                 >  1:4225900/1‑93 (MQ=255)
         tttCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCa                                                                                >  2:3340997/1‑93 (MQ=255)
         tttCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCa                                                                                <  2:1862893/93‑1 (MQ=255)
           tCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGTCGCTCTGGCGATCCATg                                                                              >  2:3129326/1‑93 (MQ=255)
           tCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATg                                                                              <  2:485344/93‑1 (MQ=255)
            cGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGa                                                                             <  1:803650/93‑1 (MQ=255)
            cGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGa                                                                             <  1:3078425/93‑1 (MQ=255)
               ttACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAaga                                                                          >  2:3932916/1‑93 (MQ=255)
                tACGTTCTCGATGAAAGCGGCGATTTCCGGAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTTACTCCCTGGCGATCCATGAagag                                                                         >  2:2618624/1‑93 (MQ=255)
                 aCGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGc                                                                        >  1:1591920/1‑93 (MQ=255)
                  cGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCa                                                                       >  1:2634021/1‑93 (MQ=255)
                   gTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCAt                                                                      >  2:3779060/1‑93 (MQ=255)
                    ttCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATc                                                                     <  1:8510/93‑1 (MQ=255)
                    ttCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATc                                                                     <  1:3124509/93‑1 (MQ=255)
                    ttCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATc                                                                     <  2:1794791/93‑1 (MQ=255)
                    ttCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATc                                                                     <  2:2381941/93‑1 (MQ=255)
                    ttCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGAATc                                                                     >  1:3413176/1‑93 (MQ=255)
                     tctcGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCt                                                                    >  1:3963479/1‑93 (MQ=255)
                     tctcGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCt                                                                    <  2:2283632/93‑1 (MQ=255)
                       tcGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGc                                                                  <  1:103559/93‑1 (MQ=255)
                       tcGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGc                                                                  >  2:4124190/1‑93 (MQ=255)
                       tcGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGc                                                                  <  2:728163/93‑1 (MQ=255)
                       tcGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGc                                                                  >  2:4748444/1‑93 (MQ=255)
                         gATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGt                                                                >  1:59913/1‑93 (MQ=255)
                         gATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGGCTTAAGCGTTGACGCCCTGTCGATCCATGAAGAGCATCTGCGt                                                                >  2:3394083/1‑93 (MQ=255)
                          aTGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTg                                                               <  2:3817708/93‑1 (MQ=255)
                          aTGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTg                                                               <  1:3627791/93‑1 (MQ=255)
                            gAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGt                                                             <  1:1073906/93‑1 (MQ=255)
                            gAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGt                                                             <  1:1073907/93‑1 (MQ=255)
                            gAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGt                                                             >  1:3993491/1‑93 (MQ=255)
                            gAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGt                                                             >  1:3746046/1‑93 (MQ=255)
                               aGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCtt                                                          >  1:1908676/1‑93 (MQ=255)
                                tcggcgATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTa                                                         >  2:2215561/2‑93 (MQ=255)
                                  ggcgATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATc                                                       >  2:3575111/1‑93 (MQ=255)
                                  ggcgATTTCCGCAAACGCGTTAACCCGGAACTGGGCGAGGTCTTAAGCGTAGACGCCCTGGCGATCCATGAAGAGCATCTGCGTCGTCTTATc                                                       >  2:3575110/1‑93 (MQ=255)
                                      aTTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCg                                                   >  2:4148680/1‑93 (MQ=255)
                                      aTTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCg                                                   <  2:2280498/93‑1 (MQ=255)
                                        ttCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAg                                                 >  2:4007077/1‑93 (MQ=255)
                                         tCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGc                                                <  1:1248654/93‑1 (MQ=255)
                                          ccGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCa                                               <  1:3571400/93‑1 (MQ=255)
                                          ccGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCa                                               >  2:3647688/1‑93 (MQ=255)
                                           cGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCTCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCAt                                              <  1:4030280/93‑1 (MQ=255)
                                            gCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCAtg                                             <  1:4513147/93‑1 (MQ=255)
                                            gCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCAtg                                             >  1:4323136/1‑93 (MQ=255)
                                              aaaCGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCAtgtg                                           <  2:1370828/93‑1 (MQ=255)
                                              aaaCGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCAtgtg                                           >  1:1933811/1‑93 (MQ=255)
                                               aaCGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTgc                                          >  2:4650684/1‑93 (MQ=255)
                                               aaCGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTgc                                          <  1:1048690/93‑1 (MQ=255)
                                               aaCGCCTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTgc                                          <  1:3974480/93‑1 (MQ=255)
                                                aCGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTgca                                         <  1:1323601/93‑1 (MQ=255)
                                                 cgcgTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTgcag                                        <  1:3595533/93‑1 (MQ=255)
                                                  gcgTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTgcagc                                       <  2:1794037/93‑1 (MQ=255)
                                                  gcgTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTgcagc                                       <  1:1721403/93‑1 (MQ=255)
                                                   cgTTAGCCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTgcagca                                      <  1:1954576/93‑1 (MQ=255)
                                                   cgTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTgcagca                                      <  2:2014700/93‑1 (MQ=255)
                                                    gTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCAt                                     >  2:2330923/1‑93 (MQ=255)
                                                    gTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCAt                                     <  2:1379343/93‑1 (MQ=255)
                                                     ttAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATa                                    <  1:704600/93‑1 (MQ=255)
                                                     ttAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATa                                    >  1:3243392/1‑93 (MQ=255)
                                                     ttAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATa                                    >  1:1035486/1‑93 (MQ=255)
                                                      tagcccGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATAc                                   <  1:3619926/90‑1 (MQ=255)
                                                      tAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATAc                                   <  2:4345414/93‑1 (MQ=255)
                                                      tAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATAc                                   >  2:1615480/1‑93 (MQ=255)
                                                        aCCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATACCg                                 >  1:3025202/1‑93 (MQ=255)
                                                        aCCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATACCg                                 >  1:591268/1‑93 (MQ=255)
                                                         cccGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATACCgg                                <  2:2534321/93‑1 (MQ=255)
                                                          ccGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATACCGGc                               >  1:4694195/1‑93 (MQ=255)
                                                          ccGGAACTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATACCGGc                               >  1:2333709/1‑93 (MQ=255)
                                                             gAACTGGTCTAGGTCTTAAGCGGTTACGCCCTTGCGATCCATGAAGAGCATCTTCGTGGTCTTATCACCGAGCATGTGCAGCATAACGGctct                            >  2:2823425/1‑93 (MQ=255)
                                                               aCTGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATACCGGCTCTCa                          >  1:2224194/1‑93 (MQ=255)
                                                                 tGGTCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATACCGGCTCTCAgc                        >  2:4431310/1‑93 (MQ=255)
                                                                   gTCGAGGTCTTAAGCGTTGGCGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATACCGGCTCTCAgcgc                      <  1:493848/93‑1 (MQ=255)
                                                                   gTCGAGGGCTTAAGCGTTGACGCCCCGGGGGTCCATGAAGAGGATCTGCGTGGGCTTATCAACCAGCATGTGGAGGATACCCGCTCTCCgcgc                      >  2:2773797/1‑93 (MQ=255)
                                                                    tCGAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATACCGGCTCTCAGcgcg                     <  2:1656833/93‑1 (MQ=255)
                                                                      gAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATACCGGCTCTCAGCGCGGt                   >  2:2534786/1‑93 (MQ=255)
                                                                      gAGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATACCGGCTCTCAGCGCGGt                   <  1:4471673/93‑1 (MQ=255)
                                                                       aGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATACCGGCTCTCAGCGCGGTg                  >  2:2354607/1‑93 (MQ=255)
                                                                       aGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATACCGGCTCTCAGCGCGGTg                  >  1:3006364/1‑93 (MQ=255)
                                                                       aGGTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATACCGGCTCTCAGCGCGGTg                  <  1:3602742/93‑1 (MQ=255)
                                                                         gTCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATACCGGCTCTCAGCGCGGTgca                >  2:2209660/1‑91 (MQ=255)
                                                                          tCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATACCGGCTCTCAGCGCGGTGAag               >  1:4285134/1‑93 (MQ=255)
                                                                          tCTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATACCGGCTCTCAGCGCGGTGAag               >  1:2799485/1‑93 (MQ=255)
                                                                           cTTAAGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATACCGGCTCTCAGCGCGGTGAaga              <  2:4563600/93‑1 (MQ=255)
                                                                              taGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATACCGGCTCTCAGCGCGGTGAAGAGAt           >  2:2882898/2‑93 (MQ=255)
                                                                               aGCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATACCGGCTCTCAGCGCGGTGAAGAGAtt          >  1:3750984/1‑93 (MQ=255)
                                                                                gCGTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATACCGGCTCTCAGCGCGGTGAAGAGATTc         <  2:2957021/93‑1 (MQ=255)
                                                                                  gtggACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATACCGGCTCTCAGCGCGGTGAAGAGATTCTg       <  2:297548/90‑1 (MQ=255)
                                                                                  gTTGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATACCGGCTCTCAGCGCGGTGAAGAGATTCTg       >  1:287251/1‑93 (MQ=255)
                                                                                   ttGACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATACCGGCTCTCAGCGCGGTGAAGAGATTCTgg      <  2:793006/93‑1 (MQ=255)
                                                                                     gACGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATACCGGCTCTCAGCGCGGTGAAGAGATTCTGGCg    <  1:1623381/93‑1 (MQ=255)
                                                                                      aCGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATACCGGCTCTCAGCGCGGTGAAGAGATTCTGGCGa   >  2:1810006/1‑93 (MQ=255)
                                                                                       cGCCCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATACCGGCTCTCAGCGCGGTGAAGAGATTCTGGCGaa  >  2:3073690/1‑93 (MQ=255)
                                                                                          |                                                                                         
GACCGGCGGTTTCGCTTACGTTCTCGATGAAAGCGGCGATTTCCGCAAACGCGTTAACCCGGAACTGGTCGAGGTCTTAAGCGTTGACGCTCTGGCGATCCATGAAGAGCATCTGCGTGGTCTTATCACCGAGCATGTGCAGCATACCGGCTCTCAGCGCGGTGAAGAGATTCTGGCGAA  >  NC_000913/3356526‑3356705

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