Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,488,212 G→C 9.7% intergenic (+68/‑334) 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,2120GC 9.7% 177.2 / 5.9 60intergenic (+68/‑334)PP_2182/PP_2183hypothetical protein/formate dehydrogenase subunit gamma
Reads supporting (aligned to +/- strand):  ref base G (27/27);  new base C (3/3);  total (30/30)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

ATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCAGAATGCGCCTGGTCGCTCGCATACAACC  >  NC_002947/2488068‑2488360
                                                                                                                                                  |                                                                                                                                                    
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                                                                                                            >  1:277747/1‑132 (MQ=255)
                                                       aaCGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGtgcatg                                                                                                            <  2:277747/132‑1 (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                                                    <  4:103044/149‑1 (MQ=255)
                                                                                              cccGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGc                                                    <  5:181863/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                                          <  1:26895/149‑1 (MQ=255)
                                                                                                        tGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGt                                          <  2:272348/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                              <  4:203830/144‑1 (MQ=39)
                                                                                                                           ttGCCGGGCCTGCCAGACGGCAGCCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCa                              >  3:203830/1‑144 (MQ=39)
                                                                                                                                    cTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCAGAATGCGCCTGGtc                >  3:194671/1‑149 (MQ=255)
                                                                                                                                     tGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCAGAATGCGCCTGGtcg               >  3:281906/1‑149 (MQ=255)
                                                                                                                                     tGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCAGAATGCGCCTGGtcg               >  1:5763/1‑149 (MQ=255)
                                                                                                                                       ccGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCAGAATGCGCCTGGtcgct             >  3:189791/1‑149 (MQ=255)
                                                                                                                                       ccGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCAGAATGCGCCTGGtcgct             >  4:195622/1‑149 (MQ=255)
                                                                                                                                          tCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCAGAATGCGCCTGGtcgctcg           >  6:280514/1‑148 (MQ=255)
                                                                                                                                              gCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGCCCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCAGAATGCGCCTGGTCGCTCGCATa       >  6:210988/1‑148 (MQ=255)
                                                                                                                                               cAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAAACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTACTCACGTGCGTCGGTTTCAGAATGCGCCTGGTCGCTCGCATACa     <  1:313691/149‑1 (MQ=255)
                                                                                                                                                  gCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGt                                          >  7:196109/1‑109 (MQ=255)
                                                                                                                                                  gCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGt                                          <  8:196109/109‑1 (MQ=255)
                                                                                                                                                  gCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCAGAATGCGCCTGGTCGCTCGCATACAAcc  >  3:285295/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGAT‑GGC‑AGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCAGAATGCGCCTGGTCGCTCGCATACAACC  >  NC_002947/2488068‑2488360

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: ‑