Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,522,791 T→A 6.0% intergenic (+25/+15) PP_4856 → / ← yhjG Dps family ferritin/protein YhjG

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,522,7910TA 6.0% 199.7 / 2.7 67intergenic (+25/+15)PP_4856/yhjGDps family ferritin/protein YhjG
Reads supporting (aligned to +/- strand):  ref base T (30/33);  new base A (2/2);  total (32/35)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.26e-01

CTACCGCGAGATCATTCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTG  >  NC_002947/5522643‑5522933
                                                                                                                                                    |                                                                                                                                               
cTACCGCGAGATCATTCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCt                                                                                                                                                 >  7:31919/1‑149 (MQ=255)
cTACCGCGAGATCATTCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCt                                                                                                                                                 >  8:232760/1‑149 (MQ=255)
 tACCGCGAGATCATTCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCtt                                                                                                                                                >  4:345041/1‑149 (MQ=255)
  aCCGCGAGATCATTCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCttt                                                                                                                                               >  4:50147/1‑149 (MQ=255)
     gcgAGATCATTCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGc                                                                                                                                            <  4:350502/149‑1 (MQ=255)
           tCATTCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGcccc                                                                                                                                      >  8:160850/1‑149 (MQ=255)
              ttCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTg                                                                                                                                   >  1:180712/1‑149 (MQ=255)
              ttCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTg                                                                                                                                   <  3:50147/149‑1 (MQ=255)
              ttCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACAACCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCAGACCTGTTGCAAGGCCTTTGATTTCTGACACCCACAAGGGGCCGCAAAAGCGGCCCCTTg                                                                                                                                   >  3:162161/1‑149 (MQ=255)
                 aaTACATCGTGGATAAAGCACCGACTACCCGTCTTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGc                                                                                                                                            <  8:293613/137‑1 (MQ=255)
                 aaTACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACTCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGc                                                                                                                                            >  7:293613/1‑137 (MQ=255)
                 aaTACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTa                                                                                                                                <  1:264519/149‑1 (MQ=255)
                  aTACTTCGGGGATAATGACCCGACTACGCGTCGTTTCTTCGAAGACATCCTGGCGCAGGAATAGGCGCACGCGTACGATATGTCCGCCGTGTTGCAATGCCTTTGATTTCTGACACCCTCAAGGGGCCGCAAAAGCGGCCCCTTgaggg                                                                                                                               <  4:162161/149‑5 (MQ=255)
                         ggggATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGt                                                                                                                         <  7:68723/148‑1 (MQ=255)
                          gggATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGc                                                                                                                       <  2:180712/149‑1 (MQ=255)
                            gATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGcccc                                                                                                                                      <  1:65830/132‑1 (MQ=255)
                            gATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGcccc                                                                                                                                      >  2:65830/1‑132 (MQ=255)
                                aaGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCa                                                                                                                  <  7:174471/148‑1 (MQ=255)
                                  gACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCAccg                                                                                                               >  3:281919/1‑149 (MQ=255)
                                   aCCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCAccac                                                                                                              >  5:308501/1‑147 (MQ=255)
                                        aCTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCg                                                                                                         >  7:70649/1‑149 (MQ=255)
                                         cTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCgg                                                                                                        >  5:333395/1‑149 (MQ=255)
                                            cccGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGgcgc                                                                                                     >  3:301246/1‑149 (MQ=255)
                                             ccGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGcc                                                                                                    >  4:148021/1‑149 (MQ=255)
                                             ccGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGcc                                                                                                    >  4:210675/1‑149 (MQ=255)
                                              cgtcgtATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCt                                                                                                   <  5:269664/149‑1 (MQ=255)
                                                       ttCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTgcctgcct                                                                                          >  3:24555/1‑149 (MQ=255)
                                                       ttCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTgcctgcct                                                                                          >  2:290123/1‑149 (MQ=255)
                                                       ttCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTgcctgcct                                                                                          >  8:300900/1‑149 (MQ=255)
                                                                aTCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATctgct                                                                                >  6:190780/1‑149 (MQ=255)
                                                                      gcgcAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGcaaca                                                                          <  8:70649/149‑1 (MQ=255)
                                                                      gcgcAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCACCGTCGGCGCCTTGCCTGCCTTCATCTGCTGcaac                                                                           <  6:308501/148‑1 (MQ=255)
                                                                       cgcAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGcaacaa                                                                         <  2:16893/149‑1 (MQ=255)
                                                                        gcAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCaacaac                                                                        >  3:70016/1‑149 (MQ=255)
                                                                        gcAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCAAAAGCGGCCCCTTgagg                                                                                                                                <  1:157904/94‑4 (MQ=38)
                                                                        gcAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCAAAAGCGGCCCCTTgagg                                                                                                                                >  2:157904/1‑91 (MQ=38)
                                                                           ggAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGgt                                                                     >  4:343587/1‑149 (MQ=255)
                                                                           ggAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGgt                                                                     <  7:160850/149‑1 (MQ=255)
                                                                           ggAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGgt                                                                     <  8:31919/149‑1 (MQ=255)
                                                                                  gagCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACt                                                              >  5:155724/1‑149 (MQ=255)
                                                                                  gagCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACt                                                              >  7:79291/1‑149 (MQ=255)
                                                                                  gagCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACt                                                              <  3:210675/149‑1 (MQ=255)
                                                                                       cgcgGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTc                                                         <  3:345041/149‑1 (MQ=255)
                                                                                          ggACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATctgctg                                                                               >  6:85437/1‑124 (MQ=255)
                                                                                          ggACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATctgctg                                                                               <  5:85437/124‑1 (MQ=255)
                                                                                               atatGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATc                                                 <  3:92908/149‑1 (MQ=255)
                                                                                                  tGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATccc                                               <  6:333395/148‑1 (MQ=255)
                                                                                                    tCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCgc                                            <  4:24555/149‑1 (MQ=255)
                                                                                                     ccGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCgct                                           <  7:300900/149‑1 (MQ=255)
                                                                                                       gACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGTTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCgctgc                                         >  5:141772/1‑149 (MQ=255)
                                                                                                          cTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTgg                                      <  8:79291/149‑1 (MQ=255)
                                                                                                             ttGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTggg                                     <  5:190780/147‑1 (MQ=255)
                                                                                                             ttGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGgc                                    <  4:281919/148‑1 (MQ=255)
                                                                                                                 aaGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCTCTGGTTCGGGTCATCCCCGCTGCTGGGCGCaa                                >  7:146821/1‑148 (MQ=255)
                                                                                                                     ccTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAgc                           <  4:301246/149‑1 (MQ=255)
                                                                                                                       tttGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAgcgc                         >  7:174608/1‑149 (MQ=255)
                                                                                                                       tttGATTTCTGACACCCTCAAGAGGCCGCTTTTGCGGCCCCTTGCTATATCCGTGCTTTCACCGCCGTCGGCGCCTTGCCAGACTTCATCTGCAGCAACAACGGAGTGCACTCGTTCGGGTCATCACCGCTGCTGGGCGCAATCAgcgc                         >  4:284213/1‑149 (MQ=255)
                                                                                                                           aTTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGGCTTGCCTGCCTTCATCTGTATGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCa                       <  1:290123/148‑1 (MQ=255)
                                                                                                                                  acacCCTCCACGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACCGTGTGCACTGTTTCGGGTCATCCCCGCTGCTTGGCGCAATCAGCGCCAGCAAccc                <  7:232760/147‑1 (MQ=255)
                                                                                                                                  acacCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCccg              >  6:334557/1‑149 (MQ=255)
                                                                                                                                   cacCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCAGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCccgc             >  4:331476/1‑149 (MQ=255)
                                                                                                                                      ccTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGgcg                                   <  1:11358/124‑1 (MQ=255)
                                                                                                                                      ccTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGgcg                                   >  2:11358/1‑124 (MQ=255)
                                                                                                                                      ccTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCccgccgc          >  1:240403/1‑149 (MQ=255)
                                                                                                                                       cTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCccgccgcc         >  3:362389/1‑149 (MQ=255)
                                                                                                                                        tCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTgg                                      <  2:185185/119‑1 (MQ=255)
                                                                                                                                        tCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTgg                                      >  1:185185/1‑119 (MQ=255)
                                                                                                                                          aaGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCccgccgccg        <  3:148021/147‑1 (MQ=255)
                                                                                                                                           aGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGcc     <  4:70016/149‑1 (MQ=255)
                                                                                                                                             gggCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCt   <  3:349302/149‑1 (MQ=255)
                                                                                                                                              ggCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTg  <  3:343587/149‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CTACCGCGAGATCATTCAATACATCGGGGATAAAGACCCGACTACCCGTCGTATCTTCGAAGACATCCTGGCGCAGGAAGAGGAGCACGCGGACGATATGTCCGACCTGTTGCAAGGCCTTTGATTTCTGACACCCTCAAGGGGCCGCTTTTGCGGCCCCTTGCTATTTCCGTGCTTTCACCGCCGTCGGCGCCTTGCCTGCCTTCATCTGCTGCAACAACGGTGTGCACTGGTTCGGGTCATCCCCGCTGCTGGGCGCAATCAGCGCCAGCAACCCCGCCGCCGGCCCTG  >  NC_002947/5522643‑5522933

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: