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.2% 362.1 / NA 119M136V (ATG→GTG) glpTsn‑glycerol‑3‑phosphate permease
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (58/60);  minor base A (1/0);  total (59/60)
Fisher's exact test for biased strand distribution p-value = 4.96e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 132 total aligned reads displayed.

TGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACATCA  >  NC_000913/2349901‑2350078
                                                                                        |                                                                                         
tgCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACcccc                                                                                       >  1:321568/1‑93 (MQ=255)
tgCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACcccc                                                                                       >  2:3053828/1‑93 (MQ=255)
 gcgcACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCt                                                                                      >  1:179498/1‑93 (MQ=255)
   gcACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTgg                                                                                    >  2:482565/1‑93 (MQ=255)
   gcACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTgg                                                                                    >  2:4270795/1‑93 (MQ=255)
   gcACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTgg                                                                                    >  2:4109443/1‑93 (MQ=255)
     acaGTTCCACGCTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCATGGaa                                                                                  <  2:2380985/93‑1 (MQ=255)
      caGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGGCACCCCACCCCCTGGaac                                                                                 >  1:3446512/1‑93 (MQ=255)
      caGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaac                                                                                 <  2:243853/93‑1 (MQ=255)
       aGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaacc                                                                                >  1:1904926/1‑93 (MQ=255)
        gatCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGAACGACCACACGGCGGCCCCACCACCCCCTGGaacca                                                                               >  2:3890802/3‑93 (MQ=255)
        gTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaacca                                                                               <  1:2668265/93‑1 (MQ=255)
          tCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaaccaac                                                                             >  1:1609747/1‑93 (MQ=255)
           ccACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaaccaacc                                                                            <  2:1055774/93‑1 (MQ=255)
            cacaCTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCg                                                                           >  1:4642955/1‑93 (MQ=255)
             acacTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGc                                                                          <  2:3001366/93‑1 (MQ=255)
              cacTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCCCCCCACCCCCTGGAACCAACCGCa                                                                         <  2:4077509/93‑1 (MQ=255)
               acTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCag                                                                        >  1:1942947/1‑93 (MQ=255)
                ctcacacAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCaga                                                                       <  2:1813086/90‑1 (MQ=255)
                cTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCaga                                                                       <  2:3169269/93‑1 (MQ=255)
                 tGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCagag                                                                      >  2:2818598/1‑93 (MQ=255)
                 tGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCagag                                                                      >  2:1163250/1‑93 (MQ=255)
                  gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAgg                                                                     >  2:4156437/1‑93 (MQ=255)
                   acacAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGa                                                                    <  1:2742719/93‑1 (MQ=255)
                     acaATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGatc                                                                  >  1:1037318/1‑91 (MQ=255)
                     acaATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGaac                                                                  >  1:3896466/1‑93 (MQ=255)
                      caATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGaaca                                                                 >  2:160876/1‑93 (MQ=255)
                        aTGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGaacaac                                                               <  1:3074186/93‑1 (MQ=255)
                         tGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAacaaca                                                              >  1:1206374/1‑93 (MQ=255)
                           ccgccACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCCCCCCCTGGAACCAACCGCAGAGGAACAACAGt                                                            <  2:4090892/93‑1 (MQ=255)
                           ccgccACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGt                                                            <  2:3387048/93‑1 (MQ=255)
                            cgccACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTa                                                           <  2:362791/93‑1 (MQ=255)
                             gccACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTAc                                                          >  1:2350576/1‑93 (MQ=255)
                             gccACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTAc                                                          <  2:3442447/93‑1 (MQ=255)
                             cccACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCCCCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTAc                                                          <  1:733392/92‑1 (MQ=255)
                               cACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGCCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACaa                                                        <  2:1169994/93‑1 (MQ=255)
                                aCGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACaaa                                                       >  1:623878/1‑93 (MQ=255)
                                aCGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACaaa                                                       >  1:3241057/1‑93 (MQ=255)
                                 cGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTacacac                                                      >  1:3320721/1‑93 (MQ=255)
                                 cGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAAc                                                      >  2:952177/1‑93 (MQ=255)
                                  gTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCCCCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACa                                                     <  2:3229871/93‑1 (MQ=255)
                                  gTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACa                                                     >  1:2506302/1‑93 (MQ=255)
                                  gTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACa                                                     >  2:234520/1‑93 (MQ=255)
                                  gTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACa                                                     >  1:1090931/1‑93 (MQ=255)
                                  gTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACa                                                     <  1:2182143/93‑1 (MQ=255)
                                   ttctctctGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACAt                                                    <  1:2590422/93‑1 (MQ=255)
                                    tctttctGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATc                                                   <  1:226030/93‑1 (MQ=255)
                                    tctttctGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATc                                                   >  1:3674915/1‑93 (MQ=255)
                                    tctttctGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATc                                                   >  1:4212577/1‑93 (MQ=255)
                                    tctttctGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGCACAAACATc                                                   >  2:477726/1‑93 (MQ=255)
                                     ctttctGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTCGAACCAACCGCAGAGGAACAACAGTACAAACATCa                                                  >  1:1052731/1‑93 (MQ=255)
                                       ttctGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCCCCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCAcc                                                <  2:1256878/93‑1 (MQ=255)
                                           gCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCaa                                            <  1:746763/93‑1 (MQ=255)
                                            cGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAAt                                           <  1:2428144/93‑1 (MQ=255)
                                              accaccAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGc                                         >  2:4262254/1‑93 (MQ=255)
                                               ccaccaGTGCACCATAGTACGACCACACGGCGGCCACCCCAACCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCt                                        >  2:2409081/1‑93 (MQ=255)
                                                caccaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTc                                       >  2:4334301/1‑93 (MQ=255)
                                                caccaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTc                                       >  2:4196829/1‑93 (MQ=255)
                                                 cccaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCg                                      <  1:970417/92‑1 (MQ=255)
                                                 accaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCg                                      >  1:147742/1‑93 (MQ=255)
                                                  ccaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGa                                     <  1:1025776/93‑1 (MQ=255)
                                                  ccaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGa                                     <  1:2534978/93‑1 (MQ=255)
                                                  ccaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGa                                     >  2:589364/1‑93 (MQ=255)
                                                  ccaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGa                                     >  1:3656541/1‑93 (MQ=255)
                                                   caGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGAc                                    <  1:2970721/93‑1 (MQ=255)
                                                        cACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGc                               <  1:1752935/93‑1 (MQ=255)
                                                         cccATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGcc                              <  1:2746345/92‑1 (MQ=255)
                                                         aCCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGcc                              >  1:3529490/1‑93 (MQ=255)
                                                          ccATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGccc                             >  2:4369359/1‑93 (MQ=255)
                                                          ccATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGccc                             <  1:1198145/93‑1 (MQ=255)
                                                          ccATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGccc                             <  1:1226247/93‑1 (MQ=255)
                                                            aTAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCAt                           >  2:2235625/1‑93 (MQ=255)
                                                            aTAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCAt                           <  2:2101947/93‑1 (MQ=255)
                                                            aTAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCAt                           >  2:4248344/1‑93 (MQ=255)
                                                             tAGTACGACCACACGGCGGCCACCCCACCCCCTGGCACCAACCCCAGAGGAACAACAGTCCAAACATCACAGCAATGCTCGACCTCGCCCATg                          >  1:4575805/1‑93 (MQ=255)
                                                                 aCGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGcac                      <  1:3286123/93‑1 (MQ=255)
                                                                 aCGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACAGCACCGCAATGCGCGACGTCGCCCATGGcac                      <  2:96609/93‑1 (MQ=255)
                                                                    aCCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACaaa                   <  1:2115128/93‑1 (MQ=255)
                                                                       acacGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGcc                >  1:2504062/1‑93 (MQ=255)
                                                                       acacGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGcc                <  1:2998102/93‑1 (MQ=255)
                                                                        cacGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGccc               <  2:3237141/93‑1 (MQ=255)
                                                                         acGGCGGCNCNCCCACCCCCTGGANCCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCa              <  2:1563725/93‑1 (MQ=255)
                                                                         acGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCa              >  2:851400/1‑93 (MQ=255)
                                                                          cggcggCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCAt             <  1:3373998/93‑1 (MQ=255)
                                                                          cggcggCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCAt             <  1:4270795/93‑1 (MQ=255)
                                                                            gcggcCCGCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACGTCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATaa           <  1:3294118/93‑1 (MQ=255)
                                                                            gcggcCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATaa           >  1:1790520/1‑93 (MQ=255)
                                                                             cggcCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATaaa          <  1:921398/93‑1 (MQ=255)
                                                                             cgcccaccccaccCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATaaa          <  2:3980871/90‑1 (MQ=255)
                                                                                 caccccacccCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAaacaac      <  2:629242/93‑1 (MQ=255)
                                                                                 caccccacccCCTGGAACCAACCGCAGAGGAACAACAGGACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAaacaac      >  1:2247044/1‑93 (MQ=255)
                                                                                 cACCCCCCCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAaacaac      <  2:3488860/93‑1 (MQ=255)
                                                                                  cccccACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAacaaca     <  1:2942683/92‑1 (MQ=255)
                                                                                  cccccACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAacaaca     <  1:186462/92‑1 (MQ=255)
                                                                                  accccaccccCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAacaaca     >  2:3926942/1‑93 (MQ=255)
                                                                                   ccccACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACAt    >  1:4314755/1‑93 (MQ=255)
                                                                                    cccACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAGGCCCATAAACAACATc   <  1:1510245/93‑1 (MQ=255)
                                                                                    cccACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACATc   <  1:4436050/93‑1 (MQ=255)
                                                                                    cccACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACATc   <  2:944054/93‑1 (MQ=255)
                                                                                    accACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACATc   <  1:744792/92‑1 (MQ=255)
                                                                                        |                                                                                         
TGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACATCA  >  NC_000913/2349901‑2350078

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 15 ≤ ATCG/ATCG < 24 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: