Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,674,777 A→G 6.4% intergenic (+139/+1) rhlE‑II → / ← yceI ATP‑dependent RNA helicase/polyisoprene‑binding acidic stress response factor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,674,7770AG 6.4% 184.4 / 2.8 62intergenic (+139/+1)rhlE‑II/yceIATP‑dependent RNA helicase/polyisoprene‑binding acidic stress response factor
Reads supporting (aligned to +/- strand):  ref base A (29/29);  new base G (2/2);  total (31/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.06e-01

GGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑C‑GCGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCA  >  NC_002947/5674632‑5674925
                                                                                                                                                   |                                                                                                                                                     
gggCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAAtt                                                                                                                                                     <  1:49032/149‑1 (MQ=255)
           aGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAGTAATTTCTGccag                                                                                                                                          <  7:287235/149‑4 (MQ=25)
            gTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTca                                                                                                                                         >  5:280275/1‑149 (MQ=255)
                  aGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCtt                                                                                                                                   >  7:55864/1‑149 (MQ=255)
                   gggggCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTc                                                                                                                                  <  5:297260/149‑1 (MQ=255)
                   gggggCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTc                                                                                                                                  <  4:284489/149‑1 (MQ=255)
                      ggCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGac                                                                                                                               <  6:100678/149‑1 (MQ=255)
                                        tCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAGTG‑GCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccaggtccag                                                                                                             >  4:149719/1‑149 (MQ=255)
                                        tCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccaggtcca                                                                                                              <  8:55864/148‑1 (MQ=255)
                                          aCGACACAAGGCCGCTCCTTCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccaggtccagg                                                                                                            <  7:321262/148‑1 (MQ=255)
                                          aCGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccaggtccagg                                                                                                            <  1:223472/148‑1 (MQ=255)
                                          aCGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATAtccaggtccagg                                                                                                            <  5:212846/148‑1 (MQ=255)
                                          aCGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAggg                                                                                                           >  2:217545/1‑149 (MQ=255)
                                              cacaAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGacac                                                                                                                             >  8:191458/1‑127 (MQ=255)
                                              cacaAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGacac                                                                                                                             <  7:191458/127‑1 (MQ=255)
                                               agaaGGCCGCTCGTGCAGATAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGc                                                                                                      <  4:173851/147‑1 (MQ=38)
                                                  aGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGaa                                                                                                   <  1:241266/149‑1 (MQ=255)
                                                    gCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGt                                                                                                 <  7:157790/149‑1 (MQ=255)
                                                       gCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAg                                                                                               >  7:158392/1‑148 (MQ=255)
                                                         tCCTGCAGAGAATTTCTTGATCTCCACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAGTAATTTCTGccag                                                                                                                                          >  6:314663/1‑100 (MQ=25)
                                                         tCCTGCAGAGAATTTCTTGATCTCCACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAGTAATTTCTGccag                                                                                                                                          <  5:314663/103‑4 (MQ=25)
                                                           cTGCAGAGAATTTCTTGACCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCAccc                                                                                          >  3:107871/1‑149 (MQ=255)
                                                            tGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGTCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGc                                                                                                      <  5:286670/136‑1 (MQ=255)
                                                            tGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGc                                                                                                      >  6:286670/1‑136 (MQ=255)
                                                              cAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCg                                                                                         <  6:280275/147‑1 (MQ=255)
                                                                gagaATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑Gcc                                                                                    >  6:29976/1‑149 (MQ=255)
                                                                       tCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGAt                                                                              <  1:263978/148‑1 (MQ=255)
                                                                          tGATCTCTACAGGAGCGGCCTTGTGTGGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAAATAATTCTGCTACACCCCATCGACACAGAAATCCAGGTCCAGGGTCTACGAAGTAGCACCCGG‑GCCTGTGATGcc                                                                           >  6:275124/1‑148 (MQ=255)
                                                                            acctctACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGaa                                                                        <  4:107871/147‑1 (MQ=255)
                                                                               tctACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAgac                                                                     >  2:254618/1‑147 (MQ=255)
                                                                                      gAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCttt                                                                                 >  7:36149/1‑130 (MQ=255)
                                                                                      gAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCttt                                                                                 <  8:36149/130‑1 (MQ=255)
                                                                                      gAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCgttcagg                                                              >  7:281934/1‑149 (MQ=255)
                                                                                          ggCCTTGTGTCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCgttcaggtt                                                            >  2:47070/1‑147 (MQ=255)
                                                                                                gtgtCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCgttcag                                                               >  2:291382/1‑138 (MQ=255)
                                                                                                 tgtCGTGATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCGTTCAGGTTCAGCgtggt                                                   >  1:276543/1‑149 (MQ=255)
                                                                                                       gATAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCCTTCAGGTTCAGCGTGGTGGTGGc                                             >  2:90398/1‑149 (MQ=255)
                                                                                                         tAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGt                                           <  2:276543/149‑1 (MQ=255)
                                                                                                         tAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGt                                           <  5:29976/149‑1 (MQ=255)
                                                                                                         tAGGGCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGt                                           <  4:227083/149‑1 (MQ=255)
                                                                                                           gggCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTGTGATGCCGAAGTCGGTCAGGTTCAGCGTGGTGGTGGCGtt                                          >  4:124505/1‑148 (MQ=255)
                                                                                                            ggCCGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTg                                         >  2:88113/1‑148 (MQ=255)
                                                                                                              ccGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGGCGCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGaa                                       <  1:254618/149‑1 (MQ=255)
                                                                                                               cGCAAAG‑CGGCCCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGc                                     >  7:294/1‑149 (MQ=255)
                                                                                                                         gcccccTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCt                            >  6:43735/4‑149 (MQ=255)
                                                                                                                         gcccccTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCt                            >  7:134400/4‑149 (MQ=255)
                                                                                                                          cccccTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTg                                         >  4:72032/3‑135 (MQ=255)
                                                                                                                          cccccTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTg                                         <  3:72032/133‑1 (MQ=255)
                                                                                                                          cccccTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTc                           <  5:244864/147‑1 (MQ=255)
                                                                                                                           ccccTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCg                          >  1:32260/2‑149 (MQ=255)
                                                                                                                             ccTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGcc                        <  1:88113/149‑1 (MQ=255)
                                                                                                                             ccTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGcc                        <  3:149719/149‑1 (MQ=255)
                                                                                                                             ccTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGc                         >  4:136336/1‑148 (MQ=255)
                                                                                                                               tCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCAc                      >  7:97115/1‑149 (MQ=255)
                                                                                                                               tCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGAACCCGG‑GCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCAc                      >  7:102005/1‑149 (MQ=255)
                                                                                                                                 ggTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCAcc                     <  8:281934/148‑1 (MQ=255)
                                                                                                                                   ttttCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCa                  >  7:106689/1‑149 (MQ=255)
                                                                                                                                   ttttCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCa                  >  3:251610/1‑149 (MQ=255)
                                                                                                                                       cAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATgg               <  1:47070/148‑1 (MQ=255)
                                                                                                                                       cAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATgg               <  8:294/148‑1 (MQ=255)
                                                                                                                                       cAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATgg               <  1:217545/148‑1 (MQ=255)
                                                                                                                                           aagCAGAAGTAATTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCAGGGGTCCt          >  2:279563/3‑149 (MQ=255)
                                                                                                                                             gcacacaTTTCTGCTGCTTCACCCCTTCGACACTGATCTCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCtt         <  4:251610/142‑1 (MQ=255)
                                                                                                                                                   aTTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCa  >  1:249480/1‑149 (MQ=255)
                                                                                                                                                   aTTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCa  <  3:136336/149‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                     
GGGCTGATAGAAGTACCGAGGGGGCTGCTTCGCAGCCCTATCACGACACAAGGCCGCTCCTGCAGAGAATTTCTTGATCTCTACAGGAGCGGCCTTGTGTCGTGATAGGGCCGCAAAG‑C‑GCGCCCTCGGTTTTCAGCTGGCAGAAATTACTTCTGCTTCACACCTTCGACACTGATATCCAGGTCCAGGGTCTGCGAAGTAGCACCCGG‑GCCTTTGATGCCGAAGTCGTTCAGGTTCAGCGTGGTGGTGGCGTTGAAGCCAGCGCGCTCGCCACCCCATGGGTCCTTGCCTTCA  >  NC_002947/5674632‑5674925

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑