Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,374,263 T→C Q50R (CAA→CGA)  sspB ← stringent starvation protein B

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,374,2630TC99.4% 1166.2 / NA 358Q50R (CAA→CGA) sspBstringent starvation protein B
Reads supporting (aligned to +/- strand):  ref base T (0/0);  major base C (176/180);  minor base A (2/0);  total (178/180)
Fisher's exact test for biased strand distribution p-value = 2.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 374 total aligned reads displayed.

GGAATGCCACCAAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTTGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCAGCAAC  >  NC_000913/3374173‑3374352
                                                                                          |                                                                                         
 gAATGCCACCAAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGcc                                                                                        <  2:5826285/93‑1 (MQ=255)
  aaTGCCACCAAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGccc                                                                                       <  2:5744529/93‑1 (MQ=255)
  aaTGCCACCAAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGccc                                                                                       <  1:8616305/93‑1 (MQ=255)
  aaTGCCACCAAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGccc                                                                                       <  1:5815412/93‑1 (MQ=255)
   aTGCCACCAAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCg                                                                                      >  1:6707788/1‑93 (MQ=255)
     gCCACCAAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTc                                                                                    >  1:6206162/1‑93 (MQ=255)
      ccaccaAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCa                                                                                   >  2:9026677/1‑93 (MQ=255)
      ccaccaAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCa                                                                                   <  1:3801653/93‑1 (MQ=255)
       caccaAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCAc                                                                                  >  1:9648659/1‑93 (MQ=255)
       caccaAAGCGCGCGTTAAAGCGCACCTCATCATGCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCAc                                                                                  <  1:7478609/93‑1 (MQ=255)
         ccaAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCAcgc                                                                                >  1:1458647/1‑93 (MQ=255)
         ccaAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCAcgc                                                                                <  2:2476671/93‑1 (MQ=255)
          caAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGACCGTCAcgcg                                                                               >  1:2816229/1‑93 (MQ=255)
           aaaGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACgcgc                                                                              <  1:4542687/93‑1 (MQ=255)
            aaGCGCGCGTTAAAGCGCTCCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCa                                                                             <  1:4333093/93‑1 (MQ=255)
            aaGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCa                                                                             <  1:7984438/93‑1 (MQ=255)
              gcgcgcGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCata                                                                           <  1:1360634/93‑1 (MQ=255)
               cgcgcgTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCatat                                                                          >  1:7296885/1‑93 (MQ=255)
                gcgcgTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATAtt                                                                         <  2:1020164/93‑1 (MQ=255)
                gcgcgTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATAtt                                                                         <  1:5206587/93‑1 (MQ=255)
                gcgcgTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATAtt                                                                         <  2:4469380/93‑1 (MQ=255)
                gcgcgTTAAAGCGCACCTCATCATTCGCCAGTTACAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATAtt                                                                         >  2:8824437/1‑93 (MQ=255)
                  gcgGTAAACCTCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTcc                                                                       <  1:3284136/93‑1 (MQ=255)
                   cgTTAAAGCACACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCa                                                                      <  1:7935665/93‑1 (MQ=255)
                     ttAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATa                                                                    >  2:7514913/1‑93 (MQ=255)
                     ttAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATa                                                                    >  1:8958329/1‑93 (MQ=255)
                     ttAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATa                                                                    >  1:3927161/1‑93 (MQ=255)
                     ttAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATa                                                                    >  1:6538802/1‑93 (MQ=255)
                        aaGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGa                                                                 >  2:9544451/1‑93 (MQ=255)
                         aGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGaa                                                                >  1:9751144/1‑93 (MQ=255)
                          gcgcACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAAc                                                               >  1:4995558/1‑93 (MQ=255)
                            gcACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCt                                                             <  2:5552564/93‑1 (MQ=255)
                              cacccatcatTCGCCATTTCCAGATTGCCCACAGCACCCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGc                                                           <  1:918972/89‑1 (MQ=255)
                              aCCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGc                                                           <  1:8091930/93‑1 (MQ=255)
                              aCCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAACGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGc                                                           >  2:3649712/1‑93 (MQ=255)
                               ccTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCa                                                          >  1:8488611/1‑93 (MQ=255)
                                cTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCAc                                                         >  2:1731195/1‑93 (MQ=255)
                                cTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCAc                                                         <  1:1893193/93‑1 (MQ=255)
                                 tcatcaTTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACg                                                        >  1:106771/1‑93 (MQ=255)
                                    tcatTCGCCAGTTCTAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCa                                                     <  2:6900506/93‑1 (MQ=255)
                                    tcatTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCACATATTCCATAGGAACCTGCACGCCa                                                     >  1:6048734/1‑93 (MQ=255)
                                     catTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAg                                                    >  1:4746478/1‑93 (MQ=255)
                                     catTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAg                                                    >  1:6934225/1‑93 (MQ=255)
                                       ttCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAggg                                                  >  2:3903550/1‑93 (MQ=255)
                                        tCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGa                                                 <  2:4053576/93‑1 (MQ=255)
                                          gCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGc                                               >  2:2063350/1‑93 (MQ=255)
                                           ccAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCg                                              >  1:7855224/1‑93 (MQ=255)
                                           ccAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCg                                              >  2:4782106/1‑93 (MQ=255)
                                           ccAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCg                                              >  2:2656891/1‑93 (MQ=255)
                                             aGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTc                                            >  1:7240604/1‑93 (MQ=255)
                                              gTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTca                                           >  2:4194750/1‑93 (MQ=255)
                                              gTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTca                                           >  2:2846951/1‑93 (MQ=255)
                                               ttCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTcac                                          <  1:2133380/93‑1 (MQ=255)
                                               ttCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTcac                                          >  1:369177/1‑93 (MQ=255)
                                                 ccAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTGACAt                                        >  2:3963784/1‑93 (MQ=255)
                                                 ccAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACAt                                        <  2:1202090/93‑1 (MQ=255)
                                                   aGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATcc                                      >  1:2681270/1‑93 (MQ=255)
                                                    gATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATcca                                     <  1:8351974/93‑1 (MQ=255)
                                                       tGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATccacca                                  <  1:707240/93‑1 (MQ=255)
                                                       tGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATccacca                                  <  2:8498335/93‑1 (MQ=255)
                                                         ccGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATccaccacc                                >  1:562240/1‑93 (MQ=255)
                                                          cGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATccaccacca                               >  1:5286713/1‑93 (MQ=255)
                                                           gACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAg                              >  2:6518276/1‑93 (MQ=255)
                                                            aCAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAgg                             <  2:6097207/93‑1 (MQ=255)
                                                            aCAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAgg                             <  2:5278409/93‑1 (MQ=255)
                                                             cAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGt                            <  1:9836216/93‑1 (MQ=255)
                                                              aGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTg                           <  2:6288648/93‑1 (MQ=255)
                                                              aGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTg                           >  2:4228120/1‑93 (MQ=255)
                                                              aGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTg                           >  2:7995468/1‑93 (MQ=255)
                                                               gCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTgc                          <  1:8794515/93‑1 (MQ=255)
                                                               gCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTgc                          <  1:7751363/93‑1 (MQ=255)
                                                               gCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTgc                          <  1:1789481/93‑1 (MQ=255)
                                                                cACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTgcg                         <  1:4877789/93‑1 (MQ=255)
                                                                  cgcgGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTgcggc                       >  1:5036222/1‑93 (MQ=255)
                                                                    cggcgCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGt                     <  2:7644430/93‑1 (MQ=255)
                                                                    cggcgCAAGGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGt                     <  2:7603108/93‑1 (MQ=255)
                                                                       cgcAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAg                  <  1:758777/93‑1 (MQ=255)
                                                                       cgcAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAg                  >  2:7775759/1‑93 (MQ=255)
                                                                       cgcAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAg                  <  1:4630773/93‑1 (MQ=255)
                                                                        gcAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGc                 <  2:9628956/93‑1 (MQ=255)
                                                                          aaTGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTg               >  2:1476586/1‑93 (MQ=255)
                                                                          aaTGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTg               >  2:5338013/1‑93 (MQ=255)
                                                                           aTGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTgg              >  2:5000929/1‑93 (MQ=255)
                                                                           aTGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTgg              >  2:3061304/1‑93 (MQ=255)
                                                                           aTGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTgg              >  2:3588559/1‑93 (MQ=255)
                                                                            tGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGt             >  1:1598028/1‑93 (MQ=255)
                                                                            tGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGt             <  2:7396232/93‑1 (MQ=255)
                                                                             gTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGtt            <  2:7324006/93‑1 (MQ=255)
                                                                             gTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGtt            >  2:8299184/1‑93 (MQ=255)
                                                                               tGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCTGCGTGAGCTGGTTAt          >  2:6787129/1‑93 (MQ=255)
                                                                                 aGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATcc        <  1:5443864/93‑1 (MQ=255)
                                                                                 aGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATcc        >  2:336527/1‑93 (MQ=255)
                                                                                   tACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCAg      >  2:8110762/1‑93 (MQ=255)
                                                                                    aCGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCAGc     <  2:7333466/93‑1 (MQ=255)
                                                                                      gATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCAGCaa   <  2:8641622/93‑1 (MQ=255)
                                                                                       aTTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATGCAGCAAc  >  2:6604529/1‑93 (MQ=255)
                                                                                       aTTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCAGCAAc  >  1:3041706/1‑93 (MQ=255)
                                                                                       aTTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCAGCAAc  >  1:9561083/1‑93 (MQ=255)
                                                                                       aTTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCAGCAAc  <  1:2338560/93‑1 (MQ=255)
                                                                                       aTTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGGTATCCCGCAAc  >  1:5517080/1‑93 (MQ=255)
                                                                                          |                                                                                         
GGAATGCCACCAAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTTGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCAGCAAC  >  NC_000913/3374173‑3374352

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