Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,542,076 C→T 12.9% intergenic (+56/+47) nadE → / ← PP_4870 NH(3)‑dependent NAD(+) synthetase/azurin

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,542,0760CT12.9% 417.7 / 47.3 164intergenic (+56/+47)nadE/PP_4870NH(3)‑dependent NAD(+) synthetase/azurin
Reads supporting (aligned to +/- strand):  ref base C (69/71);  new base T (7/14);  total (76/88)
Fisher's exact test for biased strand distribution p-value = 2.41e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.48e-01

Only 100 of 165 total aligned reads displayed.

GCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCGAACTT  >  NC_002947/5541927‑5542224
                                                                                                                                                     |                                                                                                                                                    
gCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCCACGc                                                                                                                                                      >  2:1374341/1‑150 (MQ=255)
    gggCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACCCaatg                                                                                                                                                  >  2:1571746/1‑147 (MQ=255)
      gCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTcgc                                                                                                                                                <  2:5889380/150‑1 (MQ=255)
      gCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTcgc                                                                                                                                                <  2:5549956/150‑1 (MQ=255)
        aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTCCCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATACGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                              <  2:4038140/150‑1 (MQ=255)
        aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                              <  2:1717864/150‑1 (MQ=25)
        aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                              <  2:522313/150‑1 (MQ=25)
        aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                              <  2:5334414/150‑1 (MQ=25)
        aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                              <  2:2295221/150‑1 (MQ=25)
        aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                              <  2:2294925/150‑1 (MQ=25)
        aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                              <  2:655004/150‑1 (MQ=255)
        aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                              <  2:1689240/150‑1 (MQ=255)
        aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                              <  2:4370029/150‑1 (MQ=255)
        aGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                              <  1:4956474/150‑1 (MQ=255)
          ccGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGt                                                                                                                                            <  2:4528426/150‑1 (MQ=255)
             gTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGaa                                                                                                                                         <  1:5476027/150‑1 (MQ=25)
             gTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                              >  1:4270111/1‑145 (MQ=255)
                   cAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCCCCAATAAGAGGGTGGACCCACCCACGCAATAGCGGGTGAAacccca                                                                                                                                   >  1:5968311/1‑147 (MQ=255)
                    aGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCCAATACGAGGGTGAAACCACCCACCCCATCGCGGGGGAccccccac                                                                                                                                  >  2:5288232/1‑146 (MQ=255)
                    aGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCCCGAAAACGAGGGTGAAACCAACCACGCAATCGCGGGTGAAcccccgc                                                                                                                                  >  2:184499/1‑146 (MQ=255)
                     ggCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACCAGGGTGAACCCACCGACGCATTCGCGGGTGAAACCGCTcc                                                                                                                                 >  2:2374291/1‑150 (MQ=255)
                      gCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCt                                                                                                                                <  2:6059133/150‑1 (MQ=255)
                           gTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTCACCCGCGCAAACGAGGGTGAAACCACCCACGCAATAGCGGGTGAACCCGCACCTGCacg                                                                                                                           >  2:2761148/1‑148 (MQ=255)
                           gTTCGACACCATCGCCGCCACCCCCCGAACCACCCAGCACACGCGCGAACTGCGGTTCGCTCCTTTTTACGGGGCCTGCTGGAGGGGTTTCACCCGCTAATACTAGGTTGAACCCACCGACGAATTCGCGGGTGAACCCGCTCCTGCAgg                                                                                                                           <  2:4436251/150‑1 (MQ=255)
                               gACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGTACCCACCGACGCCGTCGCGGGTGAACCCGCTCCGGCAGGGttg                                                                                                                       >  2:3635366/1‑149 (MQ=255)
                               cacacCATCGCCGCCACCAACCGCACGACCCAGCACACGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGCCGCGGGTTCACCCGCTAATACGCGGGTGAACCCACCGACGCATCGGCGGGTGACCCCGCTCCTGCAGGGTTc                                                                                                                       <  2:1645715/149‑1 (MQ=255)
                                   tcatcGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTGACCCGCGCATAAGAGGGTGGACCCCCCCACGCCTTGGCGGGTGTAACCGCGCCGGCAggggtgggcg                                                                                                                   >  1:4546799/1‑143 (MQ=255)
                                     atcGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTCACCCCCGCATACGGGGGTGAAACCACCGGCGCATTCGCGGGTGATCCCGCGCCTGCCGGGTTGGCCggc                                                                                                                 >  1:4171951/1‑148 (MQ=255)
                                      tggtcgCCACCTACCCCACGACCCCCCACACGCGCTACGTGCGGTTGGCTCCTTTTTATGGGGCCTGCCGGCGGGGGTTCCCCGGCGAATCCGCGGGTCAACCCCCCGCCGCTTGGGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCa                                                                                                                <  2:123181/148‑1 (MQ=255)
                                        gtcgCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCCCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                              >  2:254552/1‑118 (MQ=18)
                                           gCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                              <  2:5087222/115‑1 (MQ=18)
                                            ccaccaACCCCAACACCCACCACAACCGCCAACTGCCGTTCGCTCCTTTATAAGTGGCCTGCAGCAGCGGTTTCACCCGCGAAAACGAGGGTCAACCAACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCatt                                                                                                          <  2:5361157/150‑1 (MQ=37)
                                            ccACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCGTGATAAGTGGCCTGCAGGAGCGGGGTCACCCGCGCAAACGAGGGGGGACCCCCCCCCGCAATCGCGGGTGGACCCCCTCCGGGAGGGGTCGCCGGCAGCCatt                                                                                                          >  2:3985134/1‑150 (MQ=255)
                                                cTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCgg                                                                                                                                              >  1:3117100/1‑110 (MQ=255)
                                                 tACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTGACCCCCGAATAACAGGGTGAACCCACCGACGCATTCGCGGGTGTACCCGCGCCTGCAGGGGTCGCCCGCAGCCGTTGATTa                                                                                                     >  2:3445653/1‑150 (MQ=37)
                                                 tACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCgg                                                                                                                                              >  1:3980811/1‑109 (MQ=17)
                                                       caccccccGCACACGCGCGACCGGCCGTTCGCTCCTTTATATGTGGCCTGCAGGAGGGGGTTCACCCGCGAATACGTGGGTGAACCCACCGACGCATTCGCGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGAc                                                                                               <  2:728528/146‑1 (MQ=255)
                                                       caGACCCAGCACCCGCGCGAACTGCCGTTCGCTCCTTGTTAAGTGGCCTGCAGGAGCGGGTTCCCCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGAc                                                                                               <  2:4614210/149‑1 (MQ=255)
                                                             ccgCACACGCGCGACGTGCGTTTCGCTCCTTGTTACGTGGCCTGCCGGCGGGGTTCCACCCGCGAATACGCGGGTGAACCCACCGCCGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGAc                                                                                         <  1:2375413/148‑1 (MQ=25)
                                                               gCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTGCCCCCCGAATCAGCGGGTGGACCCACCGACGCAGTCGCGGGTGGACCCGCTGCGGCCGGGTGCGCCGTCAGCCGTGGCTTACTTGTCGCTGGccg                                                                                       >  2:7160398/1‑149 (MQ=255)
                                                                cacaAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGGTCACCCGCGAATACGAGGGTGTACCCACCGACGCAGTCGCGGGTGTACCCGCTCCTGCAGGGTTGGCCCTCAGCCTTTGGTTTCTTGGCGCTgtgctt                                                                                      >  2:3192185/1‑145 (MQ=255)
                                                                  caAGCGCGAACGGCCGTTCGCTCCTTTATATGTGGCCTGCAGGAGGGGGTTCACCCGCGAATACGAGGGTTAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGc                                                                                    <  2:6366829/150‑1 (MQ=38)
                                                                     gcccccaCTGCCGTTCCCTCATTTATATGGGGCCGGCTGGAGCGGTTCCACCCGCCAATACGCGGGTGACCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCtt                                                                                 <  1:1129412/144‑1 (MQ=25)
                                                                       gcgAACTGCCTTTCGCTCCTTGATAAGTGGCCTGCTGGAGGGGGTTCACCCGCGAATACGTGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtc                                                                               <  1:7347620/150‑1 (MQ=255)
                                                                       gcgAACTGCCGTTCGCTCCTTTATATGTGGCCTGCTGGAGCGGGTTCACCCGCTAATACGCGGGTTAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTtc                                                                               <  2:3393853/150‑1 (MQ=39)
                                                                              gCCGTTCGCTCCTTGACACGTGGCCTGCAGGAGCGGGTACACCCGCGACTCCGAGGGTGAACCCACCGCCGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCATTGATTACTTGACGATGACCTTGCCTTTCATCATCg                                                                        <  2:4270051/150‑1 (MQ=38)
                                                                               ccGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTAACCCCCCAAAAAGATGGTTAAACCAACCACGCAATAGCGGGTGAACCCCCACCGGCAGGGTTCGCCGGCACCCATTTATTAAATTAAGATGAACATTCCATTAATCATAga                                                                       >  2:3495908/1‑150 (MQ=255)
                                                                                cGTTCGCTCCTTGATACGTGGCCTGCAGGAGCGGGTTCACCCGCTAATACGCGGGTGAACCCACCGACGCATGCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCgag                                                                      <  2:2269096/150‑1 (MQ=255)
                                                                                           ggATAAGTGGCCGGCAGTAGCGGGGTCAGCCGCTATTACTAGGGTGAACCCACCGAGGCATGGGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGGCTTTGATTAGTTGACGATGACCTTGCCTTTCATCTTCGAGATGTGGCCCgg                                                           <  2:4759394/149‑1 (MQ=255)
                                                                                               ccggggCCTGCTGGCGCGGGTTCACCCGGTAATCCGGGGGGGACCCCCCGGCCGCATTGGGGGGTGACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAg                                                       <  1:4913335/146‑1 (MQ=255)
                                                                                                cggggccggccGGCGCGGGTTCCCCCGCGAATCCGCGGGTGCCCCCACCGACGCCGGCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAgg                                                      <  2:4678995/147‑1 (MQ=255)
                                                                                                cgccgccggccgGCGCGGGGCCCCCCGCGACTCCGCGGGTGCACCCACCGCCGCAGGCGCGGGGGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCCTCATCGAGAGGTGGCCCGGGAAgg                                                      <  2:5724388/146‑1 (MQ=255)
                                                                                                      cTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCCCCGCCGCAATCGCGGGTGAACCCCCTACTGCAGGGGTCGCCGTCGGCCTGTTATTACTATACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAATGTGCaga                                                >  2:4310640/1‑150 (MQ=255)
                                                                                                           ggcgggggTTCACCCGCTAATACGTGGTTTACCCCACCGACGCATTCGGGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGaagaag                                           <  1:2305563/145‑1 (MQ=255)
                                                                                                              gggggTTCACCCGCGAATCCGCGGGTGACCCCACCGACGCATGCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTg                                        <  1:6756091/148‑1 (MQ=255)
                                                                                                              gCGGGTACACCCGCGAATACGAGGGTCAACCCACCGACGAATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTg                                        <  2:1835592/150‑1 (MQ=255)
                                                                                                               cGGGTTCACCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGt                                       >  2:2724962/1‑150 (MQ=21)
                                                                                                                  gTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt                                    >  2:929464/1‑150 (MQ=255)
                                                                                                                  gTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGt                                    >  2:4657484/1‑150 (MQ=255)
                                                                                                                      aCCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGaagaag                                           >  1:1102443/1‑139 (MQ=21)
                                                                                                                      aCCCGCGAATGCGTCGGTGGGTTCACCCTCGTATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacc                                >  2:5134414/1‑150 (MQ=25)
                                                                                                                      aCCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcacc                                >  2:6134975/1‑150 (MQ=255)
                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCa                                                                                                                             >  2:5174832/1‑55 (MQ=255)
                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCa                                                                                                                             <  1:5174832/55‑1 (MQ=255)
                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAggc                                                                                                                          >  2:7184210/1‑57 (MQ=255)
                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCCTCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                              >  2:1085133/1‑150 (MQ=255)
                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTa                                      >  2:3141033/1‑142 (MQ=255)
                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                              >  1:3566633/1‑150 (MQ=255)
                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                              >  1:7050781/1‑150 (MQ=255)
                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                              >  1:1824776/1‑150 (MQ=255)
                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                              >  1:7329282/1‑150 (MQ=255)
                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCATTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                              >  2:3350257/1‑150 (MQ=255)
                                                                                                                        ccGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGACTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTcaccac                              >  1:2428556/1‑150 (MQ=255)
                                                                                                                          gctaaTCCGTGGGTGCCCCCCCCGCCGCATGGGCGGGGGCACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCg                            <  1:3812529/147‑1 (MQ=255)
                                                                                                                          gcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCg                            >  1:6424253/1‑150 (MQ=255)
                                                                                                                          gcgAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCg                            >  1:5278640/1‑150 (MQ=255)
                                                                                                                              aTACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCtt                        <  2:5568306/150‑1 (MQ=255)
                                                                                                                               tACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCCTCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTc                       <  2:6897519/150‑1 (MQ=255)
                                                                                                                                  gcgggTGAACCCACCTACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAg                    <  2:3717102/148‑1 (MQ=255)
                                                                                                                                  gAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGACGGTGCAGGAGAAGCTGTAGTCACCACCGGCTTCCAg                    >  2:7032291/1‑150 (MQ=255)
                                                                                                                                  gAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAg                    >  1:5508061/1‑150 (MQ=255)
                                                                                                                                  gAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAg                    >  2:4144575/1‑150 (MQ=255)
                                                                                                                                  gAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAg                    >  2:3918958/1‑150 (MQ=255)
                                                                                                                                     ggtgccccccccGCCGCATGCGGGGGTCAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCtt                 <  1:1571746/144‑1 (MQ=255)
                                                                                                                                     ggTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCtt                 >  2:1336504/1‑150 (MQ=255)
                                                                                                                                      gTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTg                >  2:5017221/1‑150 (MQ=255)
                                                                                                                                       ggaccccacccaCGCATCCGGGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTgg               <  1:2761148/146‑1 (MQ=255)
                                                                                                                                        gAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGa              <  2:1520331/150‑1 (MQ=255)
                                                                                                                                         aaCCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGaa             <  2:4162532/150‑1 (MQ=255)
                                                                                                                                         aaCCCACCCACGCATTCGCGGTTTACCCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGaa             <  1:1476290/150‑1 (MQ=255)
                                                                                                                                          gttcACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAg            >  2:6154288/4‑150 (MQ=255)
                                                                                                                                           cccACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGATGTAGTCACCACCGGCTTCCAGCTTGGAAgt           >  2:844128/1‑150 (MQ=255)
                                                                                                                                           cccACCGACGCATTCGCGGGTGAACCAGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAgt           >  2:844692/1‑150 (MQ=255)
                                                                                                                                             cACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATAATAGATATGTGGCCAGGGAATGTGCAGAAGAAGCTGTAGTCACCACCGGCTTACAGCTTGGAAgtgt         >  1:5619176/1‑150 (MQ=255)
                                                                                                                                              aCCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTc        >  2:5643396/1‑150 (MQ=255)
                                                                                                                                              aCCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTc        <  1:612642/150‑1 (MQ=255)
                                                                                                                                                cGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCGa      >  2:4509284/1‑150 (MQ=255)
                                                                                                                                                  aCGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGACGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCGAac    >  1:3852419/1‑150 (MQ=255)
                                                                                                                                                    gCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCGAactt  <  1:488346/150‑1 (MQ=255)
                                                                                                                                                    gCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCGAactt  >  1:7277720/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GCATGGGCAGCCGGTTGACCAGGCAGCGTTCGACATCATCGTCGCCACCTACCGCAAGACCCAGCACAAGCGCGAACTGCCGTTCGCTCCTTGATAAGTGGCCTGCAGGAGCGGGTTCACCCGCGAATACGAGGGTGAACCCACCGACGCATTCGCGGGTGAACCCGCTCCTGCAGGGTTCGCCGTCAGCCTTTGATTACTTGACGATGACCTTGCCTTTCATCATCGAGATGTGGCCCGGGAAGGTGCAGAAGAAGCTGTAGTCACCACCGGCTTCCAGCTTGGAAGTGTCGAACTT  >  NC_002947/5541927‑5542224

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 11 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: