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,9890TC96.7% 273.1 / ‑1.4 92M136V (ATG→GTG) glpTsn‑glycerol‑3‑phosphate permease
Reads supporting (aligned to +/- strand):  ref base T (0/3);  new base C (41/48);  total (41/51)
Fisher's exact test for biased strand distribution p-value = 2.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.93e-01

Only 100 of 102 total aligned reads displayed.

GTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACATCA  >  NC_000913/2349900‑2350078
                                                                                         |                                                                                         
gtgCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCaccccacccc                                                                                        <  1:245215/93‑1 (MQ=255)
  gcgcACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCt                                                                                      <  1:2403267/93‑1 (MQ=255)
   cgcACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTg                                                                                     <  1:568290/93‑1 (MQ=255)
    gcACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTgg                                                                                    <  2:4012100/93‑1 (MQ=255)
    gcACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTgg                                                                                    <  1:701989/93‑1 (MQ=255)
    cctcaGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTgg                                                                                    <  2:4012101/90‑1 (MQ=255)
     cacaGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATCGTACGACCACACGGCGGCCACCCCACCCCCTGGa                                                                                   >  2:3173645/1‑93 (MQ=255)
     cacaGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGa                                                                                   >  2:140399/1‑93 (MQ=255)
      acaGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaa                                                                                  >  2:3966244/1‑93 (MQ=255)
        aGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATATTACGACCACACGGCGGCCACCCCACCCCCTGGaccc                                                                                <  2:4724705/93‑1 (MQ=255)
        aGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaacc                                                                                <  1:3555160/93‑1 (MQ=255)
        aGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACAACCACACGGCGGCCACCCCACCCCCTGTaacc                                                                                >  2:1609351/1‑93 (MQ=255)
        aGTTCCACACTGACACAATGCCGCCACGTTCTGTCTGCGACCACCAGTGCACCATAGTAAGACCACACGGCGGCCACCCCACCCCCTCGaaca                                                                                >  2:1609350/1‑92 (MQ=255)
         gTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaacca                                                                               >  1:2846031/1‑93 (MQ=255)
          ttCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaaccaa                                                                              >  1:2786213/1‑93 (MQ=255)
           tCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGaaccaac                                                                             <  1:4505283/93‑1 (MQ=255)
           tCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaaccaac                                                                             <  1:163341/93‑1 (MQ=255)
           tCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaaccaac                                                                             <  1:2076314/93‑1 (MQ=255)
           tCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaaccaac                                                                             <  1:4167550/93‑1 (MQ=255)
           tCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaaccaac                                                                             <  1:3657040/93‑1 (MQ=255)
            ccACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaaccaacc                                                                            >  2:3147025/1‑93 (MQ=255)
            ccACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGaaccaacc                                                                            >  2:4445890/1‑93 (MQ=255)
             cacaCTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCg                                                                           >  2:2785098/1‑93 (MQ=255)
             cacaCTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCg                                                                           >  2:2785097/1‑93 (MQ=255)
             cacaCTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCg                                                                           <  1:2717054/93‑1 (MQ=255)
                acTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCag                                                                        <  2:2211268/93‑1 (MQ=255)
                 cTGACACAATGCCGCCTCGTTCTTTCTGCGACCACCAGTGCACCATAGTACGCCCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCaga                                                                       <  2:2052414/93‑1 (MQ=255)
                 cTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCaga                                                                       >  2:4451438/1‑93 (MQ=255)
                  tGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCagag                                                                      <  2:3482719/93‑1 (MQ=255)
                  cgACACACTGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCagag                                                                      <  1:2001376/92‑1 (MQ=255)
                   gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAgg                                                                     >  2:2918886/1‑93 (MQ=255)
                   gACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAgg                                                                     <  2:4468791/93‑1 (MQ=255)
                    acacAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGTCCACACGGCGGGCACCCCACCCCCTGGAGCCAACCGCAGAGGa                                                                    >  2:3999286/1‑93 (MQ=255)
                       caATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGaaca                                                                 <  1:4055752/93‑1 (MQ=255)
                        aaTGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGaacaa                                                                <  2:3119481/93‑1 (MQ=255)
                        aaTGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGaacaa                                                                <  1:1191258/93‑1 (MQ=255)
                         aTGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGaacaac                                                               >  1:702/1‑93 (MQ=255)
                          tGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAacaaca                                                              <  2:3361950/93‑1 (MQ=255)
                          agccgccCCGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAacaaca                                                              <  1:2533412/92‑1 (MQ=255)
                           gccgccACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAg                                                             >  2:2693309/1‑93 (MQ=255)
                           gccgccACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAg                                                             <  1:4123762/93‑1 (MQ=255)
                            ccgccACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGt                                                            <  1:2398869/93‑1 (MQ=255)
                             cgccCCGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTa                                                           <  2:2088840/93‑1 (MQ=255)
                             cgccACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCCCCCCCTGGAACCAACCGCAGAGGAACAACAGTa                                                           <  1:3496573/93‑1 (MQ=255)
                             cgccACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTa                                                           <  1:797612/93‑1 (MQ=255)
                             cgccACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTa                                                           >  1:4576877/1‑93 (MQ=255)
                             cgccACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTa                                                           >  1:794222/1‑93 (MQ=255)
                              gccACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTAc                                                          >  1:3328362/1‑93 (MQ=255)
                               ccACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACa                                                         >  1:605596/1‑93 (MQ=255)
                                   gTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCCCCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACa                                                     <  2:2884869/93‑1 (MQ=255)
                                   gTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACa                                                     <  1:3642239/93‑1 (MQ=255)
                                    ttctttctGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACac                                                    >  1:1213642/1‑92 (MQ=255)
                                    ttctttctGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACAt                                                    >  2:1782126/1‑93 (MQ=255)
                                      ctttctGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCa                                                  <  2:1762331/93‑1 (MQ=255)
                                        ttctGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACGGTACAAACATCAcc                                                <  2:1576562/93‑1 (MQ=255)
                                         tctGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCg                                               >  1:4034362/1‑93 (MQ=255)
                                           tGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCa                                             >  2:2993084/1‑93 (MQ=255)
                                             cGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAAt                                           <  2:3777638/93‑1 (MQ=255)
                                              gACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATg                                          >  1:3037497/1‑93 (MQ=255)
                                              gACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATg                                          <  2:3600278/93‑1 (MQ=255)
                                               accaccAGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGc                                         >  2:1634956/1‑93 (MQ=255)
                                                 caccaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTc                                       <  2:4043672/93‑1 (MQ=255)
                                                 caccaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTc                                       >  1:3016525/1‑93 (MQ=255)
                                                  cccaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCg                                      <  2:1331511/92‑1 (MQ=255)
                                                  accaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCg                                      >  2:4449900/1‑93 (MQ=255)
                                                   ccaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACGGTACAAACATCACCGCAATGCTCGa                                     <  1:3989620/93‑1 (MQ=255)
                                                   ccaGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGa                                     >  1:617919/1‑93 (MQ=255)
                                                    caGTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGAc                                    <  1:3720404/93‑1 (MQ=255)
                                                     aGTGCCCCATAGTACGACCTCACGGCGGCCCCCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACg                                   <  1:1053662/93‑1 (MQ=255)
                                                      gTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACgc                                  >  2:756399/1‑92 (MQ=255)
                                                      gTGCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGt                                  <  1:1806028/93‑1 (MQ=255)
                                                       ggCACCATAGTACGACCACACGGCGGTCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTc                                 <  1:2566135/92‑1 (MQ=255)
                                                        gCACCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCg                                >  1:1330551/1‑93 (MQ=255)
                                                          aCCATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGcc                              <  1:2037548/93‑1 (MQ=255)
                                                            cATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCa                            <  2:3585953/93‑1 (MQ=255)
                                                            cATAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCa                            <  1:2621701/93‑1 (MQ=255)
                                                            cATAGTACCACCCCACGGCGGCCACCCCCCCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCa                            <  1:2621702/93‑1 (MQ=255)
                                                             aTAGTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCAt                           >  1:94299/1‑93 (MQ=255)
                                                                gTACGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGc                        >  1:2016470/1‑93 (MQ=255)
                                                                   cGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGcaca                     >  2:1947514/1‑93 (MQ=255)
                                                                   cGACCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGcaaa                     >  2:3584943/1‑91 (MQ=255)
                                                                     aCCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACaaa                   <  1:2638661/93‑1 (MQ=255)
                                                                     aCCACACGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACaaa                   >  2:2841016/1‑93 (MQ=255)
                                                                       cacaCGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGc                 >  1:106464/1‑93 (MQ=255)
                                                                        acacGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGcc                >  1:1345527/1‑93 (MQ=255)
                                                                         cacGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGccc               <  1:2908388/93‑1 (MQ=255)
                                                                         cacGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGccc               <  1:1138387/93‑1 (MQ=255)
                                                                          acGGCGGCCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCa              >  2:168438/1‑93 (MQ=255)
                                                                            ggcggcCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATa            <  2:909352/93‑1 (MQ=255)
                                                                            ggcggcCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATa            <  2:1728443/93‑1 (MQ=255)
                                                                             gcggcCACCCCACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATaa           <  1:259638/93‑1 (MQ=255)
                                                                                  caccccacccCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAaacaac      <  2:428050/93‑1 (MQ=255)
                                                                                  cACCCCATCCGCTGGAACCAACCGCAAAAGAACAACAGTACAAACATCACCCCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAaacaac      <  2:1578314/93‑1 (MQ=255)
                                                                                   accccaccccCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAacaaca     >  2:1283587/1‑93 (MQ=255)
                                                                                    ccccACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACAt    <  1:4743703/93‑1 (MQ=255)
                                                                                     cccACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACATc   >  2:2892827/1‑93 (MQ=255)
                                                                                     cccACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACATc   >  2:2388388/1‑93 (MQ=255)
                                                                                     cccACCCCCTGGAACCAACCGCAGAGGAACAACAGGACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACATc   >  2:3700819/1‑93 (MQ=255)
                                                                                      ccACCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACATCa  >  2:3473540/1‑93 (MQ=255)
                                                                                      ccACCCCATGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATCCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACATCa  <  2:2705779/93‑1 (MQ=255)
                                                                                         |                                                                                         
GTGCGCACAGTTCCACACTGACACAATGCCGCCACGTTCTTTCTGCGACCACCAGTGCACCATAGTACGACCACACGGCGGCCACCCCATCCCCTGGAACCAACCGCAGAGGAACAACAGTACAAACATCACCGCAATGCTCGACGTCGCCCATGGCACAAAGCCCATAAACAACATCA  >  NC_000913/2349900‑2350078

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