Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 2,862,881 T→C L208P (CTA→CCA)  ygbM → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,862,8810TC99.7% 917.1 / NA 293L208P (CTA→CCA) ygbMhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (140/152);  minor base A (0/1);  total (140/153)
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.00e+00

Only 100 of 310 total aligned reads displayed.

TACTTTTCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACTACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGCTGGATCG  >  NC_000913/2862792‑2862970
                                                                                         |                                                                                         
tACTTTTCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGaccacc                                                                                        <  1:599466/93‑1 (MQ=255)
tACTTTTCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGaccacc                                                                                        >  2:7375362/1‑93 (MQ=255)
tACTTTTCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGaccacc                                                                                        <  2:5217337/93‑1 (MQ=255)
   ttttCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGGAATATGCGCACGTACAAATTGCCGGACCACCTGa                                                                                     >  1:5383263/1‑93 (MQ=255)
   ttttCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGa                                                                                     >  1:7439179/1‑93 (MQ=255)
   ttttCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGa                                                                                     >  2:6925887/1‑93 (MQ=255)
   ttttCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGa                                                                                     >  1:573423/1‑93 (MQ=255)
   ttttCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGa                                                                                     <  2:1753988/93‑1 (MQ=255)
    tttCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGGCGGACCACCTGAt                                                                                    >  2:1721587/1‑93 (MQ=255)
    tttCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGAt                                                                                    <  1:4557511/93‑1 (MQ=255)
    tttCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGAt                                                                                    <  2:4505624/93‑1 (MQ=255)
    tttCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGAt                                                                                    >  1:4392129/1‑93 (MQ=255)
    tttCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGAt                                                                                    <  2:7589867/93‑1 (MQ=255)
     ttCATGCGCAAAAAGTAGATGGCACCCTGACGCTTTTAATCCGCGACTACGCCGGGAAATGTGCGCACGTACAAATTGCCGGACCACCTGATc                                                                                   <  1:7038838/93‑1 (MQ=255)
     ttCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGTGAAATATGCGCACGTACAAATTGCCGGACCACCTGATc                                                                                   <  2:6089370/93‑1 (MQ=255)
     ttCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATc                                                                                   <  2:5253125/93‑1 (MQ=255)
         tacgcAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACa                                                                               <  1:6148819/91‑1 (MQ=255)
         tGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACa                                                                               <  1:4741166/93‑1 (MQ=255)
         tGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACa                                                                               <  1:3477719/93‑1 (MQ=255)
         tGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACa                                                                               <  1:246905/93‑1 (MQ=255)
            gcAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGa                                                                            <  1:5871049/93‑1 (MQ=255)
            gcAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGa                                                                            >  2:1287239/1‑93 (MQ=255)
            gcAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGa                                                                            <  2:6698619/93‑1 (MQ=255)
             cAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGaa                                                                           <  2:7886324/93‑1 (MQ=255)
              aaaaaGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAAc                                                                          >  1:1268459/1‑93 (MQ=255)
               aaaaGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAAcc                                                                         <  1:3304317/93‑1 (MQ=255)
               aaaaGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATAGACATGAAcc                                                                         >  1:6193183/1‑93 (MQ=255)
                  aGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGa                                                                      <  2:414161/93‑1 (MQ=255)
                  aGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGa                                                                      <  1:3676561/93‑1 (MQ=255)
                      gATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGAt                                                                  >  2:1031107/1‑93 (MQ=255)
                       aTGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATg                                                                 >  2:2164750/1‑93 (MQ=255)
                       aTGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATg                                                                 <  2:380568/93‑1 (MQ=255)
                        tGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATgg                                                                >  1:6056348/1‑93 (MQ=255)
                         ggCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGga                                                               <  1:879237/93‑1 (MQ=255)
                          gCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGgag                                                              <  2:6470878/93‑1 (MQ=255)
                           cAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGgaga                                                             <  1:1748298/93‑1 (MQ=255)
                            aaCCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGaa                                                            >  2:5571193/1‑93 (MQ=255)
                            aaCCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGaa                                                            >  1:2240534/1‑93 (MQ=255)
                            aaCCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGaa                                                            <  2:7237530/93‑1 (MQ=255)
                             aCCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGaaa                                                           <  2:4886886/93‑1 (MQ=255)
                              ccTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAAt                                                          >  1:1093128/1‑93 (MQ=255)
                               cTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATc                                                         >  1:3769981/1‑93 (MQ=255)
                                 gACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCaa                                                       >  1:6711216/1‑93 (MQ=255)
                                 gACGCATTTAATCCGCGACTACGCCGGGAAAAATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGaaaacca                                                       >  2:6174611/1‑89 (MQ=255)
                                   cGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACt                                                     <  1:6312670/93‑1 (MQ=255)
                                   cGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACt                                                     >  2:4368101/1‑93 (MQ=255)
                                   cGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACt                                                     >  1:1650878/1‑93 (MQ=255)
                                   cGCATTTAATCCGCGACTACGCCGGGAAATATGCACACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACt                                                     <  1:3747845/93‑1 (MQ=255)
                                     cATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTAt                                                   >  2:1385391/1‑93 (MQ=255)
                                     cATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTAt                                                   >  1:504363/1‑93 (MQ=255)
                                     cATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTAt                                                   <  1:1734086/93‑1 (MQ=255)
                                      aTTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATc                                                  <  2:5827619/93‑1 (MQ=255)
                                       tttAATCCGCGACTACGCCGGGATATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATcc                                                 <  2:4293232/93‑1 (MQ=255)
                                       tttAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATcc                                                 >  2:2361522/1‑93 (MQ=255)
                                         tAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGt                                               <  1:3562466/93‑1 (MQ=255)
                                         tAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGt                                               <  1:2806508/93‑1 (MQ=255)
                                              cgcgACTACGCCGGGAAATATGCGCTCGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTg                                          <  1:197069/93‑1 (MQ=255)
                                                 gACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTc                                       <  2:5308531/93‑1 (MQ=255)
                                                 gACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTc                                       >  1:5135724/1‑93 (MQ=255)
                                                  aCTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACGGGACGATGGAGCAATCAACTATCCGTGGCTGTTcc                                      >  1:7848745/1‑93 (MQ=255)
                                                  aCTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTcc                                      <  1:7664781/93‑1 (MQ=255)
                                                   cTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTcgg                                     >  2:3843979/1‑91 (MQ=255)
                                                    tACGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGc                                    >  2:6562059/1‑93 (MQ=255)
                                                     aCGCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGcc                                   >  1:2670054/1‑93 (MQ=255)
                                                       gCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTg                                 <  2:1849679/93‑1 (MQ=255)
                                                       gCCGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTg                                 <  2:259329/93‑1 (MQ=255)
                                                        ccGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGt                                >  2:2756610/1‑93 (MQ=255)
                                                         cGGGAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGtt                               >  1:1314517/1‑93 (MQ=255)
                                                          gggAAATATGCGCACGTGAAAATTGCCGGACCACCTGATGGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTc                              <  2:2019289/93‑1 (MQ=255)
                                                          gggAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTc                              >  2:5375882/1‑93 (MQ=255)
                                                            gAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGa                            <  2:1550202/93‑1 (MQ=255)
                                                            gAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGa                            >  1:3966853/1‑93 (MQ=255)
                                                            gAAATATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGa                            <  2:7721322/93‑1 (MQ=255)
                                                            gAAATATGCGCACGTACAAATTGCCGGACAACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGa                            <  2:4065668/93‑1 (MQ=255)
                                                            gAAATATGCGCACGGACAAATTGCCGGACCACCTGATCGACATGACCCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGa                            <  1:6361436/93‑1 (MQ=255)
                                                             aaaTATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGAt                           <  1:7588294/93‑1 (MQ=255)
                                                              aaTATGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATg                          >  2:6608226/1‑93 (MQ=255)
                                                               atatGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGa                         >  2:338736/1‑93 (MQ=255)
                                                               atatGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGa                         <  2:1437189/93‑1 (MQ=255)
                                                                tatGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAg                        >  1:2942858/1‑93 (MQ=255)
                                                                tatGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAg                        <  2:3917953/93‑1 (MQ=255)
                                                                tatGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAg                        >  2:3179191/1‑93 (MQ=255)
                                                                tatGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAg                        <  1:716209/93‑1 (MQ=255)
                                                                tatGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATACGTGGCTGTTCCGCCTGTTCGCTGAg                        >  1:4846694/1‑93 (MQ=255)
                                                                 atGCGCACGTACAAATTGCCGGACCTCCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAgg                       <  1:4939229/93‑1 (MQ=255)
                                                                  tGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGt                      <  1:7531258/93‑1 (MQ=255)
                                                                  tGCGCACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGt                      <  1:7370080/93‑1 (MQ=255)
                                                                     gcACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTggg                   <  2:5730536/93‑1 (MQ=255)
                                                                     gcACGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTggg                   >  2:881602/1‑93 (MQ=255)
                                                                       aCGTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGat                 >  1:4323138/1‑93 (MQ=255)
                                                                         gTACAAATTGCCGGACCCCCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGatat               <  2:4760318/93‑1 (MQ=255)
                                                                         gTACAAATTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGatat               <  2:7531468/93‑1 (MQ=255)
                                                                            cAAAATGCCGGCCCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAg            <  2:5946416/93‑1 (MQ=255)
                                                                             aaaTTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAgg           >  2:6264663/1‑93 (MQ=255)
                                                                             aaaTTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAgg           >  1:3148194/1‑93 (MQ=255)
                                                                              aaTTGCCGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAggg          >  1:5590554/1‑93 (MQ=255)
                                                                                   ccGGACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGCTGGa     <  2:7082399/93‑1 (MQ=255)
                                                                                     ggACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGCTGGATc   <  1:3875051/93‑1 (MQ=255)
                                                                                      gACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGCTGGATCg  <  2:2942303/93‑1 (MQ=255)
                                                                                      gACCACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGCTGGATCg  <  2:1908613/93‑1 (MQ=255)
                                                                                         |                                                                                         
TACTTTTCATGCGCAAAAAGTAGATGGCAACCTGACGCATTTAATCCGCGACTACGCCGGGAAATATGCGCACGTACAAATTGCCGGACTACCTGATCGACATGAACCGGACGATGGAGAAATCAACTATCCGTGGCTGTTCCGCCTGTTCGATGAGGTGGGATATCAGGGCTGGATCG  >  NC_000913/2862792‑2862970

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