Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,092,961 Δ1 bp 8.7% coding (709/981 nt) PP_4482 ← AraC family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,092,9610A. 8.7% 193.5 / 14.2 69coding (709/981 nt)PP_4482AraC family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (35/28);  new base . (3/3);  total (38/31)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.63e-01

GAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGTCTGCATCAGCATCTGCCGCGCCTTGTTCAACCGCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGAGTCAGCTTGGGGTGACTCGAGCCCAGGCGGTTCTGCAGCGGAATACGCTGGCGCTCGCCACCCTCGCGGATGC  >  NC_002947/5092813‑5093110
                                                                                                                                                    |                                                                                                                                                     
gagaAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGTCTGCATCAGCATCTGCCGCGCCTTGTTCAACCGCAACTCAAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGac                                                                                                                                                      <  1:1920575/150‑1 (MQ=255)
  gaAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGTCTGCATCAGCATCTGCCGCGCCTTGTTCAACCGCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGacac                                                                                                                                                    >  1:851246/1‑150 (MQ=255)
     gcaacaGCAGAGCCCGGTCTGGATGATCGACTTGCTGGTCTGCATCAGCATCTGCCGCGCCTTGTTCAACCGCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGacacaca                                                                                                                                                 <  2:1543945/147‑1 (MQ=255)
      ccACAGGACAGCCCGATCTGGATGATCGACTTGCTGGTCTGCATCAGCATCTGCCGCGCCTTGTTCAACCGCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGacacacac                                                                                                                                                >  2:721648/1‑150 (MQ=255)
           ggACAGCCCGATCTGGATGATCGACTTGCTGGTCTGCATCAGCATCTGCCGCGCCTTGTTCAACCGCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATg                                                                                                                                           <  1:46369/150‑1 (MQ=255)
             aCAGCCCGATCTGGATGATCGACTTGCTGGTCGGCAATAGCATCTGCCGCGCCGGGTTCAACCGCAACAGCAGGTAGTAATGGTCGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCCACATACACACAGGCt                                                                                                                                         >  1:1572111/1‑150 (MQ=255)
                   cGATCTGGATGATCGACTTGCTGGTCTGCATCAGCATCTGCCGCGCCTTGTTCAACCGCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGa                                                                                                                                   <  2:194656/150‑1 (MQ=255)
                    gATCTGGATGATCGACTTGCTGGTCTGCATCAGCATCTGCCGCGCCTTGTTCAACCGCAACTCCAGGTAGTACTGGCTGGGCACCCGGTTCAGGTACTGCTTGAAGATACGCTCCTACTGCCGAGGCGACACACACACGTGCTGGGCGAt                                                                                                                                  <  2:639248/150‑1 (MQ=255)
                    gATCTGGATGATCGACTTGCTGGTCTGCATCAGCATCTGCCGCGCCTTGTTCAACCGCAACTCCAGGTAGTACTGACTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGAt                                                                                                                                  >  2:309437/1‑150 (MQ=255)
                          gttgttCGACTTGCTGGTCTGCATCATCATCTGCCGCGCCTTGTTCAACCGCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTtcgtc                                                                                                                            <  2:951714/145‑1 (MQ=255)
                            tgatCGACTTGCTGGTCTGCATCAGCATCTGCCGCGCCTTGTTCAACCGCAACTCCAGGTAGTACTGGCTGGGTACTCGGTTAAGGTACTGCTAGAGGATACCCTCCAGCTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCgg                                                                                                                          >  1:1377456/1‑150 (MQ=255)
                               tCGACTTGCTGGTCTGCATCAGCATCTGCCGCGCCTTGTTCAACCGCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCggtgg                                                                                                                       >  1:1328477/1‑150 (MQ=255)
                                    ttGCTGGTCTGCATCAGCATCTGCCGCGCCTTGTTCAACCGCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGc                                                                                                                  >  2:1578783/1‑150 (MQ=255)
                                                  cAGCATCTGCCGCGCCTTGTTCAACCGCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCACCTGCCGACGCGACACACACACATGCTGGGCCATTTCGTCGGTGGTCAGCGGCTCTTCGATGtt                                                                                                    <  2:1071498/150‑1 (MQ=255)
                                                    gCATCTGCCGCGCCTTGTTCAACCGCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTgg                                                                                                  <  1:389528/150‑1 (MQ=255)
                                                       tCTGCCGCGCCTTGTTCAACCGCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCTGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCtt                                                                                               <  1:1665244/150‑1 (MQ=255)
                                                       tCTGCCGCGCCTTGTTCAACCGCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCtt                                                                                               >  2:1677185/1‑150 (MQ=255)
                                                             gcgcCTTGTTCAACCGCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACACGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCa                                                                                         <  2:1572111/150‑1 (MQ=255)
                                                              cgcCTTGTTCAACCGCAGCTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACGCATGCTGGGCGATTTCGTCGGTCGTCAGCGGCTCTTCGATGTTGGCTTCCATCAg                                                                                        <  2:1070700/150‑1 (MQ=255)
                                                                 cTTGTTCAACCGCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCaa                                                                                     >  1:338577/1‑150 (MQ=255)
                                                                    gTTCAACCGCAACTCCAGGTAGTACTGGCTGGGCACGCGCTTCAGGTACTGCTTGAAGATACGCTCCAACGGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAacac                                                                                  >  2:595293/1‑150 (MQ=255)
                                                                     ttCAACCGCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGAt                                                                                                       >  1:1891542/1‑128 (MQ=255)
                                                                     ttCAACCGCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGAt                                                                                                       <  2:1891542/128‑1 (MQ=255)
                                                                     ttCAACCGCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACg                                                                                 >  1:771325/1‑150 (MQ=255)
                                                                     ttCAACCGCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGAAacagg                                                                                 >  1:868492/1‑148 (MQ=255)
                                                                     ttCAACCGCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCG‑‑ACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTc                                                                                                          >  1:1092515/1‑123 (MQ=255)
                                                                     ttCAACCGCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCG‑‑ACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTc                                                                                                          <  2:1092515/123‑1 (MQ=255)
                                                                       cAACCGCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGc                                                                                                                                     <  2:505992/96‑1 (MQ=255)
                                                                       cAACCGCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGc                                                                                                                                     >  1:505992/1‑96 (MQ=255)
                                                                        aaCCGCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGcc                                                                              >  2:1137373/1‑150 (MQ=255)
                                                                         aCCGCAACTGCAGGTAGGACTGGCTGGGCACTCGGTTCAGGTACTGCTTGAAGATTCGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCt                                                                             <  2:861939/150‑1 (MQ=255)
                                                                         aCCGCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTCCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCt                                                                             <  2:851246/150‑1 (MQ=255)
                                                                          ccGCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTg                                                                            <  2:1377456/150‑1 (MQ=255)
                                                                           cGCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGa                                                                           <  1:721648/150‑1 (MQ=255)
                                                                            gccaCTCCAGGTAATGCTGGCTGGGCATGCGGGGCAAGTACTGCATGAAGATACGCCCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGAg                                                                          >  2:695995/4‑150 (MQ=255)
                                                                            gCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCG‑‑ACACACACATGCTGGGCGATTTCGTCggtg                                                                                                                        >  1:1289109/1‑102 (MQ=255)
                                                                            gCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCG‑‑ACACACACATGCTGGGCGATTTCGTCggtg                                                                                                                        <  2:1289109/102‑1 (MQ=255)
                                                                            gCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCG‑‑ACACACACATGCTGGGCGATTTCGTCggtg                                                                                                                        <  2:1289048/102‑1 (MQ=255)
                                                                            gCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCG‑‑ACACACACATGCTGGGCGATTTCGTCggtg                                                                                                                        >  1:1289048/1‑102 (MQ=255)
                                                                                       tagtaCTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGAGTCAGCTTgggg                                                               <  1:1677185/150‑1 (MQ=255)
                                                                                             tggctggGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGAGTCAGCTTGGGGTGACTc                                                         <  1:663909/150‑1 (MQ=255)
                                                                                               gctggGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGtt                                                                                                    >  1:1047651/1‑105 (MQ=255)
                                                                                               gctggGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGtt                                                                                                    <  2:1047651/105‑1 (MQ=255)
                                                                                               gctggGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGAGTCAGCTTGGGGTGACTCGa                                                       <  1:170988/150‑1 (MQ=255)
                                                                                                 tggGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGAGTCAGCTTGGGGTGACTCGAGc                                                     >  1:1431940/1‑150 (MQ=255)
                                                                                                  gggCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGAGTCAGCTTGGGGTGACTCGAGcg                                                    <  2:1328477/150‑2 (MQ=255)
                                                                                                        gcgGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGAGTCAGCTTGGGGTGACTCGAGCCCAGGCg                                              >  1:30920/1‑150 (MQ=255)
                                                                                                           gTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGAGTCAGCTTGGGGTGACTCGAGCCCAGGCGGtt                                           <  1:541501/150‑1 (MQ=255)
                                                                                                           gTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGAGTCAGCTTGGGGTGACTCGAGCCCAGGCGGtt                                           <  2:1431940/150‑1 (MQ=255)
                                                                                                            ttCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCTACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGAGTCAGCTTGGGGTGACTCGAGCCCAGGCGGTTc                                          >  1:981755/1‑150 (MQ=255)
                                                                                                               aGGTACTGCTTGAAGATACGCCCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGAGTCAGCTTGGGGTGACTCGAGCCCAGGCGGTTCTGc                                       >  1:1492575/1‑150 (MQ=255)
                                                                                                                ggTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGAGTCAGCTTGGGGTGACTCGAGCCCAGGCGGTTCTGCa                                      >  1:866165/1‑150 (MQ=255)
                                                                                                                 gTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGAGTCAGCTTGGGGTGACTCGAGCCCAGGCGGTTCTGCAg                                     >  2:1293587/1‑150 (MQ=255)
                                                                                                                  tACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGAGTCAGCTTGGGGTGACTCGAGCCCAGGCGGTTCTGCAGc                                    >  2:1322756/1‑150 (MQ=255)
                                                                                                                  tACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGAGTCAGCTTGGGGTGACTCGAGCCCAGGCGGTTCTGCAGc                                    <  1:309437/150‑1 (MQ=255)
                                                                                                                  tACTGCTCGAAGTCACGCTCCAACTGCCGGCGCGACACACACTCATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGAGTCAGCTTGGGGTGACTCGAGCCCAGGCGGTTCTGCAGc                                    <  2:389058/150‑1 (MQ=255)
                                                                                                                     tGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTCGGTTGCCTTGGGGTGACTCGAGCCCAGGCGGTTCTGCAGCGGa                                 >  1:960272/1‑150 (MQ=255)
                                                                                                                        ttGNAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGAGTCAGCTTGGGGTGACTCGAGCCCAGGCGGTTCTGCAGCGGAATa                              >  2:1194827/1‑150 (MQ=255)
                                                                                                                        ttGAAGATACGCTCCACCTGCCGCCGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGAGTCAGCTTGGGGTGACTCGAGCCCAGGCGGTTCTGCAGCGGAATa                              >  2:695191/1‑150 (MQ=255)
                                                                                                                        ttGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGAGTCAGCTTGGGGTGACTCGAGCCCAGGCGGTTCTGCAGCGGAATa                              >  2:1898639/1‑150 (MQ=255)
                                                                                                                        ttGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGAGTCAGCTTGGGGTGACTCGAGCCCAGGCGGTTCTGCAGCGGAATa                              >  1:1339611/1‑150 (MQ=255)
                                                                                                                             gATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGAGTCAGCTTGGGGTGACTCGAGCCCAGCGGTTCTGCAGCGGAATACGCTgg                        >  2:1351635/1‑150 (MQ=255)
                                                                                                                                aCGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAacac                                                                                  <  1:1414677/90‑1 (MQ=255)
                                                                                                                                aCGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAacac                                                                                  >  2:1414677/1‑90 (MQ=255)
                                                                                                                                aCGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGAGTCAGCTTGGGGTGACTCGAGCCCAGGCGGTTCTGCAGCGGAATACGCTGgcg                      >  1:360508/1‑150 (MQ=255)
                                                                                                                                  gCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCAGCAGCTATACGGCCTGAGTCAGCTTGCGGTGAGGCGATCCCAGGCGGTTCTGCAGCTGAATACGCTGgcgca                    >  2:1842612/1‑149 (MQ=255)
                                                                                                                                    tCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGAGTCAGCTTGGGGTGACTCGAGCCCAGGCGGTTCTGCAGCGGAATACGCTGGCGCTCg                  >  1:1777979/1‑150 (MQ=255)
                                                                                                                                    tCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTAGATGTTGGCTTCCATCAGCAACACGGCCTGAGTCAGCTTGGGGTGACTCGAGCCCAGGCGGTTCTGCAGCGGAATACGCTGGCGCTCg                  >  1:1776493/1‑150 (MQ=255)
                                                                                                                                       aaCTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGAGTCAGCTTGGGGTGACTCGAGCCCAGGCGGTTCTGCAGCGGAATACGCTGGCGCTCGccc               >  1:747480/1‑149 (MQ=255)
                                                                                                                                       aaCTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGAGTCAGCTTGGGGTGACTCGAGCCCAGGCGGTTCTGCAGCGGAATACGCTGGCGCTCGCCa               <  2:338577/150‑1 (MQ=255)
                                                                                                                                        aCTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGAGTCAGCTTGGGGTGACTCGAGCCCAGGCGGTTCTGCAGCGGAATACGCTGGCGCTCGCCAc              <  2:1339611/150‑1 (MQ=255)
                                                                                                                                                  cgACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGAGTCAGCTTGGGGTGACTCGAGCCCAGGCGGTTCTGCAGCGGAATACGCTGGCGCTCGCCACCCTCGCGGAt    <  2:771325/150‑1 (MQ=255)
                                                                                                                                                   gACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGAGTCAGCTTGGGGTGACTCGAGCCCAGGCGGTTCTGCAGCGGAATACGCTGGCGCTCGCCACCCTCGCGGATg   >  1:817558/1‑150 (MQ=255)
                                                                                                                                                    acacacacacATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGAGTCAGCTTGGGGTGACTCGAGCCCAGGCGGTTCTGCAGCGGAATACGCTGGCGCTCGCCACCCTCGCGGATgc  <  1:695995/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GAGAAGCCACAGGACAGCCCGATCTGGATGATCGACTTGCTGGTCTGCATCAGCATCTGCCGCGCCTTGTTCAACCGCAACTCCAGGTAGTACTGGCTGGGCACGCGGTTCAGGTACTGCTTGAAGATACGCTCCAACTGCCGACGCGACACACACACATGCTGGGCGATTTCGTCGGTGGTCAGCGGCTCTTCGATGTTGGCTTCCATCAGCAACACGGCCTGAGTCAGCTTGGGGTGACTCGAGCCCAGGCGGTTCTGCAGCGGAATACGCTGGCGCTCGCCACCCTCGCGGATGC  >  NC_002947/5092813‑5093110

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: