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

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,554,0881.C99.3% 841.9 / ‑1.8 276coding (2/624 nt)yhgKJRNA phosphage cyclase
Reads supporting (aligned to +/- strand):  ref base . (0/2);  new base C (119/153);  total (120/156)
Fisher's exact test for biased strand distribution p-value = 5.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.94e-01

Only 100 of 323 total aligned reads displayed.

CGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGAGGCGGCGCCGACGGGTTATCGGTCGCCGCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTT  >  NC_000913/3553996‑3554174
                                                                                               |                                                                                       
cGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCg                                                                                          >  1:2529814/1‑93 (MQ=255)
  gCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgca                                                                                        <  2:2007471/93‑3 (MQ=255)
  gCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgca                                                                                        >  1:1571286/1‑91 (MQ=255)
  gCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgca                                                                                        <  1:3131097/93‑3 (MQ=255)
  gCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgca                                                                                        <  2:5199264/93‑3 (MQ=255)
   cAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgcag                                                                                       <  1:3481472/93‑4 (MQ=255)
   cAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgcag                                                                                       <  2:3791241/93‑4 (MQ=255)
   cAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgcag                                                                                       >  1:996005/1‑90 (MQ=255)
   cAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgcag                                                                                       <  1:7639815/93‑4 (MQ=255)
    aGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgcagg                                                                                      <  1:4637534/93‑5 (MQ=255)
    aGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgcagg                                                                                      >  1:2067689/1‑89 (MQ=255)
    aGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgcagg                                                                                      >  1:7867292/1‑89 (MQ=255)
    aGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgcagg                                                                                      >  1:629414/1‑89 (MQ=255)
      ggTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcg                                                                                    >  2:6616310/1‑93 (MQ=255)
       gTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcgg                                                                                   >  1:818687/1‑93 (MQ=255)
        tAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcggg                                                                                  >  2:4006947/1‑92 (MQ=255)
          aaaaCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGgcgc                                                                                <  1:4847479/93‑1 (MQ=255)
            aaCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGCTAAAATCCGCAGGCGGCGCCg                                                                              >  2:5891119/1‑93 (MQ=255)
             aCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGa                                                                             <  1:3322931/93‑1 (MQ=255)
              ccGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAc                                                                            >  1:4063635/1‑93 (MQ=255)
              ccGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGAc                                                                            >  2:4364937/1‑93 (MQ=255)
                gTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACgg                                                                          >  1:5875338/1‑93 (MQ=255)
                gTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACgg                                                                          <  2:7108558/93‑1 (MQ=255)
                gTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACgg                                                                          >  1:1862388/1‑93 (MQ=255)
                 tGGCGTANCAGCGTGGTTTGCNGNTGNATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACggg                                                                         <  1:3640573/93‑1 (MQ=255)
                 tGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGCTAAAATCCGTAGGCGGCGCCGACggg                                                                         >  1:2626179/1‑93 (MQ=255)
                 tGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACggg                                                                         <  1:2279115/93‑1 (MQ=255)
                    cGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTa                                                                      <  2:4912845/93‑1 (MQ=255)
                      tAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATc                                                                    <  1:2985458/93‑1 (MQ=255)
                        aCAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCgg                                                                  >  2:3105270/1‑93 (MQ=255)
                        aCAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCgg                                                                  <  2:2382767/93‑1 (MQ=255)
                        aCAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGAGGATAAAATCCGCAGGCGGCGCCGATGGGTTATCgg                                                                  >  2:4183465/1‑93 (MQ=255)
                         cAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCggg                                                                 >  2:2828084/1‑92 (MQ=255)
                         cAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGt                                                                 <  1:3920346/93‑1 (MQ=255)
                          aGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtg                                                                <  2:6125962/93‑2 (MQ=255)
                           gCGTGGTTTGCTGATGAATTCCTATTTTCCCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgc                                                               <  2:3418857/93‑3 (MQ=255)
                            cGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgcc                                                              <  2:5952734/93‑4 (MQ=255)
                            cGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgcc                                                              <  1:342853/93‑4 (MQ=255)
                             gTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgccg                                                             >  1:6893634/1‑89 (MQ=255)
                              tGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGtgccgc                                                          >  2:3652426/1‑93 (MQ=255)
                               ggTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgccgcc                                                         >  1:5820635/1‑92 (MQ=255)
                                 tttGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgccgccgc                                                       <  1:5705079/93‑4 (MQ=255)
                                  ttGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACCGGGTATCGGtgccgccgct                                                      >  1:3053362/1‑89 (MQ=255)
                                      tgatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACt                                                   >  2:3329543/1‑93 (MQ=255)
                                      tgatgaATTCCTATTTTCGCCAGCAGCGGCTCCACCACCCGGCGGCTAAAATGCGCAGGCAGCGCCCACGGGTTATCGGTGCCGCCGCTCACt                                                   >  2:1355230/1‑93 (MQ=255)
                                       gatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACtt                                                  >  1:6977695/1‑93 (MQ=255)
                                       gatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACtt                                                  >  2:545271/1‑93 (MQ=255)
                                       gatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACtt                                                  >  2:530499/1‑93 (MQ=255)
                                       gatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACtc                                                  >  1:1440177/1‑92 (MQ=255)
                                       gatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGGGCCGCCGCTCACtt                                                  >  2:6837683/1‑93 (MQ=255)
                                       gatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGAGCCGACGGGTTATCGGTGCCGCCGCTCACtt                                                  >  2:1853502/1‑93 (MQ=255)
                                        atgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTc                                                 <  2:6971299/93‑1 (MQ=255)
                                        atgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGTGCCGCCGCTCACTTc                                                 >  2:6730897/1‑93 (MQ=255)
                                           aaTTCCTATGTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAac                                              <  1:1914806/93‑1 (MQ=255)
                                            aTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAaca                                             <  2:756923/93‑1 (MQ=255)
                                              tCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACAcg                                           >  1:2789640/1‑93 (MQ=255)
                                               ccTATTTTCGCCAGCAGCGGCTCCAGTCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACAcg                                           <  2:6347659/93‑1 (MQ=255)
                                                cTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACAcgcg                                         >  1:4324470/1‑93 (MQ=255)
                                                 tATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGa                                        <  2:4628966/93‑1 (MQ=255)
                                                   ttttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAg                                      >  2:562508/1‑93 (MQ=255)
                                                   ttttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAg                                      <  2:7653509/93‑1 (MQ=255)
                                                   ttttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAg                                      >  1:2178055/1‑93 (MQ=255)
                                                   ttttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAg                                      >  1:2501137/1‑93 (MQ=255)
                                                   ttttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAg                                      >  2:7380892/1‑93 (MQ=255)
                                                   ttttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGGTATCGGTGCCGCCGCTCACTTCAACACGCGAAg                                      >  1:1773624/1‑93 (MQ=255)
                                                     ttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGt                                    <  2:1164765/93‑1 (MQ=255)
                                                     ttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGt                                    <  1:641154/93‑1 (MQ=255)
                                                      tCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTc                                   <  1:4915411/93‑1 (MQ=255)
                                                       cGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTcc                                  <  2:1261524/93‑1 (MQ=255)
                                                        gCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCa                                 <  1:3766266/93‑1 (MQ=255)
                                                        gCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCa                                 <  2:6170089/93‑1 (MQ=255)
                                                         ccAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTTAACACGCGAAGGTCCAt                                >  2:479386/1‑93 (MQ=255)
                                                         ccAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCAt                                >  1:5988828/1‑93 (MQ=255)
                                                         ccAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCAt                                <  1:4455036/93‑1 (MQ=255)
                                                          cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATc                               <  2:2270082/93‑1 (MQ=255)
                                                          cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATc                               <  1:5178464/93‑1 (MQ=255)
                                                          cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATc                               <  1:7464618/93‑1 (MQ=255)
                                                          cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATc                               >  2:1423817/1‑93 (MQ=255)
                                                          cagcagCGCCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATc                               <  2:3028899/93‑1 (MQ=255)
                                                           agcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCg                              <  2:6892954/93‑1 (MQ=255)
                                                           agcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCg                              <  2:2724092/93‑1 (MQ=255)
                                                            gcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCgg                             >  2:1954379/1‑93 (MQ=255)
                                                                  gCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAAcc                       >  1:1280686/1‑93 (MQ=255)
                                                                   cTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACca                      <  2:494389/93‑1 (MQ=255)
                                                                   cTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACca                      <  2:392357/93‑1 (MQ=255)
                                                                    tCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcac                     >  2:1642190/1‑93 (MQ=255)
                                                                      cAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAg                   <  2:778142/93‑1 (MQ=255)
                                                                      cAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAg                   <  1:1967168/93‑1 (MQ=255)
                                                                      cAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAg                   <  1:6736135/93‑1 (MQ=255)
                                                                      cAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAg                   <  2:4246956/93‑1 (MQ=255)
                                                                       aGCACCCGGCGGATAAAATCCNCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgc                  <  1:1063808/93‑1 (MQ=255)
                                                                          cgcccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgcgc                <  2:2216569/91‑1 (MQ=255)
                                                                           aCCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGcgcg               >  1:6159226/1‑93 (MQ=255)
                                                                            cccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCgg              <  1:5243141/93‑1 (MQ=255)
                                                                            cccGACGGATAAAATACACAGGCGGCGCCGACGGGTTATCGGAACCGCCGCTCACTTCAACACGCGAAGGTCCATCAGCAAACCACAGcacga              <  2:7037722/93‑2 (MQ=255)
                                                                              cggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgc            <  2:6249280/93‑1 (MQ=255)
                                                                               ggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgca           >  2:3577745/1‑93 (MQ=255)
                                                                               ggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgca           >  1:4214801/1‑93 (MQ=255)
                                                                                   gATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc       >  1:7816421/1‑93 (MQ=255)
                                                                                    aTAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAcc      >  2:5109746/1‑93 (MQ=255)
                                                                                               |                                                                                       
CGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGAGGCGGCGCCGACGGGTTATCGGTCGCCGCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTT  >  NC_000913/3553996‑3554174

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