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.7% 979.6 / ‑4.7 291coding (2/624 nt)yhgKJRNA phosphage cyclase
Reads supporting (aligned to +/- strand):  ref base . (0/1);  new base C (138/152);  total (138/153)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 338 total aligned reads displayed.

CGGCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGAGGCGGCGCCGACGGGTTATCGGTCGCCGCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGTTT  >  NC_000913/3553996‑3554174
                                                                                             |                                                                                        
 ggCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgc                                                                                          <  1:5515490/93‑2 (MQ=255)
 ggCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgc                                                                                          <  1:7050424/93‑2 (MQ=255)
  gCAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgca                                                                                         <  1:1584935/93‑3 (MQ=255)
   cAGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgcag                                                                                        <  2:5068790/93‑4 (MQ=255)
    aGGGTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCgcagg                                                                                       <  1:604345/93‑5 (MQ=255)
      ggTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcg                                                                                    >  1:2244064/1‑93 (MQ=255)
      ggTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcg                                                                                    <  2:387949/93‑1 (MQ=255)
       gTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcgg                                                                                   <  2:2978734/93‑1 (MQ=255)
       gTAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcgg                                                                                   <  1:4034372/93‑1 (MQ=255)
        tAAAAACCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCAGCGGATAAAATCCGCAggcggg                                                                                  >  2:3500024/1‑92 (MQ=255)
         aaaaaCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAggcgggg                                                                                 >  2:2238973/1‑91 (MQ=255)
         aaaaaCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGgcggcg                                                                                 <  2:8776550/93‑1 (MQ=255)
         aaaaaCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGgcggcg                                                                                 <  1:3708995/93‑1 (MQ=255)
          aaaaCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGgcgc                                                                                <  2:3308136/93‑1 (MQ=255)
            aaCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCg                                                                              <  1:2074288/93‑1 (MQ=255)
            aaCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCg                                                                              >  1:6647311/1‑93 (MQ=255)
             aCCGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGa                                                                             <  2:3047134/93‑1 (MQ=255)
              ccGTGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGGTCCCGCACCCAGCGGATAAACTCCGCCGGCGGTGCCGAc                                                                            >  2:1429244/1‑93 (MQ=255)
                 tGGCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACggg                                                                         >  2:7873945/1‑93 (MQ=255)
                  ggCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGt                                                                        <  1:2785825/93‑1 (MQ=255)
                  ggCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGt                                                                        >  1:2343660/1‑93 (MQ=255)
                   gCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGtt                                                                       <  2:7706054/93‑1 (MQ=255)
                   gCGTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGtt                                                                       <  1:1378410/93‑1 (MQ=255)
                     gTAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTAt                                                                     >  2:8242851/1‑93 (MQ=255)
                      tAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATc                                                                    <  2:7347308/93‑1 (MQ=255)
                      tAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATc                                                                    >  2:4024105/1‑93 (MQ=255)
                      tAACAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATc                                                                    >  1:6955705/1‑93 (MQ=255)
                       aaCAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCg                                                                   >  2:6401684/1‑93 (MQ=255)
                        aCAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCgg                                                                  <  2:8528005/93‑1 (MQ=255)
                        aCAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCgg                                                                  >  2:5119027/1‑93 (MQ=255)
                        aCAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCCGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCgg                                                                  >  1:4556832/1‑93 (MQ=255)
                         cAGCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGGCGACGGGCTATCGGt                                                                 >  1:314489/1‑93 (MQ=255)
                          atcGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGCCGGATAAAATCCGCAGGCGGCGCCGACCGGTCATCCGTc                                                                >  2:5524086/3‑93 (MQ=255)
                           gCGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgc                                                               <  1:4130507/93‑3 (MQ=255)
                            cGTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgcc                                                              <  1:7406276/93‑4 (MQ=255)
                             gTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgccg                                                             >  1:3951556/1‑89 (MQ=255)
                             gTGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgccg                                                             >  2:1960370/1‑89 (MQ=255)
                              tGGTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCCGCGGATAAAATCCGCAGGCCGCGCCCAAGGGCTATCGGTggcgc                                                            >  2:3642937/1‑91 (MQ=255)
                               ggTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgccgcc                                                         <  1:4963481/93‑2 (MQ=255)
                               ggTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgccgcc                                                         <  1:8134876/93‑2 (MQ=255)
                                gTTTGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgccgccg                                                        <  1:5370990/93‑3 (MQ=255)
                                 tttGCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGtgccgccgc                                                       >  2:7395790/1‑90 (MQ=255)
                                    gCTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCTGCGGATAAAATCCGCAGGCGGCGCCGACGGGTGATCTGTGCCGCCGCTCa                                                     >  2:5387494/1‑93 (MQ=255)
                                     cTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCAc                                                    <  1:8500144/93‑1 (MQ=255)
                                     cTGATGAATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCAc                                                    <  1:8299627/93‑1 (MQ=255)
                                       gatgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACtt                                                  >  2:2598615/1‑93 (MQ=255)
                                       gatgaATTCCTATTTTCGCCAGCAGCAGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACtt                                                  >  2:4892558/1‑93 (MQ=255)
                                        atgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTc                                                 >  1:2719794/1‑93 (MQ=255)
                                        atgaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTc                                                 <  2:1702522/93‑1 (MQ=255)
                                          taaTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCaa                                               <  1:7948571/92‑1 (MQ=255)
                                          gaATTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCaa                                               <  2:6691527/93‑1 (MQ=255)
                                           aaTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAac                                              <  2:3275996/93‑1 (MQ=255)
                                           aaTTCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAac                                              >  2:4685155/1‑93 (MQ=255)
                                             ttCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAacac                                            >  2:1042314/1‑93 (MQ=255)
                                              tCCTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACAcg                                           <  2:7822293/93‑1 (MQ=255)
                                               ccTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAacacac                                          >  1:3387619/1‑91 (MQ=255)
                                                cTATTTTCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACAcgcg                                         >  2:6328227/1‑93 (MQ=255)
                                                  aTTTTCGCCAGCAGCGGCTCCAGCACCCGGTGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGaa                                       <  2:6122591/93‑1 (MQ=255)
                                                   ttttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGCAg                                      >  1:6097849/1‑93 (MQ=255)
                                                   ttttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAg                                      <  2:136871/93‑1 (MQ=255)
                                                   ttttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAg                                      <  1:1665635/93‑1 (MQ=255)
                                                   ttttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAg                                      >  2:5958505/1‑93 (MQ=255)
                                                    tttCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAgg                                     <  1:6440841/93‑1 (MQ=255)
                                                      tCGCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTc                                   <  1:5163041/93‑1 (MQ=255)
                                                        gCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCa                                 >  1:3846418/1‑93 (MQ=255)
                                                        gCCAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCa                                 >  2:6814897/1‑93 (MQ=255)
                                                         ccAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCAt                                <  1:5800973/93‑1 (MQ=255)
                                                         ccAGCAGCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCAt                                >  1:4864892/1‑93 (MQ=255)
                                                         ccAGCAGCGGCCCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCAt                                <  1:6310807/93‑1 (MQ=255)
                                                          cagcagCGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATc                               <  2:4557144/93‑1 (MQ=255)
                                                           agcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCg                              <  2:3918935/93‑1 (MQ=255)
                                                           agcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCg                              >  1:772915/1‑93 (MQ=255)
                                                            gcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCgg                             >  1:6345286/1‑93 (MQ=255)
                                                            gcagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCgc                             >  1:6209951/1‑92 (MQ=255)
                                                             cagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGc                            >  2:418381/1‑93 (MQ=255)
                                                             cagcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGc                            <  2:4323097/93‑1 (MQ=255)
                                                               gcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaa                          <  1:1899675/93‑1 (MQ=255)
                                                               gcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaa                          >  1:256112/1‑93 (MQ=255)
                                                               gcGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaa                          >  2:6238588/1‑93 (MQ=255)
                                                                cGGCTCCAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCaaa                         <  2:8027147/93‑1 (MQ=255)
                                                                   gtCCAGCACCCGGCGGATAAAATCCGGGTGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACca                      <  1:2394313/92‑1 (MQ=255)
                                                                    tCCACCACCAGGGGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACcac                     <  2:6499363/93‑1 (MQ=255)
                                                                      cAGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAg                   <  1:565212/93‑1 (MQ=255)
                                                                      cAGCACCCGGCGGATAAAATCCGCAGGCCGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGTAGGCCCATCGGCAAACCACAg                   >  2:5731311/1‑93 (MQ=255)
                                                                       aGCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgc                  <  1:1139978/93‑1 (MQ=255)
                                                                        gCACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgcg                 <  2:2368609/93‑1 (MQ=255)
                                                                         gaCCCGGCGGATAAAATCCGCAGCCGTCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgcgc                <  2:2919873/92‑1 (MQ=255)
                                                                         cACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgcgc                <  2:7029392/93‑1 (MQ=255)
                                                                         cACCCGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAgcgc                <  2:7492288/93‑1 (MQ=255)
                                                                           cccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCgg              >  1:7617016/1‑93 (MQ=255)
                                                                           cccGGCGGATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCgg              >  1:442373/1‑93 (MQ=255)
                                                                             cggcggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgc            <  1:1048613/93‑1 (MQ=255)
                                                                                 ggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagca        <  2:1788740/93‑1 (MQ=255)
                                                                                 ggATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGgcagca        >  1:804456/1‑93 (MQ=255)
                                                                                  gATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc       <  2:715747/93‑1 (MQ=255)
                                                                                  gATAAAATCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc       >  1:8365720/1‑93 (MQ=255)
                                                                                  gATAAAATCCGCAGGCGGCGCCGACGGGTTATCCGTGCCGCCGCTCACTTCGACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCAc       >  2:1250486/1‑93 (MQ=255)
                                                                                      aaaTCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTAACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGtt   <  2:5290677/93‑1 (MQ=255)
                                                                                       aaTCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGttt  <  2:1093667/93‑1 (MQ=255)
                                                                                       aaTCCGCAGGCGGCGCCGACGGGTTATCGGTGCCGCCGCTCACTTCAACACGCGAAGGTCCATCGGCAAACCACAGCGCGGGCAGCACCGttt  >  1:6755622/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 < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: