Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 250,911 G→C 6.4% R98G (CGC→GGC)  nfeD ← membrane bound peptidase NefD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947250,9110GC 6.4% 175.8 / 2.1 58R98G (CGC→GGC) nfeDmembrane bound peptidase NefD
Reads supporting (aligned to +/- strand):  ref base G (25/29);  new base C (2/2);  total (27/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 = 9.97e-01

ATCGTCGCCGCCTTTGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑AGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCTTGCGCCTGGGCCTG  >  NC_002947/250765‑251052
                                                                                                                                                    |                                                                                                                                              
aTCGTCGCCGCCTTTGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCgg                                                                                                                                               <  8:12320/148‑1 (MQ=255)
 tcgtcgCCGCCTTTGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGgc                                                                                                                                             <  6:80412/149‑1 (MQ=255)
    tcgCCGCCTTTGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGcgcc                                                                                                                                          <  5:55309/149‑1 (MQ=255)
    tcgCCGCCTTTGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCCGCGCCGCCAGCGCCGGcacc                                                                                                                                          <  6:190207/149‑1 (MQ=255)
         gccTTTGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGTGTAGAGAATATAGGTGCCGGCGCCCGCGCCGCCAGcgccg                                                                                                                                         <  4:102260/145‑1 (MQ=255)
         gccTTTGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGCATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTgg                                                                                                                                     >  6:110741/1‑149 (MQ=255)
         gccTTTGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGC‑‑GGCGCCCAGGTTGGTCACGGGCGCCATCGCTGCAATATGGCTGGGGTAGAGAATATAGGTGCCGGCGCCCGCGCCGCCAGcgccg                                                                                                                                         >  3:102260/1‑145 (MQ=255)
         gccTTTGGCGTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTgg                                                                                                                                     <  3:40061/149‑1 (MQ=255)
           cTTTGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCTGCGCGGGCGCCGCTgggg                                                                                                                                   <  6:350690/149‑1 (MQ=255)
            tttGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTggggg                                                                                                                                  >  2:3959/1‑149 (MQ=255)
                  tttATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGCAGGGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGAAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAAc                                                                                                                              >  1:94397/1‑149 (MQ=255)
                        gTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTgg                                                                                                                      >  4:206433/1‑149 (MQ=255)
                         tCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGc                                                                                                                     >  2:3717/1‑149 (MQ=255)
                             tGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGCCg                                                                                                                 >  6:137878/1‑149 (MQ=255)
                             tGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACg                                                                                                                 >  1:69795/1‑149 (MQ=255)
                               ggCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACgg                                                                                                                <  2:69795/148‑1 (MQ=255)
                               ggCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACgg                                                                                                                <  1:3959/148‑1 (MQ=255)
                                  gTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTAc                                                                                                            >  4:338839/1‑149 (MQ=255)
                                  gTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTAc                                                                                                            >  5:349046/1‑149 (MQ=255)
                                      cAGGCGCGCCGCCGATCCGCACCGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTg                                                                                                        <  2:94397/149‑1 (MQ=255)
                                         gcgcgcCGCCGATCCGCACCGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTg                                                                                                     <  8:135390/149‑1 (MQ=255)
                                           gcgcCGCCGATCCGCACCGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGg‑a                                                                                                  <  8:78890/149‑1 (MQ=255)
                                           gcgcCGCCGATCCGCACCGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGg‑a                                                                                                  <  3:206433/149‑1 (MQ=255)
                                             gccgccGATCCGCACCGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGg‑aga                                                                                                <  6:303213/149‑1 (MQ=255)
                                              ccgccgATCCGCACCGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGACGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGg‑agag                                                                                               <  1:98767/149‑1 (MQ=255)
                                                 ccgATCCGCACCGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑AGAGAAt                                                                                            >  6:125828/1‑149 (MQ=255)
                                                  cgATCCGCACCGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑AGAGAAtt                                                                                           <  1:143204/149‑1 (MQ=255)
                                                        gCACCGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑AGAGAATTGCCtt                                                                                      >  3:234793/1‑148 (MQ=255)
                                                         cACCGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑AGAGAATTGCCTTGa                                                                                    >  6:383202/1‑149 (MQ=255)
                                                         cACCGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑AGAGAATTGCCTTGa                                                                                    >  5:41618/1‑149 (MQ=255)
                                                           ccGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑AGAGAATTGCCTTGatc                                                                                  <  7:334438/149‑1 (MQ=255)
                                                            cGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑AGAGAATTGCCTTGatca                                                                                 >  6:367167/1‑149 (MQ=255)
                                                              gagTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGAAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑AGAGAATTGCCTTGatcatc                                                                               >  5:388555/3‑149 (MQ=255)
                                                              gCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑AGAGAATTGCCTTGatcatc                                                                               <  6:41618/149‑1 (MQ=255)
                                                                gTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑AGAGAATTGCCTTGa                                                                                    <  7:80265/142‑1 (MQ=255)
                                                                gTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑AGAGAATTGCCTTGa                                                                                    >  8:80265/1‑142 (MQ=255)
                                                                  ggc‑‑ggcGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑AGAGAATTGCCTTGATCATCTGgc                                                                           <  5:110741/149‑1 (MQ=255)
                                                                   gc‑‑ggcgCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑AGAGAATTGCCTTGATCATCTGgcg                                                                          >  5:102172/1‑149 (MQ=255)
                                                                         cgcCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑AGAGAATTGCCTTGATCATCTGGCGCATg                                                                      <  5:125828/149‑1 (MQ=255)
                                                                              aGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGg‑agaga                                                                                              >  7:141813/1‑120 (MQ=255)
                                                                              aGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGg‑agaga                                                                                              <  8:141813/120‑1 (MQ=255)
                                                                                 ttGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑AGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTc                                                              <  1:162346/149‑1 (MQ=255)
                                                                                     tCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGACAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑AGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAgc                                                          <  5:367167/149‑1 (MQ=255)
                                                                                              ccATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑AGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCggg                                                 <  5:383202/149‑1 (MQ=255)
                                                                                                   gCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑AGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTAt                                             <  2:153932/148‑1 (MQ=255)
                                                                                                     tGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑AGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCa                                          >  6:349846/1‑149 (MQ=255)
                                                                                                       cAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑AGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATg                                        <  4:234793/149‑1 (MQ=255)
                                                                                                        aaTATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑AGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGc                                       <  6:225524/149‑1 (MQ=255)
                                                                                                        aaTATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑AGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGc                                       <  1:3717/149‑1 (MQ=255)
                                                                                                           atGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑AGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGa                                    >  6:54528/1‑149 (MQ=255)
                                                                                                                  gcGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑AGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAaccac                             <  3:338839/149‑1 (MQ=255)
                                                                                                                       gagaATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑AGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGtt                        <  6:349046/149‑1 (MQ=255)
                                                                                                                       gagaATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑AGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGtt                        >  3:5879/1‑149 (MQ=255)
                                                                                                                             ataGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑AGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTgcgcgc                  >  6:13278/1‑149 (MQ=255)
                                                                                                                               aGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGGAAGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGcc                 >  1:302933/1‑149 (MQ=255)
                                                                                                                               aGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑AGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCt                >  2:344331/1‑149 (MQ=255)
                                                                                                                               aGGTGCCGGCGCTGGCGGCGACGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCAGG‑AGAGAATTGCCTTGATCATCTGGCGCATGGCACAGTCCAGCCCCCCCGGGGTATCCAAGCGGATAACCACCAATTGCGCGcc                 >  4:50555/1‑148 (MQ=255)
                                                                                                                                         gcTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑AGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCTTGCgcctggg      <  6:388555/149‑1 (MQ=255)
                                                                                                                                             gcggcgCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑GGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCTTGCgcctgggcctg  >  3:231509/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                              
ATCGTCGCCGCCTTTGGCTTTATCGTCACTGGGCGTGCCAGGCGCGCCGCCGATCCGCACCGGCGTGGC‑‑GGCGCCCAGGTTGGTCCCGGGCGCCATCGCTGCAATATGGCTGGCGTAGAGAATATAGGTGCCGGCGCTGGCGGCGCGGGCGCCGCTGGGGGCAACGTAGCTGGCGACGGGTACGGTGCTGG‑AGAGAATTGCCTTGATCATCTGGCGCATGGCACTGTCCAGCCCGCCCGGGGTATCCATGCGGATAACCACCAGTTGCGCGCCTTGCGCCTGGGCCTG  >  NC_002947/250765‑251052

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