Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 459,904 A→C 5.1% intergenic (‑50/+29) pqqB ← / ← pqqA coenzyme PQQ synthesis protein B/coenzyme PQQ synthesis protein A

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947459,9040AC 5.1% 254.4 / 4.4 79intergenic (‑50/+29)pqqB/pqqAcoenzyme PQQ synthesis protein B/coenzyme PQQ synthesis protein A
Reads supporting (aligned to +/- strand):  ref base A (41/34);  new base C (2/2);  total (43/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.51e-01

TCAGGGTGCCATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTACTACTTGG  >  NC_002947/459758‑460052
                                                                                                                                                  |                                                                                                                                                    
tCAGGGTGCCATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCa                                                                                                                                                      <  8:173951/147‑1 (MQ=255)
     gTGCCATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAg                                                                                                                                               <  7:41319/149‑1 (MQ=255)
          aTCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGt                                                                                                                                           >  3:47120/1‑148 (MQ=255)
          aTCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGt                                                                                                                                           >  4:200778/1‑148 (MQ=255)
           tCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCTCGGCCTGCCGGGGCGttt                                                                                                                                         >  4:326032/1‑147 (MQ=255)
            cACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGa                                                                                                                                         >  3:97032/1‑148 (MQ=255)
             aCGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTACGCAGAAACCGGATCGAAAAAACGAAAACGCCTCGGCCTGCCGGGGCGttttc                                                                                                                                       <  6:47716/149‑5 (MQ=255)
              cGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAAc                                                                                                                                      <  4:47120/149‑1 (MQ=255)
               ggTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCTCGGCCTGCCGGGGCGttttcgt                                                                                                                                     <  3:326032/149‑7 (MQ=255)
                    cccTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAg                                                                                                                                >  1:221830/1‑149 (MQ=255)
                      cTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGcc                                                                                                                              >  6:72285/1‑149 (MQ=255)
                          cAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTgg                                                                                                                           <  1:5651/148‑1 (MQ=255)
                                aCGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCg                                                                                                                    <  2:221830/149‑1 (MQ=255)
                                  gcagttgcagttCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCACAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAAGCCTAAGCCTGCCTTAGCGGt                                                                                                                  >  3:214216/1‑149 (MQ=255)
                                        gcagttCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACGTGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCg                                                                                                             <  4:72705/148‑1 (MQ=255)
                                          agttCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAg                                                                                                          >  5:214169/1‑149 (MQ=255)
                                           gttCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAg                                                                                                          >  2:41040/1‑148 (MQ=255)
                                              ccACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACa                                                                                                      >  7:232283/1‑149 (MQ=255)
                                               cACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACa                                                                                                      <  5:72285/148‑1 (MQ=255)
                                                aCTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATg                                                                                                    >  2:285563/1‑149 (MQ=255)
                                                aCTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATg                                                                                                    >  1:310366/1‑149 (MQ=255)
                                                    cGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGa                                                                                                <  4:285451/149‑1 (MQ=255)
                                                           ccGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGTAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTACGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCg                                                                                           >  5:226781/1‑147 (MQ=255)
                                                           ccGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCg                                                                                           >  7:278368/1‑147 (MQ=255)
                                                            cGCCACCGGCGGCGGATCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAg                                                                                        <  1:117307/149‑1 (MQ=255)
                                                               cACCGGCGGCGGAGCCGCGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCAAAGcc                                                                                      >  3:254417/1‑148 (MQ=255)
                                                                    gcggcggAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACg                                                                                >  4:171959/1‑149 (MQ=255)
                                                                    gcggcggAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACg                                                                                >  2:254357/1‑149 (MQ=255)
                                                                      ggcggAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCa                                                                              <  3:200778/149‑1 (MQ=255)
                                                                             ccGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAgt                                                                       <  3:171959/149‑1 (MQ=255)
                                                                              cGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAgtg                                                                      <  2:310366/149‑1 (MQ=255)
                                                                               gagaACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAgtgt                                                                     <  3:76001/149‑1 (MQ=255)
                                                                               gagaACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAgtg                                                                      >  6:138620/1‑148 (MQ=255)
                                                                                 gaACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTAt                                                                   <  5:41221/149‑1 (MQ=255)
                                                                                     cTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAg                                                               <  4:97032/149‑1 (MQ=255)
                                                                                     cTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAg                                                               >  8:214210/1‑149 (MQ=255)
                                                                                     cTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCAGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAg                                                               >  8:156734/1‑149 (MQ=255)
                                                                                        gATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGtt                                                            >  4:120686/1‑149 (MQ=255)
                                                                                        gATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGtt                                                            <  6:214169/149‑1 (MQ=255)
                                                                                          tGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTg                                                          <  5:138620/149‑1 (MQ=255)
                                                                                           gTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATa                                                                                  >  3:327723/1‑124 (MQ=255)
                                                                                           gTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATa                                                                                  <  4:327723/124‑1 (MQ=255)
                                                                                             aCATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTc                                                       >  2:78241/1‑149 (MQ=255)
                                                                                              cATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTc                                                       <  8:232283/148‑1 (MQ=255)
                                                                                              cATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTc                                                       <  1:78241/148‑1 (MQ=255)
                                                                                                 ggCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCcaca                                                   >  2:139517/1‑149 (MQ=255)
                                                                                                      ctctctATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGaa                                              >  3:32465/1‑149 (MQ=255)
                                                                                                       tctctATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAAt                                             >  7:97804/1‑149 (MQ=255)
                                                                                                          ctATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTa                                           >  4:56305/1‑148 (MQ=255)
                                                                                                           tATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTAc                                          >  1:105889/1‑148 (MQ=255)
                                                                                                           tATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACt                                         >  1:33111/1‑149 (MQ=255)
                                                                                                            aTGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTc                                        >  4:284505/1‑149 (MQ=255)
                                                                                                             tccAGAAACGGGATCGAAAAAACGATAACGCCCCGGCAGGCCGAGGCGTGAAACCCTTAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTcc                                       <  4:214216/147‑1 (MQ=255)
                                                                                                             tGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTc                                        <  1:254357/148‑1 (MQ=255)
                                                                                                                aGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACAt                                                                                                     <  1:66828/84‑1 (MQ=255)
                                                                                                                aGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACAt                                                                                                     >  2:66828/1‑84 (MQ=255)
                                                                                                                       ggATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGGGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAAt                               >  8:76966/1‑147 (MQ=255)
                                                                                                                       ggATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAAt                               >  3:136152/1‑147 (MQ=255)
                                                                                                                       ggATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGa                             <  8:97804/149‑1 (MQ=255)
                                                                                                                         aTCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgg                           <  1:41040/149‑1 (MQ=255)
                                                                                                                          tCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAggg                          >  6:223543/1‑149 (MQ=255)
                                                                                                                             aaaaaaCGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggttg                       >  8:262123/1‑149 (MQ=255)
                                                                                                                              aaaaaCGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggttgg                      >  5:202885/1‑149 (MQ=255)
                                                                                                                              aaaaaCGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggttgg                      >  3:234862/1‑149 (MQ=255)
                                                                                                                               aaaacgaaaatgCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggttggg                     >  3:276329/1‑149 (MQ=255)
                                                                                                                                 aacgaaaacgCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggttgggtt                   <  3:56305/149‑1 (MQ=255)
                                                                                                                                  acgaaaacgCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTAt                                                                   >  6:23199/1‑100 (MQ=255)
                                                                                                                                  acgaaaacgCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGggttgggttg                  <  4:234862/149‑1 (MQ=255)
                                                                                                                                    gaaaacgCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGgttgtt                >  5:312033/1‑149 (MQ=255)
                                                                                                                                    gaaaacgCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGgttgtt                <  3:120686/149‑1 (MQ=255)
                                                                                                                                     aaaacgCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTc               <  7:214210/149‑1 (MQ=255)
                                                                                                                                        acgCCTCGGCCTGCCGGGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGc            >  6:79731/1‑149 (MQ=255)
                                                                                                                                        acgCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGc            <  1:285563/149‑1 (MQ=255)
                                                                                                                                        acgCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGc            <  2:33111/149‑1 (MQ=255)
                                                                                                                                           ccccGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGctac         <  4:32465/149‑1 (MQ=255)
                                                                                                                                             ccGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggtt                        <  1:243087/132‑1 (MQ=255)
                                                                                                                                             ccGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAgggtt                        >  2:243087/1‑132 (MQ=255)
                                                                                                                                             ccGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGctacta       >  2:295163/1‑149 (MQ=255)
                                                                                                                                                 cAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTACTACTTg   <  6:312033/149‑1 (MQ=255)
                                                                                                                                                 cAGGCAGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTACTACTTg   <  7:156734/149‑1 (MQ=255)
                                                                                                                                                  aGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTACTACTTgg  <  4:136152/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
TCAGGGTGCCATCACGGTAGCCCTTGCAGTTGACGCAGTTGCAGTTCCACTGCGGGAACCCGCCACCGGCGGCGGAGCCGAGAACCTGGATGTACATGGCCACTCTCTATGCAGAAACCGGATCGAAAAAACGAAAACGCCCCGGCAGGCCGAGGCGTGAAACCCTAAGCCTGGCTTAGCGGTTAGCGAAGTACATGGTGACTTCGAAGCCGATACGCAGGTCAGTGTATGCAGGTTTGGTCCACATGGAATTACTCCTTCCGAATGAGGGTTGGGTTGTTCAGCTACTACTTGG  >  NC_002947/459758‑460052

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: