Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,554,118 +C coding (367/396 nt) yhgKJ ← RNA phosphage cyclase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,554,1171.C99.2% 898.6 / ‑2.1 268coding (368/396 nt)yhgKJRNA phosphage cyclase
Reads supporting (aligned to +/- strand):  ref base . (2/0);  new base C (124/142);  total (126/142)
Fisher's exact test for biased strand distribution p-value = 2.20e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.99e-01

Only 100 of 303 total aligned reads displayed.

TGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGAGGCGGCGCCGACGGGTTATCGGTCGCCGCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGGC  >  NC_000913/3554026‑3554199
                                                                                             |                                                                                  
 ggTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgccgcc                                                                                  <  1:4963481/93‑2 (MQ=255)
   tttGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgccgccgc                                                                                >  1:5661454/1‑90 (MQ=255)
   tttGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgccgccgc                                                                                >  2:7395790/1‑90 (MQ=255)
    ttGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgccgccgct                                                                               >  1:2153995/1‑89 (MQ=255)
      gCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCTGCGGATAAAATCCGCAGGCGGCGCCGACGGGTGATCTGTGCCGCCGCTCa                                                                              >  2:5387494/1‑93 (MQ=255)
      gCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCa                                                                              <  2:4558154/93‑1 (MQ=255)
       cTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCAc                                                                             <  1:8500144/93‑1 (MQ=255)
         gatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACtt                                                                           <  1:3888065/93‑1 (MQ=255)
           tgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCa                                                                         >  2:8478721/1‑93 (MQ=255)
            gaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCaa                                                                        >  2:7761845/1‑93 (MQ=255)
             aaTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAac                                                                       <  2:3275996/93‑1 (MQ=255)
             aaTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAac                                                                       >  2:4685155/1‑93 (MQ=255)
              aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAaca                                                                      >  2:6891479/1‑93 (MQ=255)
                 ccTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAacacac                                                                   >  1:3387619/1‑91 (MQ=255)
                 ccTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACAcgc                                                                   >  1:8120036/1‑93 (MQ=255)
                 ccTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACAcgc                                                                   <  2:1719285/93‑1 (MQ=255)
                  cTATTTTCGCCAGCAGCGGCTCCAGCCCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACAcgcg                                                                  <  2:6525025/93‑1 (MQ=255)
                   tATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGa                                                                 >  2:623714/1‑93 (MQ=255)
                   tATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGa                                                                 <  2:7101898/93‑1 (MQ=255)
                     ttttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGCAg                                                               >  1:6097849/1‑93 (MQ=255)
                      tttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAgg                                                              <  1:6440841/93‑1 (MQ=255)
                        tCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTc                                                            <  1:4669455/93‑1 (MQ=255)
                          gCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCa                                                          >  2:6814897/1‑93 (MQ=255)
                           ccaccaGCGGCTCCAGCGCCCGGCGGATAAAGTCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCAt                                                         <  2:4961724/93‑1 (MQ=255)
                           ccAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCAt                                                         >  1:4864892/1‑93 (MQ=255)
                           ccAGCAGCGGCCCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCAt                                                         <  1:6310807/93‑1 (MQ=255)
                            cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATc                                                        <  1:6032488/93‑1 (MQ=255)
                             agcagcGGCTCCAGCACCCGGCGGAAAAAATCCGCAGGCGGCGCCGACGGGTTATCGGAGCCGCCGCTAAATTCCACACGCGAAAGCGCATCg                                                       >  2:6133536/1‑93 (MQ=255)
                              gcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCgg                                                      >  1:8351707/1‑93 (MQ=255)
                              gcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCgg                                                      >  2:7847827/1‑93 (MQ=255)
                              gcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCgg                                                      >  1:6345286/1‑93 (MQ=255)
                              gcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCgg                                                      >  1:1922544/1‑93 (MQ=255)
                              gcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGCGCCGCCGCTCCCTTCCACACGCCACGGTCCATCgg                                                      >  1:6209952/1‑93 (MQ=255)
                               cagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGc                                                     >  2:418381/1‑93 (MQ=255)
                                 gcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaa                                                   >  2:1973746/1‑93 (MQ=255)
                                   ggCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAc                                                 >  1:1275698/1‑93 (MQ=255)
                                     cTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACca                                               <  1:8413317/93‑1 (MQ=255)
                                      tCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGGCCATCGGCAAACcac                                              >  1:7713485/1‑93 (MQ=255)
                                       ccAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcaca                                             <  1:1123817/93‑1 (MQ=255)
                                       ccAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcaca                                             <  1:2611992/93‑1 (MQ=255)
                                        cAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAg                                            <  1:7416766/93‑1 (MQ=255)
                                        cAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAg                                            >  2:7261877/1‑93 (MQ=255)
                                         aGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgc                                           >  1:1574834/1‑93 (MQ=255)
                                         aGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACAACAgc                                           >  2:2401287/1‑93 (MQ=255)
                                          gCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgcg                                          <  2:2368381/93‑1 (MQ=255)
                                          gCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgcg                                          >  1:4011752/1‑93 (MQ=255)
                                           cACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgcgc                                         <  2:2919875/93‑1 (MQ=255)
                                           cACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgcgc                                         <  2:7492288/93‑1 (MQ=255)
                                            aCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGcgcg                                        <  2:7542557/93‑1 (MQ=255)
                                            aCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGcgcg                                        >  1:4510733/1‑93 (MQ=255)
                                            aCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGcgcg                                        >  1:7997224/1‑93 (MQ=255)
                                             cccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCCCAGCGCgg                                       >  2:1309734/1‑93 (MQ=255)
                                             cccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCgg                                       <  2:6045801/93‑1 (MQ=255)
                                             cccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCgg                                       >  1:6797457/1‑93 (MQ=255)
                                             cccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCgg                                       >  1:442373/1‑93 (MQ=255)
                                              ccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGGCCATCGGCAAACCACAGCGCggg                                      >  2:6310840/1‑93 (MQ=255)
                                               cggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgc                                     <  1:1048613/93‑1 (MQ=255)
                                               cggcggATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgc                                     >  1:5234472/1‑93 (MQ=255)
                                                 gcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcag                                   >  1:3187067/1‑93 (MQ=255)
                                                   ggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagca                                 >  2:5779600/1‑93 (MQ=255)
                                                    gATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc                                <  1:379499/93‑1 (MQ=255)
                                                    gATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc                                <  2:3900920/93‑1 (MQ=255)
                                                     aTAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAcc                               >  1:912450/1‑93 (MQ=255)
                                                      tAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCg                              >  1:6967180/1‑93 (MQ=255)
                                                        aaaTCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGtt                            >  1:506733/1‑93 (MQ=255)
                                                        aaaTCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTAACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGtt                            <  2:5290677/93‑1 (MQ=255)
                                                         aaTCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGttt                           <  2:1093667/93‑1 (MQ=255)
                                                         aaTCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGAACCGttt                           >  2:6393700/1‑93 (MQ=255)
                                                            ccgcaGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgca                        <  2:4827631/89‑1 (MQ=255)
                                                            ccgcaGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgca                        <  2:1195072/89‑1 (MQ=255)
                                                             cgcaGGCGGCGCCGACGGGTTATCGGTGCCGGCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTGTgcag                       >  1:1624743/4‑93 (MQ=255)
                                                             cgcaGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgcag                       >  1:2497063/4‑93 (MQ=255)
                                                             cgcaGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgcag                       <  2:2724186/90‑1 (MQ=255)
                                                             cgcaGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgcag                       >  1:7872144/4‑93 (MQ=255)
                                                               caGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgcagca                     >  1:4124674/2‑93 (MQ=255)
                                                               caGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgcagca                     <  2:2314062/92‑1 (MQ=255)
                                                                 ggcggcGCGGACGTGTTATCGGTGCCGCCGCCTACCTGAACACGCGAAGGGCCAGCGGCAAACCACAGCGCGGGCAGCACCGGTTGCAGCAcc                   <  1:6567914/93‑1 (MQ=255)
                                                                 ggcggcGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCAcc                   >  1:8648277/1‑93 (MQ=255)
                                                                    ggcgCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGc                >  1:4334748/1‑93 (MQ=255)
                                                                    ggcgCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGc                <  2:184467/93‑1 (MQ=255)
                                                                     gcgcCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCg               <  2:1688256/93‑1 (MQ=255)
                                                                      ggcCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGt              <  2:7939579/92‑1 (MQ=255)
                                                                       gcCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTa             <  2:933996/93‑1 (MQ=255)
                                                                       gcCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTa             <  2:6604041/93‑1 (MQ=255)
                                                                       gcCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTa             <  2:3611572/93‑1 (MQ=255)
                                                                        ccGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTAc            <  2:3408281/93‑1 (MQ=255)
                                                                         cGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACa           <  2:5666861/93‑1 (MQ=255)
                                                                          gACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACaa          <  2:1487418/93‑1 (MQ=255)
                                                                          gACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACaa          >  1:5929478/1‑93 (MQ=255)
                                                                          gACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACaa          <  2:8195728/93‑1 (MQ=255)
                                                                           aCGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAAc         <  1:5521625/93‑1 (MQ=255)
                                                                           aCGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAAc         <  2:564415/93‑1 (MQ=255)
                                                                           aCGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAAc         <  2:8313414/93‑1 (MQ=255)
                                                                             gggTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACtt       >  2:4456690/1‑93 (MQ=255)
                                                                               gTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTcc     <  1:3557715/93‑1 (MQ=255)
                                                                               gTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTcc     <  1:2390983/93‑1 (MQ=255)
                                                                               gTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTcc     <  1:1464426/93‑1 (MQ=255)
                                                                               gTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTcc     <  2:8537814/93‑1 (MQ=255)
                                                                                  aTCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGgc  >  1:6632743/1‑93 (MQ=255)
                                                                                  aTCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGgc  >  2:3643687/1‑93 (MQ=255)
                                                                                             |                                                                                  
TGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGAGGCGGCGCCGACGGGTTATCGGTCGCCGCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGGC  >  NC_000913/3554026‑3554199

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