Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,240,243 G→A 7.0% intergenic (+20/‑105) PP_2837 → / → PP_2838 uronate transporter/hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,240,2430GA 7.0% 159.3 / 2.9 58intergenic (+20/‑105)PP_2837/PP_2838uronate transporter/hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (24/29);  new base A (2/2);  total (26/32)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.11e-01

CAGCGTTCATGATCAGCGCGGTGGTGGCGCTGTGCTCGGCGGTGGCGTACTTCCTGCTGCTACAGGCACCCATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAGGGGTATCTGTATGTATAAACAGTGCCTGTTGCCGCAACCTGACCTCCGCCCGTGATCGCCCATTCCGTCGACGACCCGTTCTACTACCTGC  >  NC_002947/3240096‑3240387
                                                                                                                                                   |                                                                                                                                                
caGCGTTCATGATCAGCGAGGTGGTGGCGCGGTGCGCGGCGGTGGCGTACTTCCTGCTGCTACAGGCACCCATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCtgcagtgca                                                                                                                                                <  2:1374570/150‑1 (MQ=255)
  gCGTTCATGATCAGCGCGGTGGTGGCGCTGTGCTCGGCGGTGGCGTACTTCCTGCTGCTACAGGCACCCATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGc                                                                                                                                              <  1:1279724/150‑1 (MQ=255)
    gTTCATGATCAGCGCGGTGGTGGCGCTGTGCTCGGCGGTGGCGTACTTCCTGCTGCTACAGGCACCCATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGccc                                                                                                                                            >  1:1215747/1‑150 (MQ=255)
           aTCAGCGCGGTGGTGGCGCTGTGCTCGGCGGTGGCGTACTTCCTGCTGCTACAGGCACCCATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGcccccgtagc                                                                                                                                     <  2:73226/150‑7 (MQ=255)
                cgcgGTGGTGGCGCTGTGCTCGGCGGTGGCGTACTTCCTGCTGCTACAGGCACCCGTTGGTTGTGCCAATGTCGCAGCGGTCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGt                                                                                                                                <  1:1028760/150‑1 (MQ=255)
                  cgGTGGTGGCGCTGTGCTCGGCGGTGGCGTACTTCCTGCTGCTACAGGCACCCATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTg                                                                                                                              >  2:1206306/1‑150 (MQ=255)
                     tggtggCGCTGTGCTCGGCGGTGGCGTACTTCCTGCTGCTACAGGCACCCATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATc                                                                                                                           >  2:1190537/1‑150 (MQ=255)
                       gtggCGCTGTGCTCGGCGGTGGCGTACTTCCTGCTGCTACAGGCACCCATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGt                                                                                                                         <  1:1206306/150‑1 (MQ=255)
                                  cTCGGCGGTGGCGTACTTCCTGCTGCTACAGGCACCCATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTACAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCcg                                                                                                              >  1:544130/1‑150 (MQ=255)
                                  cTCGGCGGTGGCGTACTTCCTGCTGCTACAGGCACCCATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTACAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCcg                                                                                                              >  1:543989/1‑150 (MQ=255)
                                   tCGGCGGTGGCGTACTTCCTGCTGCTACAGGCACCCATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCcgc                                                                                                             >  1:653245/1‑150 (MQ=255)
                                   tCGGCGGTGGCGTACTTCCTGCTGCTACAGGCACCCATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCcgc                                                                                                             >  1:268755/1‑150 (MQ=255)
                                     ggcggTGGCGTACTTCCTGCTGCTACAGGCACCCATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCgcgg                                                                                                           <  2:268755/150‑1 (MQ=255)
                                        ggTGGCGTACTTCCTGCTGCTACACGCTCCCATTTGCAGTGCCACGTTCGCGGCGGGCTTGATGGCTGGCGCAACTGAACAGGCGTGATCGCTACGTTGGCTGCAGTGCTGCCCCTGAACCCTGGTTCATCGTGAGGCGACCCGCGgcgc                                                                                                        <  2:52972/150‑1 (MQ=255)
                                             cGTACTTCCTGCTGCTACAGGCACCCATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTc                                                                                                   >  2:1060606/1‑150 (MQ=255)
                                              gTACTTCCTGCTGCTACAGGCACCCATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCa                                                                                                  >  1:145284/1‑150 (MQ=255)
                                               tACTTCCTGCTGCTACAGGCACCCATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCAt                                                                                                 >  2:558310/1‑150 (MQ=255)
                                                aCTTCCTGCTGCTACAGGCACCCATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACAGGGGCTGCACTGCAGcccccgta                                                                                                                                       <  1:1131482/111‑5 (MQ=255)
                                                aCTTCCTGCTGCTACAGGCACCCATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACAGGGGCTGCACTGCAGcccccgta                                                                                                                                       >  2:1131482/1‑107 (MQ=255)
                                                  ttCCTGCTGCTACAGGCACCCATTGGCAGTGCCTAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCGTGGa                                                                                              >  1:990696/1‑150 (MQ=255)
                                                   tCCTGCTGCTACAGGCACCCATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAg                                                                                             <  1:888919/150‑1 (MQ=255)
                                                      tgctgcTACAGGCACCCATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTACAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAgggg                                                                                          <  2:543989/150‑1 (MQ=255)
                                                      tgctgcTACAGGCACCCATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTACAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAgggg                                                                                          <  2:544130/150‑1 (MQ=255)
                                                       gctgctACAGGCACCCATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAGGGGt                                                                                         <  2:1215747/150‑1 (MQ=255)
                                                           ctACAGGCACCCATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAGGGGTATCt                                                                                     <  1:544438/150‑1 (MQ=255)
                                                             aCAGGCACCCATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATg                                                                                                >  2:1006727/1‑137 (MQ=255)
                                                             aCAGGCACCCATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATg                                                                                                <  1:1006727/137‑1 (MQ=255)
                                                               aGGCACCCATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAGGGGTATCtgtat                                                                                 <  2:653245/150‑1 (MQ=255)
                                                                   aCCCATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAGGGGt                                                                                         >  2:827178/1‑138 (MQ=255)
                                                                   aCCCATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAGGGGt                                                                                         <  1:827178/138‑1 (MQ=255)
                                                                    cccATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTAAGGGGGCTGCAGTGCAGCCCCTGAACACTGGTTGATCGTGAGGCGAGCCGCGGAGCATGTCATGGAGGGGGATCTGTATGtata                                                                            >  2:567795/1‑150 (MQ=255)
                                                                      cATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATTGAGGGGTATCTGTATGTATaaa                                                                          >  2:685223/1‑150 (MQ=255)
                                                                      cATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAGGGGTATCTGTATGTATaaa                                                                          >  1:67038/1‑150 (MQ=255)
                                                                        ttGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAGGGGTATCTGTATGTATAAACa                                                                        <  1:1060606/150‑1 (MQ=255)
                                                                         tGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTCCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAGGGGTATCTGTATGTATAAACAg                                                                       <  1:648429/150‑1 (MQ=255)
                                                                            cAGTGCCAAGGTCGAGGGGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAGGGGTATCTGTATGTATAAACAGTGc                                                                    <  1:558310/150‑1 (MQ=255)
                                                                               tGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAGGGGTATCTGTATGTATAAACAGTGCCTg                                                                 >  2:452564/1‑150 (MQ=255)
                                                                               tGCCAAGGCCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGCGGCGACCCGCGGCGCTTGTCATGGAGGGCTATCTGTATGTATAAACAGTGCCTg                                                                 <  1:567795/150‑1 (MQ=255)
                                                                                gCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAGGGGTATCTGTATGTATAAACa                                                                        <  1:263798/142‑1 (MQ=255)
                                                                                gCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAGGGGTATCTGTATGTATAAACa                                                                        >  2:263798/1‑142 (MQ=255)
                                                                                gCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAGGGGTATCTGTATGTATAAACAGTGCCTGt                                                                <  1:577712/150‑1 (MQ=255)
                                                                                                          tggcGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATgg                                                                                               >  1:739670/1‑93 (MQ=255)
                                                                                                          tggcGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATgg                                                                                               <  2:739670/93‑1 (MQ=255)
                                                                                                              gcAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAGGGGTATCTGTATGTATAAACAGTGCCTGTTGCCGCAACCTGACCTCCGCCCGTGATCGc                                  <  1:452564/150‑1 (MQ=255)
                                                                                                              gcAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAGGGGTATCTGTATGTATAAACAGTGCCTGTTGCCGCAACCTGACCTCCGCCCGTGATCGc                                  >  1:141319/1‑150 (MQ=255)
                                                                                                                aaCTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAGGGGTATCTGTATGTATAAACAGTGCCTGTTGCCGCAACCTGACCTCCGCCCGTGATCGccc                                <  1:556042/150‑1 (MQ=255)
                                                                                                                 aCTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAg                                                                                             <  2:861068/88‑1 (MQ=255)
                                                                                                                 aCTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAg                                                                                             >  1:861068/1‑88 (MQ=255)
                                                                                                                    gAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTTATGCGACCCGCGGCGCTTGTCATGGAGGGGTATCTGTATGTATAAACAGTGCCTGTTGCCGCAGCCTGACCTCCGCCCGTGATCGGCCATTc                            <  2:145284/150‑1 (MQ=255)
                                                                                                                     aaCAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAGGGGTATCTGTATGTATaaa                                                                          <  2:812462/103‑1 (MQ=255)
                                                                                                                     aaCAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAGGGGTATCTGTATGTATaaa                                                                          >  1:812462/1‑103 (MQ=255)
                                                                                                                      aCAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAGGGGTATCTGTATGTATAAACAGTGCCTGTTGCCGCAACCTGACCTCCGCCCGTGATCGCCCATTCCg                          <  1:280332/150‑1 (MQ=255)
                                                                                                                      aCAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAGGGGTATCTGTATGTATAAACAGTGCCTGTTGCCGCAACCTGACCTCCGCCCGTGATCGCCCATTCCg                          <  2:141319/150‑1 (MQ=255)
                                                                                                                          gCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAGTGGTATCTGTATGTATAAACAGTGCCTGTTGCCGCAACCTGACCTCCGCCCGTGATCGCCCATTCCGTcga                      <  1:1190537/150‑1 (MQ=255)
                                                                                                                             tGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAGGGGTATCTGTATGTATAAACAGTGCCTGTTGCCGCAACCTGACCTCCGCCCGTGATCGCCCATTCCGTcgacga                   >  2:78610/1‑150 (MQ=255)
                                                                                                                              gAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAGGGGTATCTGTATGTATAAACAGTGCCTGTTGCCGCAACCTGACCTCCGCCCGTGATCGCCCATTCCGTCgacgac                  <  1:357396/150‑1 (MQ=255)
                                                                                                                                 cGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAGGGGTATCTGTATGTATaaa                                                                          >  1:1213774/1‑91 (MQ=255)
                                                                                                                                 cGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAGGGGTATCTGTATGTATaaa                                                                          <  2:1213774/91‑1 (MQ=255)
                                                                                                                                            cTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAGGGGTATCTGTATGTATAAACAGTGCCTGTTGCCGCAACCTGACCTCCGCCCGTGATCGCCCATTCCGTCGACGACCCGTTCTACTACCt    >  1:1083147/1‑150 (MQ=255)
                                                                                                                                             tgcagtgcagCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCGTGGAGGGGTATCTGTATGTATAAACAGTGCCTGTTGCCGCAACCTGACCTCCGCCCGTGATCGCCCATTCCGTCGACGACCCGTCCTACTACCTg   <  2:990696/150‑1 (MQ=255)
                                                                                                                                              gcagtgcagCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAGGGGTATCTGTATGTATAAACAGTGCCTGTTGCCGCAACCTGACCTCCGCCCGTGATCGCCCATTCCGTCGACGACCCGTTCTACTACCTGc  >  1:1405303/1‑150 (MQ=255)
                                                                                                                                               cagtgcagCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAGGGGTATCTGTATGTATAAACAGTGCCTGTTGCCGCAACCTGACCTCCGCCCGTGATCGCCCATTCCGTCgacgac                  >  2:540396/1‑133 (MQ=255)
                                                                                                                                               cagtgcagCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAGGGGTATCTGTATGTATAAACAGTGCCTGTTGCCGCAACCTGACCTCCGCCCGTGATCGCCCATTCCGTCgacgac                  <  1:540396/133‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                                
CAGCGTTCATGATCAGCGCGGTGGTGGCGCTGTGCTCGGCGGTGGCGTACTTCCTGCTGCTACAGGCACCCATTGGCAGTGCCAAGGTCGAGGCGGGCGGGATGGCTGGCGCAACTGAACAGGCGTGAACGCTACGGGGGCTGCAGTGCAGCCCCTGAACCCTGGTTGATCGTGAGGCGACCCGCGGCGCTTGTCATGGAGGGGTATCTGTATGTATAAACAGTGCCTGTTGCCGCAACCTGACCTCCGCCCGTGATCGCCCATTCCGTCGACGACCCGTTCTACTACCTGC  >  NC_002947/3240096‑3240387

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 12 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 37 ≤ ATCG/ATCG < 41 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑