Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,172,208 T→C M31V (ATG→GTG)  gatC ← PTS system galactitol‑specific enzyme IIC

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,172,2080TC99.3% 465.8 / ‑4.7 147M31V (ATG→GTG) gatCPTS system galactitol‑specific enzyme IIC
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (76/70);  total (76/71)
Fisher's exact test for biased strand distribution p-value = 4.83e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.26e-01

Only 100 of 149 total aligned reads displayed.

ATCCAGCATTAAGCCAATCACAAGGCCAATGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCATGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGAAAACA  >  NC_000913/2172122‑2172297
                                                                                      |                                                                                         
aTCCAGCATTAAGCCAATCACAAGGCCAATGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTaa                                                                                     <  2:1285903/93‑1 (MQ=255)
aTCCAGCATTAAGCCAATCACAAGGCCAATGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTaa                                                                                     <  2:528997/93‑1 (MQ=255)
  ccAGCATTAAGCCAATCACAAGGCCAATGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAata                                                                                   >  1:3294191/1‑93 (MQ=255)
  caaGCATTAAGCCAATCACAAGGCCAATGACAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAata                                                                                   <  2:2837673/91‑1 (MQ=255)
   cAGCATTAAGCCAATCACAAGGCCAATGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAatat                                                                                  >  1:2295657/1‑93 (MQ=255)
      cATTAAGCCAATCACAAGGCCAATGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATAtttt                                                                               >  1:3355731/1‑93 (MQ=255)
       aTTAAGCCAATCACAAGGCCAATGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTa                                                                              <  1:3592898/93‑1 (MQ=255)
        ttAAGCCAATCACAAGGCCAATGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAg                                                                             >  2:1514971/1‑93 (MQ=255)
         tAAGCCAATCACAAGGCCAATGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGa                                                                            <  1:3232653/93‑1 (MQ=255)
            gCCAATCACAAGGCCAATGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGaaaa                                                                         <  2:181630/93‑1 (MQ=255)
             ccAATCACAAGGCCAATGCCAACAAAGCCAATCCCGATATGCAGACCTGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGaaaaa                                                                        <  1:3041419/93‑1 (MQ=255)
              cAATCACAAGGCCAATGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTTAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGaaaaaa                                                                       >  1:4686423/1‑93 (MQ=255)
              cAATCACAAGGCCAATGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGaaaaaa                                                                       <  1:1784285/93‑1 (MQ=255)
              cAATCACAAGGCCAATGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGaaaaaa                                                                       <  1:2526593/93‑1 (MQ=255)
               aaTCACAAGGCCAATGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAaat                                                                      >  2:4402844/1‑93 (MQ=255)
                 tCACAAGGCCAATGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCTCCTGCCTTCACGCCTAATATTTTAGAAAAaataa                                                                    >  1:1053075/1‑93 (MQ=255)
                 tCACAAGGCCAATGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAaataa                                                                    <  1:797089/93‑1 (MQ=255)
                 tCACAAGGCCAATGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAaataa                                                                    >  1:252499/1‑93 (MQ=255)
                 tCACAAGGCCAATGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAaataa                                                                    <  2:2190539/93‑1 (MQ=255)
                 tCACAAGGCCAATGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAaataa                                                                    <  2:1830680/93‑1 (MQ=255)
                 acacaAGGCCAATGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAaataa                                                                    <  2:4218890/92‑1 (MQ=255)
                  cacaAGGCCAATGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAaataat                                                                   <  2:3837857/93‑1 (MQ=255)
                    caAGGCCAATGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCATGCCTAATATTTTAGAAAAAATAatga                                                                 <  1:4414158/93‑1 (MQ=255)
                    caAGGCCAATGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAatga                                                                 <  2:3261847/93‑1 (MQ=255)
                     aaGGCCAATGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAatgat                                                                <  1:1039600/93‑1 (MQ=255)
                      aGGCCAATGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAatgatg                                                               <  2:3415361/93‑1 (MQ=255)
                       ggCCAATGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAAtgatga                                                              >  1:3830841/1‑93 (MQ=255)
                       ggCCAATGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAAtgatga                                                              <  1:1974274/93‑1 (MQ=255)
                       ggCCAATGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAAtgatga                                                              <  2:2924589/93‑1 (MQ=255)
                       ggCCAATGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAAtgatga                                                              <  1:4510797/93‑1 (MQ=255)
                          cAATGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACaa                                                           >  1:1241728/1‑93 (MQ=255)
                            aTGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGACTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATc                                                         <  2:4671594/93‑1 (MQ=255)
                             tGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCg                                                        <  2:2496089/93‑1 (MQ=255)
                             tGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCg                                                        <  2:3124107/93‑1 (MQ=255)
                               ccAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGgc                                                      <  2:3609611/93‑1 (MQ=255)
                                caacaaAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGgca                                                     <  1:998036/93‑1 (MQ=255)
                                caacaaAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGgca                                                     <  2:1984474/93‑1 (MQ=255)
                                caacaaAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGgca                                                     >  2:2518912/1‑93 (MQ=255)
                                 aacaaAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGgcag                                                    >  1:3283221/1‑93 (MQ=255)
                                 aacaaAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGgcag                                                    >  2:2549294/1‑93 (MQ=255)
                                 aacaaAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGgcag                                                    <  1:2163290/93‑1 (MQ=255)
                                 aacaaAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGgcag                                                    >  1:4232224/1‑93 (MQ=255)
                                 aacaaAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGgcag                                                    >  1:4352307/1‑93 (MQ=255)
                                  acaaAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGgcagc                                                   >  1:2713126/1‑93 (MQ=255)
                                   caaAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGgcagca                                                  <  2:1071584/93‑1 (MQ=255)
                                   caaAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGgcagca                                                  <  1:618475/93‑1 (MQ=255)
                                   caaAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGgcagca                                                  >  2:3421653/1‑93 (MQ=255)
                                    aaaGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCAt                                                 >  1:4458824/1‑93 (MQ=255)
                                    aaaGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCAt                                                 >  1:3323261/1‑93 (MQ=255)
                                       gCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCAc                                              >  2:3233435/1‑93 (MQ=255)
                                       gCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCAc                                              >  2:134635/1‑93 (MQ=255)
                                       gCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCAc                                              >  2:2165417/1‑93 (MQ=255)
                                       gCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCAc                                              >  1:2892003/1‑93 (MQ=255)
                                        ccAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCAcc                                             <  2:3917529/93‑1 (MQ=255)
                                         cAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCg                                            >  2:1618634/1‑93 (MQ=255)
                                             cccGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAgg                                        <  1:221863/93‑1 (MQ=255)
                                              ccGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAggg                                       >  2:3113676/1‑93 (MQ=255)
                                              ccGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAggg                                       >  1:313937/1‑93 (MQ=255)
                                               cGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGc                                      <  1:2476603/93‑1 (MQ=255)
                                               cGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGc                                      <  1:1246325/93‑1 (MQ=255)
                                                gATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGcc                                     >  2:2806725/1‑93 (MQ=255)
                                                  tatGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCCTCACCGTAGGGCCGa                                   >  2:3812153/1‑93 (MQ=255)
                                                     gCAGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGt                                <  1:2877721/93‑1 (MQ=255)
                                                       aGACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCg                              >  1:1366752/1‑93 (MQ=255)
                                                        gACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCga                             >  2:3189818/1‑93 (MQ=255)
                                                        gACCCGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCga                             <  1:2031689/93‑1 (MQ=255)
                                                           ccGCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGaa                          <  2:2970840/93‑1 (MQ=255)
                                                             gCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAata                        <  1:116074/93‑1 (MQ=255)
                                                             gCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAata                        >  1:1067645/1‑93 (MQ=255)
                                                             gCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAata                        >  1:3776974/1‑93 (MQ=255)
                                                             gCTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAata                        >  1:115276/1‑93 (MQ=255)
                                                              cTTTAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAatat                       <  2:3056185/93‑1 (MQ=255)
                                                               tttAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAAtata                      >  2:4130765/1‑93 (MQ=255)
                                                               tttAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAAtata                      >  2:3336583/1‑93 (MQ=255)
                                                               tttAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAAtata                      <  1:1723307/93‑1 (MQ=255)
                                                               tttAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAAtata                      <  1:2589256/93‑1 (MQ=255)
                                                                ttAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATaa                     <  1:356577/93‑1 (MQ=255)
                                                                ttAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATaa                     <  1:1386620/93‑1 (MQ=255)
                                                                ttAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATaa                     <  1:2710255/93‑1 (MQ=255)
                                                                 tAAAGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAAc                    <  2:1601103/93‑1 (MQ=255)
                                                                  aaaGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACg                   >  1:3944077/1‑93 (MQ=255)
                                                                   aaGCAATCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGc                  >  2:1419052/1‑93 (MQ=255)
                                                                       aaTCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGa              >  2:679037/1‑93 (MQ=255)
                                                                       aaTCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGa              >  2:2859960/1‑93 (MQ=255)
                                                                        aTCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGAc             >  1:2488122/1‑93 (MQ=255)
                                                                        aTCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGAc             >  1:2629205/1‑93 (MQ=255)
                                                                        aTCGCCTGCCTTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGAc             >  1:894636/1‑93 (MQ=255)
                                                                            cctgcctTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCt         >  1:3047167/1‑93 (MQ=255)
                                                                            cctgcctTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCt         >  2:4413002/1‑93 (MQ=255)
                                                                             ctgcctTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTg        >  2:718233/1‑93 (MQ=255)
                                                                             ctgcctTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTg        >  1:683429/1‑93 (MQ=255)
                                                                             ctgcctTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTg        >  1:4629914/1‑93 (MQ=255)
                                                                              tgcctTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGa       >  1:4190131/1‑93 (MQ=255)
                                                                               gcctTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGaa      >  2:2052863/1‑93 (MQ=255)
                                                                               gcctTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGaa      >  1:629404/1‑93 (MQ=255)
                                                                               gcctTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGaa      >  1:3954669/1‑93 (MQ=255)
                                                                                cctTCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGaaa     >  2:3644167/1‑93 (MQ=255)
                                                                                  ttCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGAAAAc   >  2:247087/1‑93 (MQ=255)
                                                                                  ttCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGAAAAc   <  1:2367798/93‑1 (MQ=255)
                                                                                   tCACGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGAAAACa  <  2:3539150/93‑1 (MQ=255)
                                                                                      |                                                                                         
ATCCAGCATTAAGCCAATCACAAGGCCAATGCCAACAAAGCCAATCCCGATATGCAGACCCGCTTTAAAGCAATCGCCTGCCTTCATGCCTAATATTTTAGAAAAAATAATGATGACAATCGGCAGCATCACCGTAGGGCCGAGGTCGAGAATATAACGCATGACTTCTGAAAACA  >  NC_000913/2172122‑2172297

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