Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,838,693 A→G I280T (ATC→ACC)  yicM ← putative transport protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,838,6930AG99.1% 697.4 / ‑3.0 218I280T (ATC→ACC) yicMputative transport protein
Reads supporting (aligned to +/- strand):  ref base A (1/1);  new base G (91/125);  total (92/126)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.80e-01

Only 100 of 223 total aligned reads displayed.

GTTAAGCCATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGATGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGTGAAGGC  >  NC_000913/3838603‑3838782
                                                                                          |                                                                                         
gTTAAGCCATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTg                                                                                         >  2:2429267/1‑93 (MQ=255)
gTTAAGCCATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTg                                                                                         >  1:3329887/1‑93 (MQ=255)
 ttAAGCCATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTgg                                                                                        <  1:1682918/93‑1 (MQ=255)
  tAAGCCATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGc                                                                                       <  2:2315128/93‑1 (MQ=255)
     gCCATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGat                                                                                    <  2:2501390/93‑1 (MQ=255)
      ccATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGatc                                                                                   <  2:2991926/93‑1 (MQ=255)
       cATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGatca                                                                                  >  2:1461794/1‑93 (MQ=255)
         tCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGatcatc                                                                                <  2:225272/93‑1 (MQ=255)
         tCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGatcatc                                                                                >  1:2396185/1‑93 (MQ=255)
         tCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGGGGCGatcatc                                                                                >  1:4085470/1‑93 (MQ=255)
          ccACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATcc                                                                               <  2:1026370/93‑1 (MQ=255)
           cACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATccc                                                                              <  1:3883942/93‑1 (MQ=255)
            ccGCCCAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCt                                                                             <  2:3972849/92‑1 (MQ=255)
              gCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGc                                                                           >  2:1221050/1‑93 (MQ=255)
               ccGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGcc                                                                          <  1:4305011/93‑1 (MQ=255)
                cGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCa                                                                         >  1:2485488/1‑93 (MQ=255)
                 gAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCAt                                                                        >  2:3059970/1‑93 (MQ=255)
                   aTCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATca                                                                      <  1:147718/93‑1 (MQ=255)
                   aTCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATca                                                                      >  1:2763138/1‑93 (MQ=255)
                   aTCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATca                                                                      >  1:2722347/1‑93 (MQ=255)
                   aTCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATca                                                                      >  2:840959/1‑93 (MQ=255)
                   aTCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATca                                                                      <  2:2883605/93‑1 (MQ=255)
                     cccGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATcacc                                                                    >  2:3401370/1‑92 (MQ=255)
                     cccGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATcaca                                                                    <  2:3657638/93‑1 (MQ=255)
                       cGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCCATCATCCCTGCCATCACAcc                                                                  <  2:2423651/93‑1 (MQ=255)
                       cGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGATGGCGATCATCCCTGCCATCACAcc                                                                  <  1:3787836/93‑1 (MQ=255)
                       cGCCAGGTTCATATACACCGGCCGCATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACAcc                                                                  <  2:2819475/93‑1 (MQ=255)
                        gCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACAccg                                                                 >  1:944607/1‑92 (MQ=255)
                          cAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCgg                                                               >  1:3706552/1‑93 (MQ=255)
                           aGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACACGGa                                                              >  2:1611456/1‑93 (MQ=255)
                             gTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACg                                                            >  1:4488299/1‑93 (MQ=255)
                               tCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGtt                                                          <  2:4416936/93‑1 (MQ=255)
                               tCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGtt                                                          <  2:775062/93‑1 (MQ=255)
                                cATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACATCCGGACGttt                                                         >  2:1177797/1‑92 (MQ=255)
                                 atatACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTga                                                        >  1:1605048/1‑92 (MQ=255)
                                 atatACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGt                                                        >  2:3525565/1‑93 (MQ=255)
                                  tataCACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGta                                                       <  2:91374/93‑1 (MQ=255)
                                   ataCCCCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGtaa                                                      <  2:1188208/93‑1 (MQ=255)
                                    ttaacCGGGGGAATATTCGGGAAGAAAGCCAACTGCCCGGGGAAAGGCATGAAGGTGGCCAGCAACCCTGCCATTCCGCCCGGACGTTGtaat                                                     <  2:1404868/90‑1 (MQ=255)
                                    taCACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGtaat                                                     >  1:3017641/1‑93 (MQ=255)
                                    taCACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGtaat                                                     <  1:3502540/93‑1 (MQ=255)
                                     acacCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGtaata                                                    <  1:2031385/93‑1 (MQ=255)
                                      cacCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGtaataa                                                   <  1:1072337/93‑1 (MQ=255)
                                      cacCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGtaataa                                                   <  2:2938548/93‑1 (MQ=255)
                                      cacCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGtaataa                                                   <  2:1730646/93‑1 (MQ=255)
                                        ccGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGc                                                 >  2:4377207/1‑93 (MQ=255)
                                         cGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCg                                                <  1:4704323/93‑1 (MQ=255)
                                         cGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCg                                                <  1:2335643/93‑1 (MQ=255)
                                         cGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCg                                                <  2:1638253/93‑1 (MQ=255)
                                            gCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGaa                                             <  1:955957/93‑1 (MQ=255)
                                               aaTATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGt                                          >  1:2195285/1‑93 (MQ=255)
                                                atatACGTGAAGACAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTa                                         <  1:1359252/93‑1 (MQ=255)
                                                atatACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTa                                         >  2:1701309/1‑93 (MQ=255)
                                                 tataCGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTAt                                        <  2:879355/93‑1 (MQ=255)
                                                 tataCGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTAt                                        >  2:3215160/1‑93 (MQ=255)
                                                  ataCGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTAtt                                       <  2:2455006/93‑1 (MQ=255)
                                                    aCGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTAtttt                                     <  1:1453658/93‑1 (MQ=255)
                                                    aCGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTAtttt                                     <  1:1948448/93‑1 (MQ=255)
                                                    aCGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTAtttt                                     <  2:2217823/93‑1 (MQ=255)
                                                    aCGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTAtttt                                     <  2:459063/93‑1 (MQ=255)
                                                      gTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGt                                   <  2:3775692/93‑1 (MQ=255)
                                                      gTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGt                                   <  2:599533/93‑1 (MQ=255)
                                                       tGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGtt                                  >  2:3158471/1‑93 (MQ=255)
                                                         aagaaAGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTc                                <  1:1419591/93‑1 (MQ=255)
                                                            aaaGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGa                             <  1:2271270/93‑1 (MQ=255)
                                                             aaGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGAt                            <  1:4613960/93‑1 (MQ=255)
                                                              aGCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATg                           <  2:4050989/93‑1 (MQ=255)
                                                               gCAAACTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGc                          >  1:1105356/1‑93 (MQ=255)
                                                                 aaaCTGCCCGGCGAACAACCTGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTAAGATGcaa                        <  2:4221645/93‑3 (MQ=255)
                                                                 aaaCTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGa                        <  1:4740033/93‑1 (MQ=255)
                                                                 aaaCTGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGa                        <  1:1051177/93‑1 (MQ=255)
                                                                     tGCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGt                    >  2:1105651/1‑93 (MQ=255)
                                                                      gCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGtt                   <  1:3141860/93‑1 (MQ=255)
                                                                      gCCCGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGtt                   <  1:3196854/93‑1 (MQ=255)
                                                                        ccGGCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCg                 >  2:2149747/1‑93 (MQ=255)
                                                                          ggCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGcc               >  1:445033/1‑93 (MQ=255)
                                                                           gCGAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCt              <  1:1213985/93‑1 (MQ=255)
                                                                            ggAAAGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTg             <  2:3954816/92‑1 (MQ=255)
                                                                             gAAAGACATGAAGGTGGCGATCATCCCTGCCATCGCACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTgg            <  1:4504357/93‑1 (MQ=255)
                                                                              aaaGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGc           >  2:546229/1‑93 (MQ=255)
                                                                              aaaGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGAAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGc           >  2:1252992/1‑93 (MQ=255)
                                                                               aaGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCa          >  2:1451012/1‑93 (MQ=255)
                                                                               aaGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCa          >  1:244330/1‑93 (MQ=255)
                                                                               aaGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGCATTTTGTTTCTGATGCGAGGGTTCGCCTGGCa          >  2:3604464/1‑93 (MQ=255)
                                                                                aGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAg         <  2:3551565/93‑1 (MQ=255)
                                                                                aGACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAg         >  2:73793/1‑93 (MQ=255)
                                                                                aGACACGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAg         <  1:1504355/93‑1 (MQ=255)
                                                                                 gACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGt        <  2:823447/93‑1 (MQ=255)
                                                                                 gACATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGt        >  2:267141/1‑93 (MQ=255)
                                                                                   cATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGTGa      <  2:4571252/93‑1 (MQ=255)
                                                                                   cATGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGTGa      >  1:1156025/1‑93 (MQ=255)
                                                                                    aTGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGTGaa     <  1:2858739/93‑1 (MQ=255)
                                                                                    aTGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGTGaa     <  1:1709692/93‑1 (MQ=255)
                                                                                     tGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGTGAAg    >  2:1920609/1‑93 (MQ=255)
                                                                                     tGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGTGAAg    >  1:3796101/1‑93 (MQ=255)
                                                                                     tGAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGTGAAg    >  2:1398654/1‑93 (MQ=255)
                                                                                      gAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGTGAAgg   <  1:426091/93‑1 (MQ=255)
                                                                                      gAAGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGGTCGCCTGGCAGTGAAgg   >  1:886168/1‑93 (MQ=255)
                                                                                       aaGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGTGAAGGc  >  1:568678/1‑93 (MQ=255)
                                                                                       aaGGTGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGTGAAGGc  >  1:1834431/1‑93 (MQ=255)
                                                                                          |                                                                                         
GTTAAGCCATCCACGCCGAATCCCGCCAGGTTCATATACACCGGGCGAATATACGTGAAGAAAGCAAACTGCCCGGCGAAAGACATGAAGATGGCGATCATCCCTGCCATCACACCCGGACGTTGTAATAAGCGGAAAGTATTTTGTTTCTGATGCGAGGGTTCGCCTGGCAGTGAAGGC  >  NC_000913/3838603‑3838782

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