Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,488,205 C→A 9.4% intergenic (+61/‑341) PP_2182 → / → PP_2183 hypothetical protein/formate dehydrogenase subunit gamma

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,488,2050CA 9.4% 176.3 / 6.7 61intergenic (+61/‑341)PP_2182/PP_2183hypothetical protein/formate dehydrogenase subunit gamma
Reads supporting (aligned to +/- strand):  ref base C (29/26);  new base A (3/3);  total (32/29)
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.49e-01

CCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCAGAATGCGCCTGGTCGCTCG  >  NC_002947/2488060‑2488351
                                                                                                                                                   |                                                                                                                                                  
ccGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTgg                                                                                                                                                 >  8:233981/1‑149 (MQ=255)
  gACGAAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGCCACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCCTTGCCGGGGCTGCCGTCTGGcc                                                                                                                                               >  3:170225/1‑148 (MQ=255)
   aCGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAg                                                                                                                                              >  3:98181/1‑149 (MQ=255)
        aTCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCg                                                                                                                                         <  1:225085/149‑1 (MQ=255)
             gCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAAt                                                                                                                                    <  3:180594/149‑1 (MQ=255)
                aaGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGa                                                                                                                                 <  7:226631/149‑1 (MQ=255)
                     aaaGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATAGGCGTACCCAATTCCGCGGCTCGCGGCTGGCAGGCGCGGCAAAGGAGTa                                                                                                                              >  1:272348/1‑149 (MQ=255)
                     aaaGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGTCCATTGCCGGGCCTGCCAGACGGCAGCCCCGGc                                                                                                                                       >  4:263865/1‑138 (MQ=38)
                     aaaGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGTCCATTGCCGGGCCTGCCAGACGGCAGCCCCGGc                                                                                                                                       <  3:263865/138‑1 (MQ=38)
                     aaaGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAgg                                                                                                                            >  7:33330/1‑149 (MQ=255)
                     aaaGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCACGTAGACATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTCCCGTCTGGCAGGCCCGGCAATGAAATAgg                                                                                                                            >  8:299393/1‑149 (MQ=255)
                           ccaggcccaggcTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCggg                                                                                                                      >  3:262636/1‑148 (MQ=255)
                           ccaggcccaggcTCCGCCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCATTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTCCCGCGGCTGCCGTCTGGCAGGCCCGGCAATGCACTAGCTTCGGt                                                                                                                      >  4:244635/1‑149 (MQ=255)
                                cccaggcTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCga                                                                                                                 <  8:33330/149‑1 (MQ=255)
                                 ccaggcTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCgag                                                                                                                >  7:296806/1‑149 (MQ=255)
                                             ccccGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGtgcat                                                                                                    >  2:223731/1‑149 (MQ=255)
                                                cgccgccgACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGtgcatgc                                                                                                  >  3:172009/1‑148 (MQ=255)
                                                cgccgccgACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGtgcatgc                                                                                                  >  4:77432/1‑148 (MQ=255)
                                                 gccgccgACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGCCTGCCAGACGGCAGCCCCGGCAATgg                                                                                                                                  <  6:177416/115‑1 (MQ=37)
                                                 gccgccgACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGCCTGCCAGACGGCAGCCCCGGCAATgg                                                                                                                                  >  5:177416/1‑115 (MQ=37)
                                                      cgACAACGCAACGGGGCATTCTTCCGACGGCGCGCAGTCGCGTGCTGCCCCGTATCCATGCTGGAT‑GGC‑AGCCCCTTGCCGGGGCTGCCGTCTTGCAGGCCCTGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATg                                                                                           <  4:170225/149‑1 (MQ=255)
                                                       gACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATgg                                                                                          >  3:144429/1‑149 (MQ=255)
                                                           ccGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTc                                                                                      <  7:233981/149‑1 (MQ=255)
                                                             gCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTc                                                                                      >  3:103044/1‑147 (MQ=255)
                                                             gCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCaa                                                                                    <  1:223731/149‑1 (MQ=255)
                                                               aaCGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTGTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGa                                                                                  <  5:198130/149‑1 (MQ=255)
                                                               aaCGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGtgcatg                                                                                                   <  2:277747/132‑1 (MQ=255)
                                                               aaCGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGtgcatg                                                                                                   >  1:277747/1‑132 (MQ=255)
                                                                       aGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCt                                                                           >  5:316130/1‑148 (MQ=255)
                                                                         tCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGAt                                                                        <  4:172009/149‑1 (MQ=255)
                                                                             ccGACTGAGCTCAGTCGCGTTCTGCCCCGTATCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCgg                                                                     <  3:244635/148‑1 (MQ=255)
                                                                               gACTGAGCTCAGTCGCGTTCTGCCCGGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCaa                                                                  <  7:299393/149‑1 (MQ=255)
                                                                                    aGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGt                                                             >  8:281002/1‑149 (MQ=255)
                                                                                     gCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTg                                                            >  6:221644/1‑149 (MQ=255)
                                                                                     gCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAAACCCACCTGATCGGCAAGCCGTg                                                            >  2:313691/1‑149 (MQ=255)
                                                                                       tCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGc                                                          <  7:281002/149‑1 (MQ=255)
                                                                                           tCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAgat                                                      >  4:281797/1‑149 (MQ=255)
                                                                                                   tGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCgg                                              >  2:26895/1‑149 (MQ=255)
                                                                                                     ccccGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGTCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCgggg                                            <  6:11285/149‑1 (MQ=255)
                                                                                                      cccGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGc                                           <  5:181863/149‑1 (MQ=255)
                                                                                                      cccGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGc                                           <  4:103044/149‑1 (MQ=255)
                                                                                                      cccGCAGCCCTGCTGGAC‑GCC‑AGCCCAGTGCCGGGGCGGCCGTCTGGCAGGCCCGGCAATGGACTAGGTGCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGc                                           <  4:262636/149‑1 (MQ=255)
                                                                                                       ccGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCg                                          <  4:144429/149‑1 (MQ=255)
                                                                                                         gTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCTACGATGGGTTCAAGACCCACCTGATCGGCAAGCCGTGGGAGATGATGCCGGGGCGt                                         >  3:44861/1‑148 (MQ=255)
                                                                                                           aGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCTGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTaa                                       <  5:221644/148‑1 (MQ=255)
                                                                                                            gCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATc                                     <  3:77432/149‑1 (MQ=255)
                                                                                                               aTGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACg                                  >  3:44939/1‑149 (MQ=255)
                                                                                                                tGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGt                                 <  2:272348/149‑1 (MQ=255)
                                                                                                                tGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGt                                 <  1:26895/149‑1 (MQ=255)
                                                                                                                     gAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGTGTAATCACGTGCGt                             <  5:161483/148‑1 (MQ=255)
                                                                                                                     gAT‑GGC‑AGCCCATTGCCGGGGCGTCCGTCTGGGAGGCCCGGCAATGGACTAGGTGCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGt                             <  8:296806/148‑1 (MQ=255)
                                                                                                                          gC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGtt                        <  3:281797/149‑1 (MQ=255)
                                                                                                                             aGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTc                      <  6:316130/149‑1 (MQ=255)
                                                                                                                                   ttGCCGGGCCTGCCAGACGGCAGCCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCa                     >  3:203830/1‑144 (MQ=39)
                                                                                                                                   ttGCCGGGCCTGCCAGACGGCAGCCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCa                     <  4:203830/144‑1 (MQ=39)
                                                                                                                                            cTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCAGAATGCGCCTGGtc       >  3:194671/1‑149 (MQ=255)
                                                                                                                                             tGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCAGAATGCGCCTGGtcg      >  3:281906/1‑149 (MQ=255)
                                                                                                                                             tGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCAGAATGCGCCTGGtcg      >  1:5763/1‑149 (MQ=255)
                                                                                                                                               ccGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCAGAATGCGCCTGGtcgct    >  4:195622/1‑149 (MQ=255)
                                                                                                                                               ccGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCAGAATGCGCCTGGtcgct    >  3:189791/1‑149 (MQ=255)
                                                                                                                                                  tCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCAGAATGCGCCTGGtcgctcg  >  6:280514/1‑148 (MQ=255)
                                                                                                                                                   |                                                                                                                                                  
CCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCAGAATGCGCCTGGTCGCTCG  >  NC_002947/2488060‑2488351

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