Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,488,197 G→C 10.5% intergenic (+53/‑349) 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,1970GC10.5% 158.4 / 9.6 57intergenic (+53/‑349)PP_2182/PP_2183hypothetical protein/formate dehydrogenase subunit gamma
Reads supporting (aligned to +/- strand):  ref base G (24/27);  new base C (3/3);  total (27/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 = 7.68e-01

CGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCA  >  NC_002947/2488049‑2488332
                                                                                                                                                      |                                                                                                                                       
cGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCgggg                                                                                                                                         <  1:137069/149‑1 (MQ=255)
cGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCgggg                                                                                                                                         <  5:21090/149‑1 (MQ=255)
     ccAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGcc                                                                                                                                    >  8:226631/1‑149 (MQ=255)
         gACCGACGCAATCTGGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTTAGCGCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGccagac                                                                                                                                <  6:7476/149‑5 (MQ=255)
           ccGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTgg                                                                                                                              >  8:233981/1‑149 (MQ=255)
             gACGAAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGCCACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCCTTGCCGGGGCTGCCGTCTGGcc                                                                                                                            >  3:170225/1‑148 (MQ=255)
              aCGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAg                                                                                                                           >  3:98181/1‑149 (MQ=255)
                   aTCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCg                                                                                                                      <  1:225085/149‑1 (MQ=255)
                        gCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAAt                                                                                                                 <  3:180594/149‑1 (MQ=255)
                           aaGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGa                                                                                                              <  7:226631/149‑1 (MQ=255)
                                aaaGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATAGGCGTACCCAATTCCGCGGCTCGCGGCTGGCAGGCGCGGCAAAGGAGTa                                                                                                           >  1:272348/1‑149 (MQ=255)
                                aaaGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGTCCATTGCCGGGCCTGCCAGACGGCAGCCCCGGc                                                                                                                    >  4:263865/1‑138 (MQ=38)
                                aaaGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGTCCATTGCCGGGCCTGCCAGACGGCAGCCCCGGc                                                                                                                    <  3:263865/138‑1 (MQ=38)
                                aaaGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAgg                                                                                                         >  7:33330/1‑149 (MQ=255)
                                aaaGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCACGTAGACATGCTGGATGGCAGCCCATTGCCGGGGCTCCCGTCTGGCAGGCCCGGCAATGAAATAgg                                                                                                         >  8:299393/1‑149 (MQ=255)
                                      ccaggcccaggcTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCggg                                                                                                   >  3:262636/1‑148 (MQ=255)
                                      ccaggcccaggcTCCGCCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCATTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTCCCGCGGCTGCCGTCTGGCAGGCCCGGCAATGCACTAGCTTCGGt                                                                                                   >  4:244635/1‑149 (MQ=255)
                                           cccaggcTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCga                                                                                              <  8:33330/149‑1 (MQ=255)
                                            ccaggcTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCgag                                                                                             >  7:296806/1‑149 (MQ=255)
                                                        ccccGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGtgcat                                                                                 >  2:223731/1‑149 (MQ=255)
                                                           cgccgccgACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGtgcatgc                                                                               >  3:172009/1‑148 (MQ=255)
                                                           cgccgccgACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGtgcatgc                                                                               >  4:77432/1‑148 (MQ=255)
                                                            gccgccgACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGCCTGCCAGACGGCAGCCCCGGCAATgg                                                                                                               <  6:177416/115‑1 (MQ=37)
                                                            gccgccgACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGCCTGCCAGACGGCAGCCCCGGCAATgg                                                                                                               >  5:177416/1‑115 (MQ=37)
                                                                 cgACAACGCAACGGGGCATTCTTCCGACGGCGCGCAGTCGCGTGCTGCCCCGTATCCATGCTGGATGGCAGCCCCTTGCCGGGGCTGCCGTCTTGCAGGCCCTGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATg                                                                        <  4:170225/149‑1 (MQ=255)
                                                                  gACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATgg                                                                       >  3:144429/1‑149 (MQ=255)
                                                                      ccGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTc                                                                   <  7:233981/149‑1 (MQ=255)
                                                                        gCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTc                                                                   >  3:103044/1‑147 (MQ=255)
                                                                        gCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCaa                                                                 <  1:223731/149‑1 (MQ=255)
                                                                          aaCGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTGTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGa                                                               <  5:198130/149‑1 (MQ=255)
                                                                          aaCGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGtgcatg                                                                                >  1:277747/1‑132 (MQ=255)
                                                                          aaCGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGtgcatg                                                                                <  2:277747/132‑1 (MQ=255)
                                                                                  aGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCt                                                        >  5:316130/1‑148 (MQ=255)
                                                                                    tCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGAt                                                     <  4:172009/149‑1 (MQ=255)
                                                                                        ccGACTGAGCTCAGTCGCGTTCTGCCCCGTATCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCgg                                                  <  3:244635/148‑1 (MQ=255)
                                                                                          gACTGAGCTCAGTCGCGTTCTGCCCGGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCaa                                               <  7:299393/149‑1 (MQ=255)
                                                                                               aGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGt                                          >  8:281002/1‑149 (MQ=255)
                                                                                                gCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTg                                         >  6:221644/1‑149 (MQ=255)
                                                                                                gCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAAACCCACCTGATCGGCAAGCCGTg                                         >  2:313691/1‑149 (MQ=255)
                                                                                                  tCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGc                                       <  7:281002/149‑1 (MQ=255)
                                                                                                      tCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAgat                                   >  4:281797/1‑149 (MQ=255)
                                                                                                              tGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCgg                           >  2:26895/1‑149 (MQ=255)
                                                                                                                ccccGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGTCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCgggg                         <  6:11285/149‑1 (MQ=255)
                                                                                                                 cccGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGc                        <  5:181863/149‑1 (MQ=255)
                                                                                                                 cccGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGc                        <  4:103044/149‑1 (MQ=255)
                                                                                                                 cccGCAGCCCTGCTGGACGCCAGCCCAGTGCCGGGGCGGCCGTCTGGCAGGCCCGGCAATGGACTAGGTGCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGc                        <  4:262636/149‑1 (MQ=255)
                                                                                                                  ccGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCg                       <  4:144429/149‑1 (MQ=255)
                                                                                                                    gTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCTACGATGGGTTCAAGACCCACCTGATCGGCAAGCCGTGGGAGATGATGCCGGGGCGt                      >  3:44861/1‑148 (MQ=255)
                                                                                                                      aGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCTGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTaa                    <  5:221644/148‑1 (MQ=255)
                                                                                                                       gCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATc                  <  3:77432/149‑1 (MQ=255)
                                                                                                                          aTGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACg               >  3:44939/1‑149 (MQ=255)
                                                                                                                           tGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGt              <  1:26895/149‑1 (MQ=255)
                                                                                                                           tGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGt              <  2:272348/149‑1 (MQ=255)
                                                                                                                                gATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGTGTAATCACGTGCGt          <  5:161483/148‑1 (MQ=255)
                                                                                                                                gATGGCAGCCCATTGCCGGGGCGTCCGTCTGGGAGGCCCGGCAATGGACTAGGTGCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGt          <  8:296806/148‑1 (MQ=255)
                                                                                                                                     gCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGtt     <  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)
                                                                                                                                                      |                                                                                                                                       
CGCAGCCAGGACCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCA  >  NC_002947/2488049‑2488332

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: