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,0350TC99.2% 1343.9 / ‑2.7 391R89R (CGT→CGCyhhYhypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (1/2);  new base C (180/208);  total (181/210)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.87e-01

Only 100 of 403 total aligned reads displayed.

TGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAAC  >  NC_000913/3578945‑3579124
                                                                                          |                                                                                         
 gtagaCGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTgg                                                                                        >  2:3241537/3‑93 (MQ=255)
  gagaCGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGga                                                                                       >  1:3977160/1‑93 (MQ=255)
   agaCGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGgaa                                                                                      <  2:1593203/93‑1 (MQ=255)
    gaCGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGgcag                                                                                     >  2:3651062/1‑93 (MQ=255)
      cgtcgtGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATTTGTGTCGACTCTCGCTGgaagaa                                                                                   >  2:4728058/1‑93 (MQ=255)
      cgtcgtGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGgaagaa                                                                                   <  2:3086503/93‑1 (MQ=255)
       ttcgtGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAAc                                                                                  <  2:1338366/92‑1 (MQ=255)
        tcgtGGGCCATCTCACCATTGACGTGCAACAGCGCCCCCGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAAcc                                                                                 <  2:622050/93‑1 (MQ=255)
        tcgtGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAAcc                                                                                 >  1:203371/1‑93 (MQ=255)
        gcgtGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAAcc                                                                                 <  1:8033788/92‑1 (MQ=255)
         cgtGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACcg                                                                                <  2:5983208/93‑1 (MQ=255)
          gtGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACcgc                                                                               >  1:985910/1‑93 (MQ=255)
          gtGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTAGAAGAACcgc                                                                               >  1:1051033/1‑93 (MQ=255)
             ggCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCgcggc                                                                            <  2:823758/93‑1 (MQ=255)
             ggCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCgcggc                                                                            >  1:5702251/1‑93 (MQ=255)
               gcATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGt                                                                          <  1:1305197/92‑1 (MQ=255)
               ccATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGt                                                                          <  1:7793827/93‑1 (MQ=255)
               ccATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGt                                                                          <  2:4044000/93‑1 (MQ=255)
                cATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTc                                                                         <  1:3778819/93‑1 (MQ=255)
                 aTCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCg                                                                        >  1:3746209/1‑93 (MQ=255)
                  tctcACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGTGGCGTCGc                                                                       >  1:8205333/1‑93 (MQ=255)
                     cACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAg                                                                    >  1:4634502/1‑93 (MQ=255)
                     cACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAg                                                                    >  1:487914/1‑93 (MQ=255)
                     cACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAg                                                                    >  2:2651229/1‑93 (MQ=255)
                      aCCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAgc                                                                   >  1:269165/1‑93 (MQ=255)
                          ttGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGccc                                                               >  2:6147034/1‑93 (MQ=255)
                          ttGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGccc                                                               <  1:686287/93‑1 (MQ=255)
                          ttGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGccc                                                               <  1:1387352/93‑1 (MQ=255)
                          ttGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGccc                                                               <  1:8693217/93‑1 (MQ=255)
                           tGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCt                                                              <  2:7376840/93‑1 (MQ=255)
                            gACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTg                                                             <  2:5472148/93‑1 (MQ=255)
                            gACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTg                                                             >  2:6320993/1‑93 (MQ=255)
                            gACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTg                                                             >  1:5912831/1‑93 (MQ=255)
                             aCGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGACCCGCGGCGTCGCCAGCGCCCTGa                                                            <  2:412767/93‑1 (MQ=255)
                             aCGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGa                                                            <  1:2484644/93‑1 (MQ=255)
                              cGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGAt                                                           <  2:4354482/93‑1 (MQ=255)
                                tGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGc                                                         <  2:847042/93‑1 (MQ=255)
                                tGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGc                                                         <  2:4760162/93‑1 (MQ=255)
                                    aCAGCGCCGACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCgaga                                                     <  2:5852442/93‑1 (MQ=255)
                                    aCAGCGCCCACTCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGCCGCCAGCGCCCTGATGCgaga                                                     <  2:8654033/93‑1 (MQ=255)
                                    aCAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCgaga                                                     >  1:2842812/1‑93 (MQ=255)
                                    aCAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCgaga                                                     >  1:3588741/1‑93 (MQ=255)
                                      aGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCgagaga                                                   >  2:4530303/1‑93 (MQ=255)
                                       gcgcCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAgat                                                  <  1:1463138/93‑1 (MQ=255)
                                         gcCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAgatga                                                <  1:5321458/93‑1 (MQ=255)
                                          cccACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAgatgat                                               <  2:5638750/93‑1 (MQ=255)
                                             aCGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGa                                            >  2:679998/1‑93 (MQ=255)
                                             aCGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGa                                            >  1:8358040/1‑93 (MQ=255)
                                             aCGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGa                                            <  1:1932876/93‑1 (MQ=255)
                                               gccgcAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGaaa                                          <  2:1141571/93‑1 (MQ=255)
                                                ccgcAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAAt                                         <  1:7075405/93‑1 (MQ=255)
                                                ccgcAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAAt                                         <  2:3357002/93‑1 (MQ=255)
                                                ccgcAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAAt                                         <  2:2121729/93‑1 (MQ=255)
                                                 cgcAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAAtg                                        <  2:6463112/93‑1 (MQ=255)
                                                  gcAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAAtgt                                       >  1:6125523/1‑93 (MQ=255)
                                                  gcAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAAtgt                                       <  2:763462/93‑1 (MQ=255)
                                                   cAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAAtgtg                                      >  1:6226166/1‑93 (MQ=255)
                                                    aGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGc                                     <  1:4479587/93‑1 (MQ=255)
                                                    aGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGc                                     <  2:7747048/93‑1 (MQ=255)
                                                     gTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCg                                    >  1:8175642/1‑93 (MQ=255)
                                                      tCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGa                                   <  1:2288389/93‑1 (MQ=255)
                                                        aTGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACa                                 >  1:7148461/1‑93 (MQ=255)
                                                         tGTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACaa                                >  1:7006326/1‑93 (MQ=255)
                                                          gTTGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAAc                               <  1:3300097/93‑1 (MQ=255)
                                                           ttGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACt                              <  2:2912372/93‑1 (MQ=255)
                                                           ttGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACt                              <  1:5985994/93‑1 (MQ=255)
                                                           ttGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACt                              <  1:6723387/93‑1 (MQ=255)
                                                            tgtcGATTTTGGTGTCTTTGTCGGCTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTg                             <  1:2339123/90‑1 (MQ=255)
                                                            tGCCGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTg                             <  2:5721213/93‑1 (MQ=255)
                                                              cggATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGt                           <  2:5301064/91‑1 (MQ=255)
                                                              ccGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGt                           <  1:1793463/93‑1 (MQ=255)
                                                              ccGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGt                           <  1:3204419/93‑1 (MQ=255)
                                                              ccGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGt                           >  1:8301932/1‑93 (MQ=255)
                                                              ccGATTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGt                           <  2:1683992/93‑1 (MQ=255)
                                                                 aTTTTGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGCCACCTGGTTGc                        >  1:7387959/1‑93 (MQ=255)
                                                                  ttttGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCg                       >  1:2210515/1‑93 (MQ=255)
                                                                  ttttGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCg                       >  1:8516028/1‑93 (MQ=255)
                                                                   tttGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGGCAACTGGTTGCgg                      <  1:3450974/93‑1 (MQ=255)
                                                                   tttGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCgg                      <  2:7289901/93‑1 (MQ=255)
                                                                    ttGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCggg                     <  2:1249345/93‑1 (MQ=255)
                                                                    ttGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCggg                     <  2:8060333/93‑1 (MQ=255)
                                                                     tGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGt                    <  1:8661255/93‑1 (MQ=255)
                                                                     tGGTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGt                    <  1:651014/93‑1 (MQ=255)
                                                                       gTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAg                  <  1:5158271/93‑1 (MQ=255)
                                                                       gTATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAg                  <  1:5556547/93‑1 (MQ=255)
                                                                        tATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGa                 <  2:7209393/93‑1 (MQ=255)
                                                                        tATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGa                 >  1:4806965/1‑93 (MQ=255)
                                                                        tATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGa                 <  2:525799/93‑1 (MQ=255)
                                                                        tATCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGa                 >  2:2350999/1‑93 (MQ=255)
                                                                          tCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATc               <  2:704102/93‑1 (MQ=255)
                                                                          tCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATc               <  2:459527/93‑1 (MQ=255)
                                                                          tCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATc               >  2:1959663/1‑93 (MQ=255)
                                                                          tCTGTGTCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATc               <  2:8358040/93‑1 (MQ=255)
                                                                          tCTGTGTCGACTCTCGCTGGAAGAACAGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGCTGCGGGGAGATc               >  1:7694654/1‑93 (MQ=255)
                                                                              tgtCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCAt           >  1:5124131/1‑93 (MQ=255)
                                                                              tgtCGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCAt           <  2:3861347/93‑1 (MQ=255)
                                                                                 cGACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGa        <  1:7591517/93‑1 (MQ=255)
                                                                                  gACTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGaa       <  2:7544886/93‑1 (MQ=255)
                                                                                   aCTCTCGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAAc      >  2:2242629/1‑93 (MQ=255)
                                                                                      ctcGCTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTaa   >  1:4275660/1‑93 (MQ=255)
                                                                                          |                                                                                         
TGGAGACGTCGTGGGCCATCTCACCATTGACGTGCAACAGCGCCCACGCCGCAGTCATGTTGCCGATTTTGGTATCTGTGTCGACTCTCGTTGGAAGAACCGCGGCGTCGCCAGCGCCCTGATGCGAGAGATGATTGAAATGTGCGACAACTGGTTGCGGGTAGATCGCATTGAACTAAC  >  NC_000913/3578945‑3579124

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