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.6% 820.6 / ‑5.4 252Q50R (CAA→CGA) sspBstringent starvation protein B
Reads supporting (aligned to +/- strand):  ref base T (0/1);  new base C (130/121);  total (130/122)
Fisher's exact test for biased strand distribution p-value = 4.84e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.86e-01

Only 100 of 260 total aligned reads displayed.

GGAATGCCACCAAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTTGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCAGCAAC  >  NC_000913/3374173‑3374352
                                                                                          |                                                                                         
ggAATGCCACCAAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGc                                                                                         >  2:3331776/1‑93 (MQ=255)
ggAATGCCACCAAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGc                                                                                         >  1:424553/1‑93 (MQ=255)
ggAATGCCACCAAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGc                                                                                         >  1:5416612/1‑93 (MQ=255)
 gAATGCCACCAAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGcc                                                                                        <  1:351054/93‑1 (MQ=255)
 gAATGCCACCAAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGcc                                                                                        <  2:5484032/93‑1 (MQ=255)
   aTGCCACCAAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCg                                                                                      >  1:1679320/1‑93 (MQ=255)
   aTGCCACCAAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCg                                                                                      >  1:5363327/1‑93 (MQ=255)
   aTGCCACCAAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCg                                                                                      <  1:4997746/93‑1 (MQ=255)
     gCCACCAAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTc                                                                                    >  2:2657767/1‑93 (MQ=255)
      ccaccaAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCa                                                                                   >  1:1642383/1‑93 (MQ=255)
      ccaccaAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCa                                                                                   >  1:6332144/1‑93 (MQ=255)
       caccaAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCAc                                                                                  >  1:2794034/1‑93 (MQ=255)
       caccaAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCAc                                                                                  >  2:157876/1‑93 (MQ=255)
        accaAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCAcg                                                                                 >  2:2038867/1‑93 (MQ=255)
        accaAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCAcg                                                                                 >  1:2306200/1‑93 (MQ=255)
          caAAGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCAcgcg                                                                               >  1:703909/1‑93 (MQ=255)
           aaaGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACgcgc                                                                              >  1:6297381/1‑93 (MQ=255)
            aaGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCa                                                                             >  2:4134852/1‑93 (MQ=255)
            aaGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCa                                                                             >  2:2934961/1‑93 (MQ=255)
             aGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCat                                                                            >  2:4373860/1‑93 (MQ=255)
             aGCGCGCGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCat                                                                            <  1:2515208/93‑1 (MQ=255)
              gcgcgcGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCata                                                                           >  2:2475413/1‑93 (MQ=255)
              gcgcgcGTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCata                                                                           >  2:4064938/1‑93 (MQ=255)
                  gcgTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTcc                                                                       >  1:182936/1‑93 (MQ=255)
                   cgTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCa                                                                      >  2:689990/1‑93 (MQ=255)
                    gTTAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCAt                                                                     <  1:1954871/93‑1 (MQ=255)
                     ttAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATa                                                                    <  1:1255436/93‑1 (MQ=255)
                     ttAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATa                                                                    >  1:6210025/1‑93 (MQ=255)
                     ttAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATa                                                                    <  2:5271144/93‑1 (MQ=255)
                      tAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAg                                                                   <  2:1727248/93‑1 (MQ=255)
                      tAAAGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAg                                                                   <  2:4785363/93‑1 (MQ=255)
                       aaaGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAgg                                                                  >  1:5183263/1‑93 (MQ=255)
                        aaGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGa                                                                 >  2:3255937/1‑93 (MQ=255)
                        aaGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGa                                                                 <  1:5712610/93‑1 (MQ=255)
                         aGCGCACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGaa                                                                >  2:5097436/1‑93 (MQ=255)
                           cgcACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAAcc                                                              <  2:5970448/93‑1 (MQ=255)
                           cgcACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAAcc                                                              <  1:3629594/93‑1 (MQ=255)
                            gcACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCt                                                             >  1:2782710/1‑93 (MQ=255)
                            gcACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCt                                                             <  1:1505840/93‑1 (MQ=255)
                            gcACCTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCt                                                             <  1:931520/93‑1 (MQ=255)
                               ccTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCa                                                          <  2:6155566/93‑1 (MQ=255)
                               ccTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCa                                                          >  2:6023134/1‑93 (MQ=255)
                               ccTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCa                                                          <  2:4826185/93‑1 (MQ=255)
                                cTCATCATTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCAc                                                         <  2:4293156/93‑1 (MQ=255)
                                  catcatTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGc                                                       >  1:3314511/1‑93 (MQ=255)
                                   atcatTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGcc                                                      <  2:6414101/93‑1 (MQ=255)
                                    tcatTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCa                                                     >  2:1386386/1‑93 (MQ=255)
                                     catTCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAg                                                    <  1:4337837/93‑1 (MQ=255)
                                       ttCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAggg                                                  <  1:5574368/93‑1 (MQ=255)
                                       ttCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAggg                                                  <  2:6160923/93‑1 (MQ=255)
                                       ttCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAggg                                                  <  2:4896191/93‑1 (MQ=255)
                                        tCGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGa                                                 >  1:122773/1‑93 (MQ=255)
                                         cGCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAg                                                <  1:2204173/93‑1 (MQ=255)
                                          gCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGc                                               >  2:5949665/1‑93 (MQ=255)
                                          gCCAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCACGGAGc                                               >  2:2990853/1‑93 (MQ=255)
                                           ccAGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCg                                              >  2:6457316/1‑93 (MQ=255)
                                             aGTTCCAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTc                                            >  1:2060766/1‑93 (MQ=255)
                                                 ccAGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCTGCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACAt                                        <  2:785362/93‑1 (MQ=255)
                                                   aGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATcc                                      <  1:1116177/93‑1 (MQ=255)
                                                   aGATTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCACGGAGCGTCACATcc                                      >  1:4473524/1‑93 (MQ=255)
                                                     aTTGCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATccac                                    <  1:6487783/93‑1 (MQ=255)
                                                        gCCGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATccaccac                                 >  2:1924128/1‑93 (MQ=255)
                                                         ccGACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATccaccacc                                >  2:4508317/1‑93 (MQ=255)
                                                           gACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAg                              >  2:5446776/1‑93 (MQ=255)
                                                           gACAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAg                              >  1:3702698/1‑93 (MQ=255)
                                                            aCAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAgg                             >  1:6104719/1‑93 (MQ=255)
                                                             cAGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGt                            <  2:4178554/93‑1 (MQ=255)
                                                              aGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTg                           >  1:4153891/1‑93 (MQ=255)
                                                              aGCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTg                           <  1:307568/93‑1 (MQ=255)
                                                               gCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTgc                          <  2:3185802/93‑1 (MQ=255)
                                                               gCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTgc                          >  2:3763789/1‑93 (MQ=255)
                                                               gCACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTgc                          >  2:6611794/1‑93 (MQ=255)
                                                                cACGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTgcg                         >  2:2626630/1‑93 (MQ=255)
                                                                 aCGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTgcgg                        >  2:6091164/1‑93 (MQ=255)
                                                                 aCGCGGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTgcgg                        <  2:3082040/93‑1 (MQ=255)
                                                                  cgcgGCGCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTgcggc                       <  2:2183820/93‑1 (MQ=255)
                                                                    cggcgCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGt                     <  1:5067217/93‑1 (MQ=255)
                                                                     ggcgCAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTg                    <  2:5264865/93‑1 (MQ=255)
                                                                      gcgcAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGa                   <  2:5744489/93‑1 (MQ=255)
                                                                        gcAATGTTGAGTACGATTTGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGc                 <  1:5851246/93‑1 (MQ=255)
                                                                         cAATGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCt                <  1:1567988/93‑1 (MQ=255)
                                                                           aTGTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTgg              >  1:2658630/1‑93 (MQ=255)
                                                                             gTTGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGtt            <  2:1218063/93‑1 (MQ=255)
                                                                              ttGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTa           >  2:942117/1‑93 (MQ=255)
                                                                              ttGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTa           <  2:5903644/93‑1 (MQ=255)
                                                                               tGAGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTAt          <  1:5149623/93‑1 (MQ=255)
                                                                                 aGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATcc        >  1:1821678/1‑93 (MQ=255)
                                                                                 aGTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATcc        <  2:2219469/93‑1 (MQ=255)
                                                                                  gTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCa       >  2:1536259/1‑93 (MQ=255)
                                                                                  gTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCa       >  1:663044/1‑93 (MQ=255)
                                                                                  gTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCa       <  1:5438987/93‑1 (MQ=255)
                                                                                  gTACGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCa       <  1:4770669/93‑1 (MQ=255)
                                                                                     cGATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCAGCa    <  1:3890621/93‑1 (MQ=255)
                                                                                      gATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCAGCaa   >  1:4068746/1‑93 (MQ=255)
                                                                                      gATTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCAGCaa   >  2:817048/1‑93 (MQ=255)
                                                                                      gATTCGCCCGTCACACGCATATTCCATAGGAACCTGCAGGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGGTGGTTATCCAGCaa   >  1:6131199/1‑93 (MQ=255)
                                                                                       aTTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCAGCAAc  <  1:6688885/93‑1 (MQ=255)
                                                                                       aTTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCAGCAAc  <  2:4660588/93‑1 (MQ=255)
                                                                                       aTTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCAGCAAc  >  2:6277192/1‑93 (MQ=255)
                                                                                       aTTCGCCCGTCACGCGCATATTCCATAGGAACCTGCACGCCAGGGAGCGTCACATCCACCACCAGGTGCGGCGTGAGCTGGTTATCCAGCAAc  <  1:3809463/93‑1 (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 < 28 ≤ ATCG/ATCG < 33 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: