Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,970,692 T→C C187R (TGC→CGC)  rffG → dTDP‑glucose 4,6‑dehydratase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,970,6920TC95.9% 566.0 / 8.7 195C187R (TGC→CGC) rffGdTDP‑glucose 4,6‑dehydratase
Reads supporting (aligned to +/- strand):  ref base T (6/2);  new base C (110/77);  total (116/79)
Fisher's exact test for biased strand distribution p-value = 4.77e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.20e-01

Only 100 of 230 total aligned reads displayed.

AGCAGCCCTTATTCCGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCGGTA  >  NC_000913/3970602‑3970781
                                                                                          |                                                                                         
 gcagcCCTTATTCCGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCt                                                                                        <  2:1003571/93‑1 (MQ=255)
  cagcCCTTATTCCGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTc                                                                                       <  1:3403010/93‑1 (MQ=255)
   agcCCTTATTCCGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCg                                                                                      >  2:63377/1‑93 (MQ=255)
      ccTTATTCCGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAAt                                                                                   >  1:462673/1‑93 (MQ=255)
      ccTTATTCCGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAAt                                                                                   >  2:2266094/1‑93 (MQ=255)
       cTTATTCCGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATa                                                                                  >  1:18365/1‑93 (MQ=255)
        ttATTCCGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATaa                                                                                 >  1:1123680/1‑93 (MQ=255)
         tATTCCGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAAc                                                                                <  2:4165399/93‑1 (MQ=255)
         tATTCCGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAAc                                                                                >  1:3451457/1‑93 (MQ=255)
          aTTCCGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACt                                                                               >  1:2342585/1‑93 (MQ=255)
           ttCCGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCATCAACCGCTCGAATAACTa                                                                              <  1:2614031/93‑1 (MQ=255)
           ttCCGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTa                                                                              <  1:1018285/93‑1 (MQ=255)
           ttCCGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTa                                                                              >  2:3896164/1‑93 (MQ=255)
           ttCCGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTa                                                                              <  2:427044/93‑1 (MQ=255)
           ttCCGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTa                                                                              <  1:4418874/93‑1 (MQ=255)
           ttCCGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTa                                                                              >  2:1225791/1‑93 (MQ=255)
            tCCGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTAc                                                                             <  1:1525504/93‑1 (MQ=255)
             ccGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACg                                                                            >  2:456709/1‑93 (MQ=255)
              cgcgTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACgg                                                                           >  2:2589018/1‑93 (MQ=255)
                cgTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGcc                                                                         >  1:330247/1‑93 (MQ=255)
                 gTCAAAAGCCAGCAGCGACCATCTGGTGTGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGccc                                                                        <  1:658593/93‑1 (MQ=255)
                 gTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGccc                                                                        >  2:3041767/1‑93 (MQ=255)
                  tCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCt                                                                       <  1:2905777/93‑1 (MQ=255)
                  tCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCt                                                                       >  1:693763/1‑93 (MQ=255)
                   cAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTGTCACCAACCGCTCGAATAACTACGGCCCtt                                                                      <  1:2420032/93‑1 (MQ=255)
                    aaaaGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTa                                                                     >  2:2494274/1‑93 (MQ=255)
                    aaaaGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTa                                                                     >  1:3253696/1‑93 (MQ=255)
                    aaaaGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTa                                                                     >  2:4641603/1‑93 (MQ=255)
                       aGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCa                                                                  >  2:1673938/1‑93 (MQ=255)
                         ccAGCAGCGCCCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTGTCACCAGCCGCTCGAATAACTACGGCCCTTACctct                                                                <  2:2313145/93‑1 (MQ=255)
                           agcagcGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACttt                                                              >  1:433584/1‑93 (MQ=255)
                           agcagcGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACttt                                                              >  2:742970/1‑93 (MQ=255)
                            gcagcGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTc                                                             <  2:3095632/93‑1 (MQ=255)
                            gcagcGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTc                                                             >  1:1143947/1‑93 (MQ=255)
                            gcagcGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTc                                                             <  1:3640802/93‑1 (MQ=255)
                            gcagcGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTc                                                             <  2:3769767/93‑1 (MQ=255)
                            gcagcGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTc                                                             <  1:195512/93‑1 (MQ=255)
                             cagcGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTcc                                                            <  2:2935788/93‑1 (MQ=255)
                              agcGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCg                                                           <  1:3195891/93‑1 (MQ=255)
                               gcGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCgg                                                          <  1:3878509/93‑1 (MQ=255)
                                cGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGa                                                         >  2:1470642/1‑93 (MQ=255)
                                cGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGa                                                         <  2:3473506/93‑1 (MQ=255)
                                 gACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGaa                                                        <  2:3272967/93‑1 (MQ=255)
                                 gACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGaa                                                        >  1:1711628/1‑93 (MQ=255)
                                 gACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGaa                                                        >  2:1119933/1‑93 (MQ=255)
                                  aCCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGaaa                                                       >  2:914301/1‑93 (MQ=255)
                                   ccATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGaaaa                                                      >  2:2577496/1‑93 (MQ=255)
                                   ccATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGaaaa                                                      <  2:129720/93‑1 (MQ=255)
                                    cATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGaaaaa                                                     >  1:2503136/1‑93 (MQ=255)
                                     aTCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAAc                                                    <  1:2289398/93‑1 (MQ=255)
                                     aTCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAAc                                                    >  2:1075775/1‑93 (MQ=255)
                                      tCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACt                                                   >  2:4326842/1‑93 (MQ=255)
                                      tCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACt                                                   >  1:1199717/1‑93 (MQ=255)
                                       cTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTg                                                  <  1:4362268/93‑1 (MQ=255)
                                        tGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGa                                                 <  1:3671963/93‑1 (MQ=255)
                                        tGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGa                                                 <  1:3820693/93‑1 (MQ=255)
                                          gTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATc                                               <  2:1632692/93‑1 (MQ=255)
                                            gcgcgcCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATccc                                             <  1:958533/93‑1 (MQ=255)
                                                  cTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCcgat                                       >  2:1204036/1‑93 (MQ=255)
                                                   tGGCTGCGCACTTATGCTCTGCGGACGTTTATCACCAACCTCTCGAATAACTCCGGCCTTTCCCACTTCCCGGAAAAACTGATCCCGCcgatg                                      <  2:1817830/93‑1 (MQ=255)
                                                    ggCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCcgatga                                     <  1:271018/93‑1 (MQ=255)
                                                     gCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCCgatgat                                    >  1:739075/1‑93 (MQ=255)
                                                      cTGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCCGATGATc                                   >  2:2139651/1‑93 (MQ=255)
                                                       tGCGGACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCCGATGATcc                                  >  1:3769272/1‑93 (MQ=255)
                                                        gCGGACCTATGGTCTGCCGACGCTTATCACCACCCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCCGATGATCCt                                 <  1:359423/93‑1 (MQ=255)
                                                         cGGCCCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTCCCGAAAAAACTGATCCCGCCGATGATCCTc                                <  2:3678423/93‑1 (MQ=255)
                                                           gACCTATGGTCTGCCGACGTTTATCACCACCCGCTCAAATAACTACGGCCCTTCCCACTTCCCGGAAAAACTGATCCCGCCGATGATCCTCaa                              <  2:2828241/93‑1 (MQ=255)
                                                           gACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCCGATGATCCTCaa                              >  1:2286651/1‑93 (MQ=255)
                                                           gACCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCCGATGATCCTCaa                              >  2:3356984/1‑93 (MQ=255)
                                                           gACCTATGGTCGGCCGACGCTTATCCCCACCCGCTCGAATAACTACGGCCCTTACCACTTCCCGGAAAAACTGATCCCGCCGATGATCCTCaa                              <  1:143879/93‑1 (MQ=255)
                                                            aCCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCCGATGATCCTCAAc                             >  2:1556904/1‑93 (MQ=255)
                                                            aCCTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCCGATGATCCTCAAc                             >  1:2192986/1‑93 (MQ=255)
                                                             ccTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCCGATGATCCTCAAcg                            >  1:1894550/1‑93 (MQ=255)
                                                             ccTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCCGATGATCCTCAAcg                            >  2:829830/1‑93 (MQ=255)
                                                             ccTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCCGATGATCCTCAAcg                            >  2:1866930/1‑93 (MQ=255)
                                                             ccTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCCGATGATCCACAAcg                            >  1:2440098/1‑93 (MQ=255)
                                                             ccTAGGGTCTGCCGACGTTTATCCCCACCCGCTCGAATAACTACGGCCCTTCCCCCTTCCGGAAAAAACTGATCCCGCCGATGATCCTCAAcg                            <  2:3378667/93‑1 (MQ=255)
                                                              cTATGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCCGATGATCCTCAAcgc                           >  1:3791781/1‑93 (MQ=255)
                                                                aTGGTCTGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCCGATGATCCTCAACgcgc                         <  1:3908398/93‑1 (MQ=255)
                                                                  ggTCTGCCGACTCTTATCCCAAACCGCCCGAATAACTACGGCCCTTACCACTTCCCGGAAAAACTGATCCCGCCGATGATCCTCAACGCGCTg                       <  1:537774/93‑1 (MQ=255)
                                                                      tGCCGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCCGATGATCCTCAACGCGCTGGCgg                   >  1:1341721/1‑93 (MQ=255)
                                                                         cGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCCGATGATCCTCAACGCGCTGGCgggca                >  1:3088979/1‑91 (MQ=255)
                                                                         cGACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCCGATGATCCTCAACGCGCTGGCGGGTa                >  1:4666345/1‑93 (MQ=255)
                                                                          gACGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCCGATGATCCTCAACGCGCTGGCGGGGaa               >  2:3612823/1‑93 (MQ=255)
                                                                            cctttATCACCACCCGCTCAAATACCTACGGCCCTTCCCTCTTTCGGGAAAAACTGACCCCGCCGATGATCCTCAACGCGCTGGCGGGTAAAt             <  2:4431929/90‑1 (MQ=255)
                                                                            cGCTTATCACCACCCCCCGAAATAACTGCGCCCTTTCCCTCTTACCGAAAAACCTGATCCCGCCGATGATCCTCAACGCGCTGGCGGGTAAAt             <  2:2700577/93‑1 (MQ=255)
                                                                            cGCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCCGATGATCCTCAACGCGCTGGCGGGTAAAt             >  2:1509626/1‑93 (MQ=255)
                                                                             gCTTATCACCAACCGCTCGAATAACTACGTCCCTTACCACTTTCCGGAAAAACTGACCCCGCCGATGATCCTCAACGCGCTGGCGGGTAAATc            <  1:1286381/93‑1 (MQ=255)
                                                                             gCTTATCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCCGATGATCCTCAACGCGCTGGCGGGTAAATc            >  1:1048514/1‑93 (MQ=255)
                                                                              cTTATCCCCAACCGCTCAAATAACTCGGCCCCTTCCCACTTCCCGAAAAAACTGATCCCGCCGATGATCCTCAACGCGCTGGCGGGTAAATCg           <  2:4490024/93‑1 (MQ=255)
                                                                               ttATCCCCAACCGCTCGAATAACTACGGCCCTTACCACTTCCCGGAAAAACTGACCCCGCCGATGATCCTCAACGCGCTGGCGGGTAAATCGc          <  2:2032235/93‑1 (MQ=255)
                                                                               taaccACCAACCGCTCGAATAACAACGCCCCTTCCCTCTTTCGGAAAAACCGGACCCCGCCGATGATCCTCAACGCGCTGGCGGGTAAATCGc          <  2:4139289/89‑1 (MQ=255)
                                                                                tATCACCAACCGCTCTAATAACTACGGCCCTTACCACTTTCCGAAAAAACGGATCCCGCCGATGATCCTCAACGCGCTGGCGGGTAAATCGCt         <  1:1407164/93‑1 (MQ=255)
                                                                                  tCACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCCGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGc       >  1:3960390/1‑93 (MQ=255)
                                                                                   cACCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCCGATGATCCTCAACGCGCTGGCGGGGAAATCGCTGcc      >  1:2601994/1‑93 (MQ=255)
                                                                                    aCCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCCGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCg     >  2:2214871/1‑93 (MQ=255)
                                                                                    aCCAACCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCCGATGATCCTCAACGCGCTGGCGGGGAAATCGCTGCCg     >  1:4004236/1‑93 (MQ=255)
                                                                                       aaCTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCCGATGATCCTCAACGCGCTGGCGGGAAAATCGCTGCCGGta  >  2:3934954/1‑93 (MQ=255)
                                                                                       aaCCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCCGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCGGta  >  2:4061820/1‑93 (MQ=255)
                                                                                       aaCCGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCCGATGATCCTCAACGAGCTGGCGGGGAAATAGTTGCCGGtt  >  2:4214393/1‑92 (MQ=255)
                                                                                          |                                                                                         
AGCAGCCCTTATTCCGCGTCAAAAGCCAGCAGCGACCATCTGGTGCGCGCCTGGCTGCGGACCTATGGTCTGCCGACGCTTATCACCAACTGCTCGAATAACTACGGCCCTTACCACTTTCCGGAAAAACTGATCCCGCTGATGATCCTCAACGCGCTGGCGGGTAAATCGCTGCCGGTA  >  NC_000913/3970602‑3970781

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: