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.C98.5% 1252.1 / 7.3 335coding (368/396 nt)yhgKJRNA phosphage cyclase
Reads supporting (aligned to +/- strand):  ref base . (3/2);  new base C (155/175);  total (158/177)
Fisher's exact test for biased strand distribution p-value = 6.70e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.26e-01

Only 100 of 380 total aligned reads displayed.

GCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGAGGCGGCGCCGACGGGTTATCGGTCGCCGCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGGCGCT  >  NC_000913/3554023‑3554202
                                                                                                |                                                                                     
 cGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgccgcc                                                                                     <  1:4693848/96‑2 (MQ=255)
 agTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgccgcc                                                                                     <  1:2085248/95‑2 (MQ=255)
    ggTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgccgccgct                                                                                  >  2:8366213/1‑92 (MQ=255)
    ggTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgccgccgct                                                                                  <  1:4899077/96‑5 (MQ=255)
      tttGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCa                                                                                 <  1:12150735/96‑1 (MQ=255)
         gCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGTGCCGCCGCTCACtt                                                                              >  1:7169867/1‑96 (MQ=255)
          cggatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGACGGCGGCGACGGTGTATCGGTGCCGCCGCTCACTTc                                                                             >  1:10259785/3‑96 (MQ=255)
            gatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCaa                                                                           >  2:5690844/1‑96 (MQ=255)
            gatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCaa                                                                           >  1:5826639/1‑96 (MQ=255)
             atgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAac                                                                          <  2:8060489/96‑1 (MQ=255)
               gaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAacac                                                                        <  2:2608209/96‑1 (MQ=255)
                 aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGTCGACGGGTTATCGGTGCCGCCGCTCACTTCAACAcgc                                                                      >  2:12232100/1‑96 (MQ=255)
                 aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACAcgc                                                                      <  1:6337667/96‑1 (MQ=255)
                 aTTCCTATTTTCGCCAGCAGCGCCTCGAGCACCCGGCGCATAAAATCCGCAGGCGGCGCCGCCGGGTTATCGGTGCCGCCGCTCACTTCAACAcgc                                                                      <  1:12450072/96‑1 (MQ=255)
                  ttCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACAcgcg                                                                     <  1:9436926/96‑1 (MQ=255)
                      tATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAgg                                                                 >  2:10495964/1‑96 (MQ=255)
                      tATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAgg                                                                 <  2:959315/96‑1 (MQ=255)
                      tATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGTGCCGCCGCTCACTTCAACACGCGAAgg                                                                 >  2:11518000/1‑96 (MQ=255)
                       aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGt                                                                <  1:12908208/96‑1 (MQ=255)
                       aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGt                                                                <  2:5899062/96‑1 (MQ=255)
                       aTTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGt                                                                >  1:8106563/1‑96 (MQ=255)
                           tCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCAt                                                            <  2:10822357/96‑1 (MQ=255)
                            cGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGCCGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATc                                                           <  1:2517514/96‑1 (MQ=255)
                            cGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATc                                                           <  1:6615410/96‑1 (MQ=255)
                             gCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCg                                                          <  1:9199598/96‑1 (MQ=255)
                             gCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCg                                                          >  2:6977185/1‑96 (MQ=255)
                              ccAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCgg                                                         >  1:13181904/1‑96 (MQ=255)
                              ccAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCgg                                                         <  2:1632156/96‑1 (MQ=255)
                               cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGc                                                        <  2:11112095/96‑1 (MQ=255)
                               cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGc                                                        >  2:7937651/1‑96 (MQ=255)
                                 gcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaa                                                      >  1:4354928/1‑96 (MQ=255)
                                  cagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaaa                                                     >  1:7722399/1‑96 (MQ=255)
                                   agcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAc                                                    >  2:10268060/1‑96 (MQ=255)
                                   agcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAc                                                    >  1:1075131/1‑96 (MQ=255)
                                   agcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGTGCCGCCGCTCACTTCAACACGCGGAGGTCCATCGGGAAAc                                                    >  1:2666795/1‑96 (MQ=255)
                                    gcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAcc                                                   <  1:11600048/96‑1 (MQ=255)
                                     cGGCTCCAGCCCCCGGCGGATAAACTCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACca                                                  <  1:389975/96‑1 (MQ=255)
                                     cGGATCCAGCACCCGGCGGATAAAAACCGCAGGCGGCGCAGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACca                                                  <  1:1229798/96‑1 (MQ=255)
                                      ggCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAcccc                                                 <  1:11076942/96‑3 (MQ=255)
                                      ggCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcac                                                 >  2:5535685/1‑96 (MQ=255)
                                      ggCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcac                                                 >  2:8617049/1‑96 (MQ=255)
                                      ggCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGTGCCGCCGCTCACTTCAACACGCGAAAGTCCATCGGCAAACcac                                                 >  1:1850661/1‑96 (MQ=255)
                                       gCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcaca                                                >  2:11484595/1‑96 (MQ=255)
                                       gCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcaca                                                >  2:5073448/1‑96 (MQ=255)
                                         tCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgc                                              <  2:7806313/96‑1 (MQ=255)
                                         tCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgc                                              <  2:1781248/96‑1 (MQ=255)
                                         tCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGCCCATCGGCAAACCACAgc                                              >  1:4645251/1‑96 (MQ=255)
                                         tCCAGCACCCGGCGGATAAAACCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgc                                              <  2:1214023/96‑1 (MQ=255)
                                         tCCAACACCCGGCGGATACAATCCGCAGGCGGCGCCGACAGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgc                                              >  1:10137656/1‑96 (MQ=255)
                                            aGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACAACAGcgcg                                           <  2:2169892/96‑1 (MQ=255)
                                             gcccccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCgg                                          <  2:4761088/93‑1 (MQ=255)
                                             gCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCgg                                          <  1:10364858/96‑1 (MQ=255)
                                             gCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCgg                                          <  2:12723398/96‑1 (MQ=255)
                                             gCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCCCgg                                          <  2:10109740/96‑1 (MQ=255)
                                               gcccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgc                                        <  1:852849/95‑1 (MQ=255)
                                               aCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgc                                        <  1:6809235/96‑1 (MQ=255)
                                                cccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAACACAGCGCGGgca                                       >  2:626345/1‑96 (MQ=255)
                                                 ccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcag                                      <  1:1553999/96‑1 (MQ=255)
                                                  cggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTTAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagc                                     <  2:14028606/96‑1 (MQ=255)
                                                   ggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagca                                    >  2:9785578/1‑96 (MQ=255)
                                                   ggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagca                                    <  2:9211370/96‑1 (MQ=255)
                                                   ggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagca                                    <  2:13859306/96‑1 (MQ=255)
                                                   ggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagca                                    <  1:9667558/96‑1 (MQ=255)
                                                    gcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc                                   <  2:11884741/96‑1 (MQ=255)
                                                    gcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc                                   >  1:13126183/1‑96 (MQ=255)
                                                    gcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc                                   >  2:3665680/1‑96 (MQ=255)
                                                     cggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAcc                                  >  2:1335019/1‑96 (MQ=255)
                                                     cggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAcc                                  <  2:2342112/96‑1 (MQ=255)
                                                         tAAAGTCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGttt                              <  1:12950587/96‑1 (MQ=255)
                                                         tAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGttt                              <  2:7408961/96‑1 (MQ=255)
                                                          aaaaTCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCTGCAAACCACAGCGCGGGCAGCACCGCTTg                             >  1:4763804/1‑96 (MQ=255)
                                                           aaaTCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgc                            >  1:8832035/1‑96 (MQ=255)
                                                            aaTCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgca                           <  2:13585788/96‑1 (MQ=255)
                                                               ccgcaGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTgcagca                        >  1:11907634/5‑96 (MQ=255)
                                                                  caGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCa                     <  2:13234407/95‑1 (MQ=255)
                                                                   aGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAg                    >  1:3025059/1‑96 (MQ=255)
                                                                   aGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGCCCATCGCCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAg                    >  1:12502463/1‑96 (MQ=255)
                                                                    ggcggcGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCTAACCACAGCGCGGGCAGCACCGTTTGCagcagcagc                   >  1:12696742/1‑96 (MQ=255)
                                                                    ggcggcGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGc                   <  2:3217858/96‑1 (MQ=255)
                                                                      cggcgCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCgc                 <  1:11374964/96‑2 (MQ=255)
                                                                      cggcgCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGt                 >  1:7134002/1‑96 (MQ=255)
                                                                      cggcgCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGt                 >  1:3691975/1‑96 (MQ=255)
                                                                      cggcgCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGt                 <  2:7169867/96‑1 (MQ=255)
                                                                      cggcgCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGt                 >  2:284137/1‑96 (MQ=255)
                                                                       ggcgCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTa                >  1:5320476/1‑96 (MQ=255)
                                                                        gcgcCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTAc               >  2:5280762/1‑96 (MQ=255)
                                                                          gcCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACaa             >  1:11530888/1‑96 (MQ=255)
                                                                          gcCGACGGGGTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGTCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACaa             >  1:1583758/1‑96 (MQ=255)
                                                                           ccGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAAc            >  1:7427614/1‑96 (MQ=255)
                                                                             gACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACtt          >  2:4641899/1‑96 (MQ=255)
                                                                              aCGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTc         <  2:5477323/96‑1 (MQ=255)
                                                                               ctggTTACCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTcc        <  1:13965458/94‑1 (MQ=255)
                                                                               cGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACACCTTcc        <  1:333966/96‑1 (MQ=255)
                                                                               cGGGTTAGCGGTGCCGACGATCACTTCAACCCGCGAAGGTCCCTCGGCAAACCACAGCGCGGGCAGCCCCGTTTGCAGCACCAGCGTACAACTTcc        <  1:6494080/96‑1 (MQ=255)
                                                                                gggTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCg       >  1:6999745/1‑96 (MQ=255)
                                                                                gggTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCg       <  1:10012440/96‑1 (MQ=255)
                                                                                 ggTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCgg      >  2:8786360/1‑96 (MQ=255)
                                                                                  gTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGgc     <  2:4909749/96‑1 (MQ=255)
                                                                                  gTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGgc     <  2:3714954/96‑1 (MQ=255)
                                                                                    tATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGgcgc   >  2:4032178/1‑96 (MQ=255)
                                                                                                |                                                                                     
GCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGAGGCGGCGCCGACGGGTTATCGGTCGCCGCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTTGCAGCACCAGCGTACAACTTCCGGCGCT  >  NC_000913/3554023‑3554202

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