Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,850,058 G→A G118G (GGC→GGT ilvB ← acetolactate synthase I large subunit

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,850,0580GA99.6% 905.4 / NA 249G118G (GGC→GGTilvBacetolactate synthase I large subunit
Reads supporting (aligned to +/- strand):  ref base G (0/0);  major base A (142/106);  minor base T (0/1);  total (142/107)
Fisher's exact test for biased strand distribution p-value = 4.30e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 262 total aligned reads displayed.

CGGGAGTTCTTCGATATGTCTGACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGT  >  NC_000913/3849968‑3850147
                                                                                          |                                                                                         
cGGGAGTTCTTCGATATGTCTGACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTAcc                                                                                         >  2:2184444/1‑93 (MQ=255)
 gggAGTTCTTCGATATGTCTGACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCg                                                                                        >  1:1049402/1‑93 (MQ=255)
 gggAGTTCTTCGATATGTCTGACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCg                                                                                        >  2:1976509/1‑93 (MQ=255)
  ggAGTTCTTCGATATGTCTGACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGa                                                                                       >  1:1227410/1‑93 (MQ=255)
  ggAGTTCTTCGATATGTCTGACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGa                                                                                       >  2:455893/1‑93 (MQ=255)
  ggAGTTCTTCGATATGTCTGACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGa                                                                                       >  1:1857009/1‑93 (MQ=255)
  ggAGTTCTTCGATATGTCTGACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGa                                                                                       >  2:180884/1‑93 (MQ=255)
   gAGTTCTTCGATATGTCTGACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGat                                                                                      >  1:1607972/1‑93 (MQ=255)
   gAGTTCTTCGATATGTCTGACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGat                                                                                      >  2:5494422/1‑93 (MQ=255)
      ttcttcGATATGTCTGACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCTACTTCCTGGAAGGCGTCGGTACCGatcat                                                                                   >  2:2275696/1‑93 (MQ=255)
      ttcttcGATATGTCTGACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGatcat                                                                                   <  2:2440701/93‑1 (MQ=255)
       tcttcGATATGTCTGACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGatcatc                                                                                  >  1:3551412/1‑93 (MQ=255)
       tcttcGATATGTCTGACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGatcatc                                                                                  >  2:3143286/1‑93 (MQ=255)
       tcttcGATATGTCTGACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGatcatc                                                                                  >  2:5808924/1‑93 (MQ=255)
           cGATATGTCTGACCAGATAGTTGTGTTTGGTGATGGGCATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAgg                                                                              <  1:914337/93‑1 (MQ=255)
            gATATGTCTGACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGc                                                                             >  1:5848899/1‑93 (MQ=255)
             atatGTCTGACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCg                                                                            <  2:1722447/93‑1 (MQ=255)
             atatGTCTGACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCg                                                                            >  1:4545525/1‑93 (MQ=255)
              tatGTCTGACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCgg                                                                           >  1:2257440/1‑93 (MQ=255)
              tatGTCTGACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCgg                                                                           >  2:3328365/1‑93 (MQ=255)
              tatGTCTGACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCgg                                                                           <  2:414041/93‑1 (MQ=255)
              tatGTCTGACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCgg                                                                           <  1:4610288/93‑1 (MQ=255)
               atGTCTGACCAGATAGTTGTGTTTGGTGATGGGGATAGCGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCTAGGCggg                                                                          >  1:5987217/1‑93 (MQ=255)
               atGTCTGACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCggg                                                                          >  2:1194967/1‑93 (MQ=255)
               atGTCTGACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCggg                                                                          >  2:45130/1‑93 (MQ=255)
               atGTCTGACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCggg                                                                          >  2:5331596/1‑93 (MQ=255)
                tGTCTGACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGa                                                                         >  1:4361627/1‑93 (MQ=255)
                  tCTGACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAAc                                                                       >  2:1421967/1‑93 (MQ=255)
                   cTGACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAAcc                                                                      <  2:1674708/93‑1 (MQ=255)
                   cTGACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGGAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGCAcc                                                                      >  1:1785077/1‑93 (MQ=255)
                     gACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTg                                                                    <  2:3365921/93‑1 (MQ=255)
                     gACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTg                                                                    <  1:4059197/93‑1 (MQ=255)
                      aCCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGa                                                                   <  1:5033331/93‑1 (MQ=255)
                       ccAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGAc                                                                  <  2:5715047/93‑1 (MQ=255)
                        ctgATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGAcc                                                                 >  2:3971854/3‑93 (MQ=255)
                        cAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGAcc                                                                 >  1:5423778/1‑93 (MQ=255)
                         aGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCa                                                                <  1:328528/93‑1 (MQ=255)
                         aGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCa                                                                >  1:2193256/1‑93 (MQ=255)
                           aTAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAgg                                                              >  2:582804/1‑92 (MQ=255)
                            tAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTg                                                             <  1:3043834/93‑1 (MQ=255)
                            tAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTg                                                             >  1:5293094/1‑93 (MQ=255)
                            tAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTg                                                             >  2:4672580/1‑93 (MQ=255)
                              gTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGAt                                                           <  1:3815265/93‑1 (MQ=255)
                                tgtgTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGc                                                         <  2:4267215/93‑1 (MQ=255)
                                tgtgTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGc                                                         >  2:5287018/1‑93 (MQ=255)
                                  tgtTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCaa                                                       >  1:37049/1‑93 (MQ=255)
                                     ttGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATc                                                    <  2:92989/93‑1 (MQ=255)
                                     ttGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATc                                                    <  1:4174318/93‑1 (MQ=255)
                                       ggTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAg                                                  >  1:3212508/1‑93 (MQ=255)
                                        gTGATGGGGATAGAGATGCCGTAGGTGTCCCCTTCCTGGAATGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGc                                                 <  1:466895/93‑1 (MQ=255)
                                          gATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGCACCTGACCAGTGATGCAAATCAGCgg                                               >  2:3166734/1‑93 (MQ=255)
                                           aTGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCggg                                              >  1:2584797/1‑93 (MQ=255)
                                            tGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGa                                             <  1:3151033/93‑1 (MQ=255)
                                            tGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGa                                             <  2:1062064/93‑1 (MQ=255)
                                             ggggATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGAt                                            <  1:803598/93‑1 (MQ=255)
                                             ggggATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGAt                                            <  2:5065345/93‑1 (MQ=255)
                                              gggATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATg                                           >  2:2853991/1‑93 (MQ=255)
                                              gggATAGAGATGCCGTAGGGGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGGCCAGGGATGCAAATCAGCGGGATg                                           >  2:1572787/1‑93 (MQ=255)
                                                gATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGa                                         <  1:6283151/93‑1 (MQ=255)
                                                gATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGa                                         <  1:5071432/93‑1 (MQ=255)
                                                gATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGa                                         <  2:2963519/93‑1 (MQ=255)
                                                  tAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGt                                       >  2:4610173/1‑93 (MQ=255)
                                                  tAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGt                                       <  1:64944/93‑1 (MQ=255)
                                                  tAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGt                                       <  1:242682/93‑1 (MQ=255)
                                                    gagaTGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTcc                                     >  2:3867735/1‑93 (MQ=255)
                                                     agaTGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCa                                    <  1:2044822/93‑1 (MQ=255)
                                                          ccGTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCcgc                               <  1:1260779/93‑1 (MQ=255)
                                                          ccGTAGGTGCCCACTGCCGGGAAGGCGTCGGTACCGATCATCGAGGCGGGACCCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCcgc                               <  2:5392920/93‑1 (MQ=255)
                                                            gTAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCgcgc                             >  2:5627345/1‑93 (MQ=255)
                                                             tAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGTCCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCa                            <  1:4651649/93‑1 (MQ=255)
                                                             tAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCa                            <  2:701883/93‑1 (MQ=255)
                                                             tAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCa                            <  2:2112893/93‑1 (MQ=255)
                                                             tAGGTGTCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGCGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCa                            >  2:4018594/1‑93 (MQ=255)
                                                                  gtCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGc                       >  2:6037122/1‑93 (MQ=255)
                                                                  gtCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGc                       >  2:4801439/1‑93 (MQ=255)
                                                                   tCCACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCCGCa                      >  2:4946160/1‑93 (MQ=255)
                                                                    ccACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCaa                     <  2:4111909/93‑1 (MQ=255)
                                                                    ccACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCaa                     <  2:4995955/93‑1 (MQ=255)
                                                                    ccACTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCaa                     <  1:5254486/93‑1 (MQ=255)
                                                                       cTTCCTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATgg                  >  1:274283/1‑93 (MQ=255)
                                                                          ccTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCgg               <  2:3434408/93‑1 (MQ=255)
                                                                          ccTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCgg               >  1:2977429/1‑93 (MQ=255)
                                                                           cTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGt              >  2:1563470/1‑93 (MQ=255)
                                                                           cTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGt              >  1:1455326/1‑93 (MQ=255)
                                                                           cTGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGt              >  2:2653737/1‑93 (MQ=255)
                                                                            tGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTc             >  1:490759/1‑93 (MQ=255)
                                                                            tGGAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAAAGGCGGTc             >  1:112684/1‑93 (MQ=255)
                                                                            tGGAAGGCATCGGTACCGATCGTCGAGGCGATAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGGGCATCGGCAATGGCGGtt             <  2:4088821/93‑2 (MQ=255)
                                                                             ggAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCa            >  1:1689659/1‑93 (MQ=255)
                                                                              gAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCAc           >  1:6136681/1‑93 (MQ=255)
                                                                              gAAGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCAc           >  1:2388855/1‑93 (MQ=255)
                                                                               aaGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCAcc          <  2:797488/93‑1 (MQ=255)
                                                                                aGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCa         >  2:4731940/1‑93 (MQ=255)
                                                                                aGGCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCa         <  2:4415504/93‑1 (MQ=255)
                                                                                 ggCGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGGATCGGCAATGGCGGTGCCCAg        >  2:5183470/1‑93 (MQ=255)
                                                                                   cGTCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGt      <  1:3505025/93‑1 (MQ=255)
                                                                                     tCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAgggtc    >  2:3547715/1‑90 (MQ=255)
                                                                                     tCGGTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTa    <  1:973502/93‑1 (MQ=255)
                                                                                       ggTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGt  <  1:1923813/93‑1 (MQ=255)
                                                                                       ggTACCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGGTAGt  >  2:3735403/1‑93 (MQ=255)
                                                                                          |                                                                                         
CGGGAGTTCTTCGATATGTCTGACCAGATAGTTGTGTTTGGTGATGGGGATAGAGATGCCGTAGGTGTCCACTTCCTGGAAGGCGTCGGTGCCGATCATCGAGGCGGGAACCTGACCAGTGATGCAAATCAGCGGGATGGAGTCCAGCCGCGCATCGGCAATGGCGGTCACCAGGTTAGT  >  NC_000913/3849968‑3850147

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