Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,488,203 G→A 10.0% intergenic (+59/‑343) 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,2030GA10.0% 159.1 / 6.9 60intergenic (+59/‑343)PP_2182/PP_2183hypothetical protein/formate dehydrogenase subunit gamma
Reads supporting (aligned to +/- strand):  ref base G (28/26);  new base A (3/3);  total (31/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 = 1.00e+00

CCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCAGAATGCGCCTGGTCGCT  >  NC_002947/2488060‑2488349
                                                                                                                                                 |                                                                                                                                                  
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                                                                                                 <  2:277747/132‑1 (MQ=255)
                                                               aaCGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGtgcatg                                                                                                 >  1:277747/1‑132 (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                               <  2:272348/149‑1 (MQ=255)
                                                                                                                tGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGt                               <  1:26895/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)
                                                                                                                                            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)
                                                                                                                                                 |                                                                                                                                                  
CCGACGCAATCGTGCGAAGTCAAAGGCCCAGGCCCAGGCTCCGGCCCCCGCCGCCGACACCGCAACGGTGCAGTCTACCGACTGAGCTCAGTCGCGTTCTGCCCCGTAGCCATGCTGGATGGCAGCCCATTGCCGGGGCTGCCGTCTGGCAGGCCCGGCAATGGACTAGGTTCGGTTTCGAGAGCCTGGTGCATGCAACGATGGTTTCAAGACCCACCTGATCGGCAAGCCGTGGCAGATGATGCCGGGGCGTAATCACGTGCGTCGGTTTCAGAATGCGCCTGGTCGCT  >  NC_002947/2488060‑2488349

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: