Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,763,547 G→C 7.9% intergenic (+24/+1651) PP_3326 → / ← PP_3328 hypothetical protein/ring‑cleaving dioxygenase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,763,5470GC 7.9% 193.2 / 7.7 63intergenic (+24/+1651)PP_3326/PP_3328hypothetical protein/ring‑cleaving dioxygenase
Reads supporting (aligned to +/- strand):  ref base G (26/32);  new base C (3/2);  total (29/34)
Fisher's exact test for biased strand distribution p-value = 6.54e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.18e-01

CGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCATGATGGGCACTTGCAT  >  NC_002947/3763401‑3763694
                                                                                                                                                  |                                                                                                                                                   
cGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACg                                                                                                                                                  <  1:208032/150‑1 (MQ=255)
    aTGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCg                                                                                                                                              <  1:16521/150‑1 (MQ=255)
      gAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCGTCATGGCCGAg                                                                                                                                            <  2:1058816/150‑1 (MQ=255)
       aaGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCGTCATGGCCGAgg                                                                                                                                           <  2:476825/150‑1 (MQ=255)
        aGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGt                                                                                                                                          >  2:1000787/1‑150 (MQ=255)
               aaCATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGa                                                                                                                                                    <  2:730561/133‑1 (MQ=255)
               aaCATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGa                                                                                                                                                    >  1:730561/1‑133 (MQ=255)
                aCATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGt                                                                                                                                  >  1:593190/1‑150 (MQ=255)
                 catcatCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGtt                                                                                                                                 <  1:1141043/150‑1 (MQ=255)
                  atcatcGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCATCATAGCTTGTGCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTg                                                                                                                                <  2:495457/150‑1 (MQ=255)
                    catcGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTTTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGa                                                                                                                              >  2:941768/1‑150 (MQ=255)
                     atcGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAg                                                                                                                             <  1:570059/150‑1 (MQ=255)
                           gCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTc                                                                                                                       <  1:455105/150‑1 (MQ=255)
                            cccAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCg                                                                                                                      <  1:170735/150‑1 (MQ=255)
                                    aCCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGt                                                                                                              <  2:99481/150‑1 (MQ=255)
                                      ccTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGt                                                                                                              <  1:1076642/148‑1 (MQ=255)
                                      ccTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGt                                                                                                              >  2:1076642/1‑148 (MQ=255)
                                      ccTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTCAGTATGCCATACCGCTAGTAAACTAAGCCACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGc                                                                                                            <  1:797043/150‑1 (MQ=255)
                                        tGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCt                                                                                                          <  2:257400/150‑1 (MQ=255)
                                        tGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCt                                                                                                          <  2:823517/150‑1 (MQ=255)
                                        tGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCt                                                                                                          >  1:823517/1‑150 (MQ=255)
                                        tGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGACGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCt                                                                                                          >  1:269937/1‑150 (MQ=255)
                                               ggCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCa                                                                                                   <  1:376657/150‑1 (MQ=255)
                                                     cGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATcc                                                                                             <  1:178998/150‑1 (MQ=255)
                                                               gTGGCCCAGCAGCGGCTGTTCCGTTCATCGTCACCCCGCGGTAAGCCGTTTCACTAGTAGACTATGCTACATACCTGGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGAc                                                                                   <  2:141269/150‑1 (MQ=255)
                                                               gTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGAc                                                                                   <  2:503336/150‑1 (MQ=255)
                                                                    ccAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTg                                                                              >  2:967386/1‑150 (MQ=255)
                                                                     cagcagAGCCTGTTAAGTTCCCTGTCGTTCCTGGGTAAGACGCACCACGAGTAAACTACGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTgg                                                                             <  2:269937/150‑1 (MQ=255)
                                                                          gagCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACaa                                                                        >  1:454665/1‑150 (MQ=255)
                                                                          gagCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACaa                                                                        >  2:567995/1‑150 (MQ=255)
                                                                               tGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCg                                                                   >  1:49536/1‑150 (MQ=255)
                                                                                     cttcaCTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCg                                                             <  1:600201/149‑1 (MQ=255)
                                                                                        caCTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGACTGCGTGCCTGCAACCAGCATCCACACCATCACCATTTGCACAAGCCCGGCTGCggag                                                          >  2:600201/1‑148 (MQ=255)
                                                                                         aCTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCca                                                         >  2:832732/1‑150 (MQ=255)
                                                                                          cTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCcac                                                        >  2:164496/1‑150 (MQ=255)
                                                                                                     gggTAAGCCGCACAAAGAGTAAACTAAGCGACATAGCTCGGCCATGACGGCCGAGGTATGACAGTTGCAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGa                                             >  2:308551/1‑150 (MQ=255)
                                                                                                      ggTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCAt                                                                                               <  1:868921/99‑1 (MQ=255)
                                                                                                      ggTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCAt                                                                                               >  2:868921/1‑99 (MQ=255)
                                                                                                             cGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATg                                     <  2:199772/150‑1 (MQ=255)
                                                                                                              gCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATgg                                    >  2:663912/1‑150 (MQ=255)
                                                                                                                     gAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCAtccc                             <  1:567995/150‑4 (MQ=255)
                                                                                                                     gAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTc                             >  1:561075/1‑150 (MQ=255)
                                                                                                                     gAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTc                             <  1:967386/150‑1 (MQ=255)
                                                                                                                          aaCTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAac                                           >  1:930874/1‑131 (MQ=255)
                                                                                                                          aaCTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAac                                           <  2:930874/131‑1 (MQ=255)
                                                                                                                             tAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACCCGGCAATCTAAGGCGATGGGCATGTCGATTCCCt                     <  2:49536/150‑1 (MQ=255)
                                                                                                                             tAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCt                     >  1:562029/1‑150 (MQ=255)
                                                                                                                             caactgccATACCTCGGCCGTCATGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCt                     >  2:519378/8‑150 (MQ=39)
                                                                                                                                gCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCa                  <  2:562029/150‑1 (MQ=255)
                                                                                                                                 cGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCat                 >  2:176279/1‑150 (MQ=255)
                                                                                                                                  gACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCatg                >  2:641755/1‑150 (MQ=255)
                                                                                                                                   aCATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCatga               >  1:1279799/1‑150 (MQ=255)
                                                                                                                                   aCATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGAGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCAtt                                                                               >  1:107782/1‑86 (MQ=255)
                                                                                                                                   aCATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGAGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCAtt                                                                               <  2:107782/86‑1 (MQ=255)
                                                                                                                                       aCCTCGGCCGTCATGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTAACTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCATGATggg           >  2:816107/1‑150 (MQ=255)
                                                                                                                                       aCCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGc                                                                            <  2:158422/85‑1 (MQ=255)
                                                                                                                                       aCCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGc                                                                            >  1:158422/1‑85 (MQ=255)
                                                                                                                                       aCCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCATGATggg           <  1:832732/150‑1 (MQ=255)
                                                                                                                                        ccTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCATGATGGGc          <  1:941768/150‑1 (MQ=255)
                                                                                                                                         cTCGGCCGTCATGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCATGATGGGCa         >  1:480684/1‑150 (MQ=255)
                                                                                                                                             gCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCa                                                                                       <  2:367644/68‑1 (MQ=255)
                                                                                                                                             gCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCa                                                                                       >  1:367644/1‑68 (MQ=255)
                                                                                                                                                aTGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCAGTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCATGATGGGCACTTGCAt  <  1:663912/150‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                   
CGCAATGAAGGACTCAACATCATCGCGGCCCAGGCTACCCTGCCAATGGCAGACGGCAGCCCCGTGGCCCAGCAGAGCCTGTTCAGTTCACTGTCGTTCCTGGGTAAGCCGCACCACGAGTAAACTAAGCGACATACCTCGGCCATGACGGCCGAGGTATGGCAGTTGGAGGCGTTCGGAATGAGTGCCTGCAACCAGCATCCGCACCATGACCATTGGCACAAGCCCGGCTGCGGCCACACGGCAATCGAACACGATGGGCATGTCGATTACCTGCATGATGGGCACTTGCAT  >  NC_002947/3763401‑3763694

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