Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,579,035 T→C R89R (CGT→CGC yhhY → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,579,0350TC89.3% 785.9 / 75.7 299R89R (CGT→CGCyhhYhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (13/19);  new base C (117/150);  total (130/169)
Fisher's exact test for biased strand distribution p-value = 8.51e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.22e-01

Only 100 of 300 total aligned reads displayed.

TGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGT  >  NC_000913/3578945‑3579127
                                                                                          |                                                                                            
tGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTg                                                                                            >  1:1239876/1‑93 (MQ=255)
tGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTg                                                                                            >  1:3597123/1‑93 (MQ=255)
tGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTg                                                                                            <  2:2488057/93‑1 (MQ=255)
   agaCGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGgaa                                                                                         <  2:330772/93‑1 (MQ=255)
   agaCGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGgaa                                                                                         <  2:2140711/93‑1 (MQ=255)
    gaCGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGgaag                                                                                        >  2:661392/1‑93 (MQ=255)
     aCGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGgaaga                                                                                       >  2:1168125/1‑93 (MQ=255)
      cgtcgtGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGgaagaa                                                                                      <  1:1207777/93‑1 (MQ=255)
      cgtcgtGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGgaagaa                                                                                      <  2:1025067/93‑1 (MQ=255)
       gtcgtGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAAc                                                                                     <  2:6108111/93‑1 (MQ=255)
       gtcgtGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAAc                                                                                     >  2:6682128/1‑93 (MQ=255)
        tcgtGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAAcc                                                                                    <  2:2596076/93‑1 (MQ=255)
         cgtGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACcg                                                                                   <  1:4964262/93‑1 (MQ=255)
         cgtGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACcg                                                                                   <  1:2019379/93‑1 (MQ=255)
          gtGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACcgc                                                                                  <  1:1820619/93‑1 (MQ=255)
           tGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACcgcg                                                                                 >  1:2635447/1‑93 (MQ=255)
            gggCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCgcgg                                                                                >  1:5552262/1‑93 (MQ=255)
            gggCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCgcgg                                                                                <  2:5053107/93‑1 (MQ=255)
              gCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCgcggcg                                                                              <  1:5329055/93‑1 (MQ=255)
               ccATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGt                                                                             <  1:3036880/93‑1 (MQ=255)
                cATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTc                                                                            <  1:2258110/93‑1 (MQ=255)
                  tctcACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGc                                                                          <  2:3258115/93‑1 (MQ=255)
                  tctcACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGc                                                                          <  2:7526297/93‑1 (MQ=255)
                   ctcACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGcc                                                                         <  2:5503964/93‑1 (MQ=255)
                   ctcACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGcc                                                                         <  1:810988/93‑1 (MQ=255)
                      aCCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCTATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAgc                                                                      <  2:4240354/93‑1 (MQ=255)
                      aCCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAgc                                                                      >  1:1974303/1‑93 (MQ=255)
                      aCCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAgc                                                                      >  2:4451381/1‑93 (MQ=255)
                        cATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAgcgc                                                                    >  2:2224135/1‑93 (MQ=255)
                         aTTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGcc                                                                   >  1:6798657/1‑93 (MQ=255)
                           tGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCt                                                                 >  2:2910471/1‑93 (MQ=255)
                           tGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCt                                                                 <  1:5951305/93‑1 (MQ=255)
                            gACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTg                                                                >  2:3005610/1‑93 (MQ=255)
                            gACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTg                                                                >  2:2686300/1‑93 (MQ=255)
                             aCGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGa                                                               <  1:2053982/93‑1 (MQ=255)
                             aCGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGa                                                               <  1:6442993/93‑1 (MQ=255)
                              cGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGAt                                                              <  1:46972/93‑1 (MQ=255)
                              cGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGAt                                                              <  1:1038961/93‑1 (MQ=255)
                              cGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGAt                                                              <  2:5778972/93‑1 (MQ=255)
                              cGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGAt                                                              <  1:6145619/93‑1 (MQ=255)
                               gTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATg                                                             >  2:7168450/1‑93 (MQ=255)
                                tGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGc                                                            >  2:5218234/1‑93 (MQ=255)
                                tGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTACCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGc                                                            <  1:7371928/93‑1 (MQ=255)
                                  cAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCga                                                          <  1:2764364/93‑1 (MQ=255)
                                   aaCAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCgag                                                         <  1:6926395/93‑1 (MQ=255)
                                    aCAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCgaga                                                        <  2:3503787/93‑1 (MQ=255)
                                    aCAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCgaga                                                        >  2:38006/1‑93 (MQ=255)
                                     cAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCgagag                                                       <  1:482550/93‑1 (MQ=255)
                                       gcgcCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAgat                                                     <  2:6408194/93‑1 (MQ=255)
                                         gcCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAgatga                                                   <  2:368319/93‑1 (MQ=255)
                                         gcCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAgatga                                                   <  2:1454603/93‑1 (MQ=255)
                                         gcCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAgatga                                                   <  2:6007497/93‑1 (MQ=255)
                                         gcCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAgatga                                                   <  1:3714804/93‑1 (MQ=255)
                                          cccACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAgatgat                                                  <  2:7292807/93‑1 (MQ=255)
                                          cccACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAgatgat                                                  >  2:5941211/1‑93 (MQ=255)
                                           ccACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGAtt                                                 <  2:128981/93‑1 (MQ=255)
                                               gccgcAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGaaa                                             <  2:3197298/93‑1 (MQ=255)
                                                 cgcAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAAtg                                           <  2:1771042/93‑1 (MQ=255)
                                                  gcAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAAtgt                                          >  2:4090124/1‑93 (MQ=255)
                                                  gcAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAAtgt                                          <  2:6179516/93‑1 (MQ=255)
                                                  gcAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAAtgt                                          >  1:915453/1‑93 (MQ=255)
                                                  gcAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAAtgt                                          >  1:5794505/1‑93 (MQ=255)
                                                  gcAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAAtgt                                          >  2:2104732/1‑93 (MQ=255)
                                                   cAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAAtgtg                                         <  1:2134049/93‑1 (MQ=255)
                                                    aGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGc                                        <  1:3240271/93‑1 (MQ=255)
                                                      tCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGa                                      <  1:4020709/93‑1 (MQ=255)
                                                       cATTTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGAc                                     >  1:6969063/1‑93 (MQ=255)
                                                       cATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGAc                                     >  1:6239473/1‑93 (MQ=255)
                                                       cATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGAc                                     <  1:4418682/93‑1 (MQ=255)
                                                       cATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGAc                                     <  1:1698959/93‑1 (MQ=255)
                                                         tGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACaa                                   >  1:1431821/1‑93 (MQ=255)
                                                          gTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAAc                                  <  1:6092801/93‑1 (MQ=255)
                                                           ttGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACt                                 <  2:7387098/93‑1 (MQ=255)
                                                           ttGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACt                                 <  1:4203477/93‑1 (MQ=255)
                                                            tGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTg                                <  2:1340100/93‑1 (MQ=255)
                                                            tGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTg                                <  2:200715/93‑1 (MQ=255)
                                                            tGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTg                                >  1:5187608/1‑93 (MQ=255)
                                                             gCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTgg                               <  1:3852527/93‑1 (MQ=255)
                                                             gCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTgg                               <  2:241374/93‑1 (MQ=255)
                                                              ccGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGt                              <  1:660398/93‑1 (MQ=255)
                                                               cGATTTTGGTATCTGTGTCGGCTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGtt                             >  2:6879512/1‑93 (MQ=255)
                                                               cGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGtt                             >  2:1609182/1‑93 (MQ=255)
                                                               cGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGtt                             <  2:5155941/93‑1 (MQ=255)
                                                               cGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGtt                             <  2:2893941/93‑1 (MQ=255)
                                                                  ttttGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCg                          >  2:411254/1‑93 (MQ=255)
                                                                  ttttGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCg                          <  1:5618909/93‑1 (MQ=255)
                                                                   tttGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCgg                         <  2:6434388/93‑1 (MQ=255)
                                                                     tGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGt                       <  1:3434810/93‑1 (MQ=255)
                                                                     tGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGt                       <  1:2872557/93‑1 (MQ=255)
                                                                      ggTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTa                      <  2:1412612/93‑1 (MQ=255)
                                                                      ggTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTa                      <  2:56583/93‑1 (MQ=255)
                                                                        tATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGa                    <  1:7435527/93‑1 (MQ=255)
                                                                          tCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATc                  >  2:649123/1‑93 (MQ=255)
                                                                             gtgtCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCa               >  2:3915090/1‑93 (MQ=255)
                                                                                 cGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGa           >  2:166128/1‑93 (MQ=255)
                                                                                   aCTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAAc         <  1:4178888/93‑1 (MQ=255)
                                                                                   aCTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAAc         <  1:5734791/93‑1 (MQ=255)
                                                                                      ctcGCTGGAAGAACCGCGGCGTCGTCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTaa      <  2:4911064/93‑1 (MQ=255)
                                                                                      ctcGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTaa      <  2:4382439/93‑1 (MQ=255)
                                                                                       tcGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAAc     >  2:4615211/1‑93 (MQ=255)
                                                                                          |                                                                                            
TGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAACCGT  >  NC_000913/3578945‑3579127

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