Predicted mutation
evidence seq id position mutation annotation gene description
RA NC_000913 3,554,113 Δ1 bp coding (372/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,1130C.99.6% 1120.0 / ‑3.3 252coding (372/396 nt)yhgKJRNA phosphage cyclase
Reads supporting (aligned to +/- strand):  ref base C (0/0);  major base . (119/132);  minor base A (0/1);  total (119/133)
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.96e-01

Only 100 of 281 total aligned reads displayed.

TGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGAGGCGGCGCCGACGGGTTATCGGTCGCCGCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGGC  >  NC_000913/3554026‑3554199
                                                                                         |                                                                                        
tGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgccgc                                                                                    <  1:3075529/93‑1 (MQ=255)
tGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGtgccgc                                                                                    >  2:3652426/1‑93 (MQ=255)
 ggTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgccgcc                                                                                   >  1:5820635/1‑92 (MQ=255)
   tttGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGCAGGGTTATCGGtgccgccgc                                                                                 >  1:2192139/1‑90 (MQ=255)
   tttGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgccgccgc                                                                                 >  1:6415666/1‑90 (MQ=255)
   tttGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgccgccgc                                                                                 <  1:5705079/93‑4 (MQ=255)
    ttGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgccgccgct                                                                                >  1:1552409/1‑89 (MQ=255)
    ttGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgccgccgct                                                                                >  2:6845863/1‑89 (MQ=255)
      gCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCa                                                                              <  1:3432080/93‑1 (MQ=255)
        tgatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACt                                                                            >  2:1355231/1‑93 (MQ=255)
         gatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACtt                                                                           >  2:6657810/1‑93 (MQ=255)
         gatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACtt                                                                           >  1:6977695/1‑93 (MQ=255)
          atgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTc                                                                          <  2:1098504/93‑1 (MQ=255)
             aaTTCCTATGTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAac                                                                       <  1:1914806/93‑1 (MQ=255)
              aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAaca                                                                      <  1:206063/93‑1 (MQ=255)
               ttCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAacac                                                                     >  2:634146/1‑93 (MQ=255)
                tCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACAcg                                                                    >  1:2789640/1‑93 (MQ=255)
                  cTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACAcgcg                                                                  >  1:4552477/1‑93 (MQ=255)
                   tATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGa                                                                 >  1:2801326/1‑93 (MQ=255)
                    aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCACAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGaa                                                                <  2:6998001/93‑1 (MQ=255)
                    aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGaa                                                                <  2:5045316/93‑1 (MQ=255)
                     ttttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAg                                                               >  1:2501137/1‑93 (MQ=255)
                     ttttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAg                                                               <  1:4161315/93‑1 (MQ=255)
                       ttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGt                                                             <  1:1702058/93‑1 (MQ=255)
                         cGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTcc                                                           <  2:127929/93‑1 (MQ=255)
                         cGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTcc                                                           <  1:2173119/93‑1 (MQ=255)
                         cGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACATCAACACGCGAAGGTcc                                                           >  2:2177777/1‑93 (MQ=255)
                          gCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCa                                                          <  2:6712922/93‑1 (MQ=255)
                           ccAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCAt                                                         <  1:3660486/93‑1 (MQ=255)
                           ccAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCAt                                                         <  1:4455036/93‑1 (MQ=255)
                           ccAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCAt                                                         >  1:5988828/1‑93 (MQ=255)
                           ccAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCAt                                                         <  1:1784586/93‑1 (MQ=255)
                            cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATc                                                        <  2:2270083/93‑1 (MQ=255)
                            cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATc                                                        >  2:7620422/1‑93 (MQ=255)
                            cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATc                                                        >  1:5142002/1‑93 (MQ=255)
                            cagcagCGCCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATc                                                        <  2:3028899/93‑1 (MQ=255)
                             agcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCg                                                       <  1:4442882/93‑1 (MQ=255)
                              gcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCgg                                                      <  2:4040150/93‑1 (MQ=255)
                                  cGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaaa                                                  <  1:5214176/93‑1 (MQ=255)
                                  cGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaaa                                                  >  2:4185207/1‑93 (MQ=255)
                                  cGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaaa                                                  >  1:7930754/1‑93 (MQ=255)
                                    gCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAcc                                                <  2:522844/93‑1 (MQ=255)
                                     cTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACca                                               <  2:392357/93‑1 (MQ=255)
                                      tCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcac                                              >  2:1642190/1‑93 (MQ=255)
                                      tCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcac                                              <  2:7495904/93‑1 (MQ=255)
                                      tCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcac                                              >  1:1095504/1‑93 (MQ=255)
                                      tCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCAGCTCACTTCAACACACGAAGGTCCATCGGCAcaccac                                              >  2:2227861/1‑93 (MQ=255)
                                        cAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAg                                            >  2:3760894/1‑93 (MQ=255)
                                        cAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAg                                            <  2:4246956/93‑1 (MQ=255)
                                        cAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAg                                            <  2:778142/93‑1 (MQ=255)
                                        cAGCACCCGGCGGATAAAACCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACAACAg                                            >  1:2537186/1‑93 (MQ=255)
                                         aGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgc                                           <  1:679531/93‑1 (MQ=255)
                                          gCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgcg                                          >  2:3347426/1‑93 (MQ=255)
                                           cACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgcgc                                         <  2:5519168/93‑1 (MQ=255)
                                            aCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGcgcg                                        >  1:6159226/1‑93 (MQ=255)
                                            aCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGcgcg                                        <  2:3476468/93‑1 (MQ=255)
                                            aCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGcgcg                                        >  2:3269450/1‑93 (MQ=255)
                                             cccGGCGGATAAAATCCGCGGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCgg                                       <  1:7498065/93‑1 (MQ=255)
                                               cggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgc                                     >  1:7330571/1‑93 (MQ=255)
                                               cggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgc                                     <  2:4483008/93‑1 (MQ=255)
                                                ggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgca                                    >  1:894754/1‑93 (MQ=255)
                                                ggcggATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGAGCCGCCGCTCACTTCAACACGCGAAGCTCCATCGGCAAACCACAGCGCGAgca                                    >  2:1985500/1‑93 (MQ=255)
                                                   ggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagcc                                 >  2:5935142/1‑92 (MQ=255)
                                                   caacaaaaTCCGCAGGCGGCGCAGTCGCGTTACCGGTGCCGCCTCTCACTTAAACACGCGAGGGTCCATCGGCAATCCACAGCGCGGgcagca                                 <  2:4913526/89‑1 (MQ=255)
                                                    gATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc                                >  1:7816421/1‑93 (MQ=255)
                                                    gATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc                                >  2:6155472/1‑93 (MQ=255)
                                                        aaaTCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGtt                            <  2:5780946/93‑1 (MQ=255)
                                                        aaaTCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATAGGCAGACCACAGCGCGGGCAGCACCGtt                            >  2:6913949/1‑93 (MQ=255)
                                                        aaaTCCGCAGGCGGCGCCGACGGGTTAACGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGtt                            <  2:4395044/93‑1 (MQ=255)
                                                             ccgcaGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgca                        <  2:1251205/89‑1 (MQ=255)
                                                             ccgcaGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgca                        <  2:274503/89‑1 (MQ=255)
                                                             ccgcaGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgca                        <  2:2867518/89‑1 (MQ=255)
                                                              cgcaGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgcag                       >  2:816429/4‑93 (MQ=255)
                                                                  ggcggcGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCAcc                   >  1:3954473/1‑93 (MQ=255)
                                                                  ggcggcGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGACCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCAcc                   >  2:2345278/1‑93 (MQ=255)
                                                                  gcccgcgcCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCAcc                   <  1:3214744/89‑1 (MQ=255)
                                                                   gcggcgCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCa                  <  1:4902714/93‑1 (MQ=255)
                                                                   gcggcgCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCa                  <  1:5572046/93‑1 (MQ=255)
                                                                   gcggcgCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCa                  <  2:1856904/93‑1 (MQ=255)
                                                                    gggcgCGGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAACCCACAGCGCGGGCAGCACCGTTTGCAGCACCAg                 <  2:4768354/92‑1 (MQ=255)
                                                                    cggcgCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAg                 <  1:7239525/93‑1 (MQ=255)
                                                                    cggcgCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAg                 <  1:2340973/93‑1 (MQ=255)
                                                                    cggcgCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAg                 >  2:6216293/1‑93 (MQ=255)
                                                                      gcgcCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCg               >  2:3153078/1‑93 (MQ=255)
                                                                      gcgcCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCg               <  2:5650285/93‑1 (MQ=255)
                                                                        gcCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTa             <  2:5373216/93‑1 (MQ=255)
                                                                         ccGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTAc            >  1:1012390/1‑93 (MQ=255)
                                                                          cGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACa           >  2:6732934/1‑93 (MQ=255)
                                                                            aCGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAAc         >  1:3844993/1‑93 (MQ=255)
                                                                             cGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACt        <  1:2631619/93‑1 (MQ=255)
                                                                               ggTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTc      <  1:181294/93‑1 (MQ=255)
                                                                                gTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTcc     <  2:7257889/93‑1 (MQ=255)
                                                                                gTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTcc     <  1:1304569/93‑1 (MQ=255)
                                                                                gTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTcc     <  1:4276555/93‑1 (MQ=255)
                                                                                gTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTcc     <  2:629414/93‑1 (MQ=255)
                                                                                gTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTcc     >  2:516989/1‑93 (MQ=255)
                                                                                gTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGAACCGTTTGCAGCACCAGCGCACAACTTcc     >  1:6378674/1‑93 (MQ=255)
                                                                                  tATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCgg   >  1:6952354/1‑93 (MQ=255)
                                                                                   aTCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGgc  <  2:2713971/93‑1 (MQ=255)
                                                                                   aTCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGgc  >  1:7847066/1‑93 (MQ=255)
                                                                                         |                                                                                        
TGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGAGGCGGCGCCGACGGGTTATCGGTCGCCGCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGGC  >  NC_000913/3554026‑3554199

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