Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,349,989 T→C M136V (ATG→GTG)  glpT ← sn‑glycerol‑3‑phosphate permease

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,349,9890TC99.5% 661.8 / ‑5.1 209M136V (ATG→GTG) glpTsn‑glycerol‑3‑phosphate permease
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (118/90);  total (118/91)
Fisher's exact test for biased strand distribution p-value = 4.35e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.80e-01

Only 100 of 224 total aligned reads displayed.

TGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACATCA  >  NC_000913/2349899‑2350078
                                                                                          |                                                                                         
  tgCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACcccc                                                                                       >  2:337753/1‑93 (MQ=255)
  tgCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACcccc                                                                                       <  1:7303316/93‑1 (MQ=255)
   gcgcACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCt                                                                                      <  1:1344359/93‑1 (MQ=255)
      cacaGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGa                                                                                   >  1:6428463/1‑93 (MQ=255)
       acaGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaa                                                                                  >  2:5447978/1‑93 (MQ=255)
        caGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaac                                                                                 >  2:2150069/1‑93 (MQ=255)
         aGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaacc                                                                                <  2:5144389/93‑1 (MQ=255)
          gTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaacca                                                                               >  2:9495106/1‑93 (MQ=255)
          gTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaacca                                                                               >  1:8892453/1‑93 (MQ=255)
             ccACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaaccaacc                                                                            >  2:2626333/1‑93 (MQ=255)
                cacTGCCACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCa                                                                         <  1:6315670/93‑1 (MQ=255)
                cacTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCa                                                                         >  2:1629746/1‑93 (MQ=255)
                cacTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCa                                                                         <  2:10021094/93‑1 (MQ=255)
                cacTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCa                                                                         <  1:9464284/93‑1 (MQ=255)
                cacTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCa                                                                         <  1:9181707/93‑1 (MQ=255)
                cacTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCa                                                                         <  1:7155874/93‑1 (MQ=255)
                cacTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCa                                                                         >  1:3639916/1‑93 (MQ=255)
                cacTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCa                                                                         >  2:7745795/1‑93 (MQ=255)
                 acTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCag                                                                        <  2:9059181/93‑1 (MQ=255)
                 acTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCag                                                                        <  1:6304835/93‑1 (MQ=255)
                  cTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCaga                                                                       >  2:3960296/1‑93 (MQ=255)
                  cTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCaga                                                                       <  1:8465314/93‑1 (MQ=255)
                  cTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCaga                                                                       >  1:5433586/1‑93 (MQ=255)
                    gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAgg                                                                     <  2:7078568/93‑1 (MQ=255)
                    gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAgg                                                                     <  1:7493369/93‑1 (MQ=255)
                     acacAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGa                                                                    <  1:2978062/93‑1 (MQ=255)
                      cccaATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGaa                                                                   <  1:2674825/91‑1 (MQ=255)
                       acaATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCCGAGGaac                                                                  >  2:6694275/1‑93 (MQ=255)
                        caATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGaaca                                                                 >  1:8818287/1‑93 (MQ=255)
                        caATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGaaca                                                                 >  2:2759781/1‑93 (MQ=255)
                         aaTGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCCCCCCCCCCCCTGGAACCAACCGCAGAGGaacaa                                                                <  2:4707018/93‑1 (MQ=255)
                          aTGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGaacaac                                                               >  2:768164/1‑93 (MQ=255)
                          aTGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGaacaac                                                               <  2:9297034/93‑1 (MQ=255)
                           agccgccCCCTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCTCCCCCTGGAACCAACCGCAGAGGAacaaca                                                              <  2:7274200/92‑1 (MQ=255)
                            gccgccACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAg                                                             <  1:561114/93‑1 (MQ=255)
                            gccgccACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACACCCTGGAACCAACCGCAGAGGAACAACAg                                                             >  2:3862465/1‑93 (MQ=255)
                              cgccACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCGCCCCTCCCCCTGTAACCAACCGCAGAGGAACAACAGTa                                                           <  1:4796335/93‑1 (MQ=255)
                               gccACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTAc                                                          <  1:3013505/93‑1 (MQ=255)
                               gccACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTAc                                                          <  2:1724685/93‑1 (MQ=255)
                                 cACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACaa                                                        <  1:3731819/93‑1 (MQ=255)
                                  aCGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACaaa                                                       >  2:4202049/1‑93 (MQ=255)
                                  aCGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACaaa                                                       <  2:2317402/93‑1 (MQ=255)
                                   cGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAAc                                                      >  1:1985918/1‑93 (MQ=255)
                                   cGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAAc                                                      >  2:9089466/1‑93 (MQ=255)
                                    gTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACa                                                     <  1:2599423/93‑1 (MQ=255)
                                     ttctttctGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACAt                                                    <  2:435878/93‑1 (MQ=255)
                                     ttctttctGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACAt                                                    >  2:7351336/1‑93 (MQ=255)
                                      tctttctGCGACCGCCAGTGCACCATAGTACGACCACACGGCGGCCTCCCCGCCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATc                                                   <  1:9150646/93‑1 (MQ=255)
                                      tctttctGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATc                                                   <  1:5825747/93‑1 (MQ=255)
                                      tctttctGCGACCACCAGTCCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATc                                                   <  1:2915262/93‑1 (MQ=255)
                                       ctttctGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCa                                                  >  1:6876196/1‑93 (MQ=255)
                                        tttctGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCAc                                                 >  2:9072128/1‑93 (MQ=255)
                                          tctGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCg                                               >  2:4423853/1‑93 (MQ=255)
                                            tGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCa                                             >  2:7834578/1‑93 (MQ=255)
                                             gCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCaa                                            <  1:275650/93‑1 (MQ=255)
                                              cGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAAt                                           <  2:1090126/93‑1 (MQ=255)
                                              cGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAAt                                           >  1:5335924/1‑93 (MQ=255)
                                              cGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAAt                                           <  1:1230237/93‑1 (MQ=255)
                                              cGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAAt                                           <  1:6700081/93‑1 (MQ=255)
                                               gACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATg                                          <  2:8622078/93‑1 (MQ=255)
                                               gACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATg                                          <  1:3910318/93‑1 (MQ=255)
                                                 ccaccaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCt                                        >  1:1812955/1‑93 (MQ=255)
                                                   accaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCg                                      >  1:4096554/1‑93 (MQ=255)
                                                   accaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCg                                      >  2:5365109/1‑93 (MQ=255)
                                                   accaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGAAATGCTCg                                      >  1:7832748/1‑93 (MQ=255)
                                                    ccaGTGCACCATAGTACGGCCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGa                                     <  2:5251789/93‑1 (MQ=255)
                                                    ccaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCgc                                     >  2:4939233/1‑92 (MQ=255)
                                                    ccaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGa                                     >  2:1883633/1‑93 (MQ=255)
                                                    ccaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGa                                     <  2:4895432/93‑1 (MQ=255)
                                                    ccaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGa                                     <  2:4867703/93‑1 (MQ=255)
                                                    ccaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGa                                     <  1:6182979/93‑1 (MQ=255)
                                                     caGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGAc                                    <  2:2459832/93‑1 (MQ=255)
                                                      aGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACg                                   >  1:597596/1‑93 (MQ=255)
                                                       gagCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGt                                  >  2:6554556/3‑93 (MQ=255)
                                                            ccATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGccc                             >  2:9638469/1‑93 (MQ=255)
                                                             cATAGTACGCCCCCACGGCGACCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCa                            <  2:375983/93‑1 (MQ=255)
                                                             cATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCa                            <  1:5089257/93‑1 (MQ=255)
                                                             cATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCa                            <  2:3527173/93‑1 (MQ=255)
                                                             cATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCa                            <  2:1617627/93‑1 (MQ=255)
                                                             cATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCa                            <  2:251382/93‑1 (MQ=255)
                                                               tAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATg                          >  2:7272571/1‑93 (MQ=255)
                                                                aGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATgg                         >  1:4182708/1‑93 (MQ=255)
                                                                 gTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGc                        <  1:9351253/93‑1 (MQ=255)
                                                                 gTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGc                        <  1:871672/93‑1 (MQ=255)
                                                                  tACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGca                       >  1:465799/1‑93 (MQ=255)
                                                                  tACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGca                       <  1:2341766/93‑1 (MQ=255)
                                                                  tACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGca                       >  1:1635252/1‑93 (MQ=255)
                                                                      aCCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACaaa                   >  1:9921838/1‑93 (MQ=255)
                                                                      aCCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACaaa                   <  1:240247/93‑1 (MQ=255)
                                                                          cacGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGccc               <  1:4370631/93‑1 (MQ=255)
                                                                           acGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCa              >  1:4841355/1‑93 (MQ=255)
                                                                           acGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCa              >  1:2178618/1‑93 (MQ=255)
                                                                                 gcCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAaaca        <  2:5708386/93‑1 (MQ=255)
                                                                                 gcCACCCCACCCCCTGGAACCAACCACAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAaaca        >  2:8056925/1‑93 (MQ=255)
                                                                                  ccaccccaccCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAaacaa       <  2:6229248/93‑1 (MQ=255)
                                                                                    accccaccccCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAacaaca     <  1:1473701/93‑1 (MQ=255)
                                                                                     ccccACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACAt    <  2:7316879/93‑1 (MQ=255)
                                                                                     ccccACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACAt    >  2:1512350/1‑93 (MQ=255)
                                                                                     ccccACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCAATAAACAACAt    >  2:7755388/1‑93 (MQ=255)
                                                                                       ccACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACATCa  >  2:4175896/1‑93 (MQ=255)
                                                                                          |                                                                                         
TGTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACATCA  >  NC_000913/2349899‑2350078

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