Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,865,969 A→G 8.0% A303A (GCT→GCC) ‡ PP_3415 ← LacI family transcriptional regulator

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,865,9690AG 8.0% 173.2 / 5.5 62A303A (GCT→GCC) ‡PP_3415LacI family transcriptional regulator
Reads supporting (aligned to +/- strand):  ref base A (30/27);  new base G (2/3);  total (32/30)
Fisher's exact test for biased strand distribution p-value = 6.67e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.95e-01

CTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTC  >  NC_002947/3865821‑3866112
                                                                                                                                                     |                                                                                                                                               
cTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCGTGTGCCGCGCCCCGCCCCACCGCa                                                                                                                                                 <  7:185611/149‑1 (MQ=255)
  gtgGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCACGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCGTTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGc                                                                                                                                               >  7:250383/1‑149 (MQ=255)
   tgGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGc                                                                                                                                               <  3:51033/148‑1 (MQ=255)
           gcgcgAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGgtgt                                                                                                                                      <  4:52553/149‑1 (MQ=255)
           gcgcgAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGgtgt                                                                                                                                      >  4:123210/1‑149 (MQ=255)
           gcgcgAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGgtgt                                                                                                                                      <  2:337413/149‑1 (MQ=255)
            cgcgAGAAGGCCTTGGCGGGCACTCTCAGGTTCTCTCGCGCACCATCAGCTCAAAGACCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCACCTGTGCCGCGCCCCGCCCCACCGCGGGCTGCGgtggg                                                                                                                                     <  4:142677/149‑3 (MQ=255)
                   aGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATgg                                                                                                                              <  8:29843/149‑1 (MQ=255)
                     gacatGGCGGGCACTCTCAGGTGCTCTCGCCCACCAACACCACAAAGCCCAGATCCTGCTGCGCATTCGTCACCAGCTTCCCATCCAGCAACGCCACCACCGCCTGTGCCGCCCCCCGCCCCACCGCATCCCGCCGTGTGCGAATGGAg                                                                                                                            >  4:311929/5‑149 (MQ=255)
                             gggCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCg                                                                                                                    >  7:264284/1‑149 (MQ=255)
                                  ctctCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCa                                                                                                               <  5:95157/149‑1 (MQ=255)
                                   tctcAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCAc                                                                                                              >  5:341176/1‑149 (MQ=255)
                                   tctcAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCAc                                                                                                              >  5:73213/1‑149 (MQ=255)
                                    ctcAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCAc                                                                                                              >  3:130795/1‑148 (MQ=255)
                                       aGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCAcc                                                                                                             <  3:311929/146‑1 (MQ=255)
                                             tctcGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGcc                                                                                                    >  2:1715/1‑149 (MQ=255)
                                             tctcGCGCACCATCAGCTCAAAGCCAACAACCTGCTGCGCAATCGTCACCCGCTTCCCATCCAGCAACGCAAGCAGCCCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATCCGcc                                                                                                    >  5:271765/1‑149 (MQ=255)
                                                 gcgcACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAgg                                                                                                >  8:207777/1‑149 (MQ=255)
                                                    cACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCg                                                                                             >  1:166469/1‑149 (MQ=255)
                                                     aCCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCgg                                                                                            <  3:123210/149‑1 (MQ=255)
                                                       cATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCgg                                                                                            >  1:143206/1‑147 (MQ=255)
                                                         tCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGgaga                                                                                        <  7:207777/149‑1 (MQ=255)
                                                              tCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGCTCGtt                                                                                   >  5:273411/1‑149 (MQ=255)
                                                              tCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGtt                                                                                   >  8:116834/1‑149 (MQ=255)
                                                              tCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCGCCCCACCGCGGGCTGCGgtgggg                                                                                                                                    >  4:350511/1‑96 (MQ=255)
                                                              tCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCGCCCCACCGCGGGCTGCGgtgggg                                                                                                                                    <  3:350511/99‑4 (MQ=255)
                                                               cAAAGCCCAGATCCTGCTGCGCATTCGGCACCCGCTTCCCTTCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGTTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGtt                                                                                   <  6:271765/148‑1 (MQ=255)
                                                                aaaGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTTCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTg                                                                                  <  1:276732/148‑1 (MQ=255)
                                                                aaaGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCTCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGtt                                                                                   <  7:314236/147‑1 (MQ=255)
                                                                aaaGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCGGGCTGCGgtgggg                                                                                                                                    <  7:77074/98‑4 (MQ=255)
                                                                aaaGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCGGGCTGCGgtgggg                                                                                                                                    >  8:77074/1‑95 (MQ=255)
                                                                aaaGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTg                                                                                  >  5:147809/1‑148 (MQ=255)
                                                                  aGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGaaa                                                                               >  8:176709/1‑149 (MQ=255)
                                                                   gCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGaaa                                                                               >  1:96571/1‑149 (MQ=255)
                                                                                    cATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTgg                                                             >  8:284187/1‑149 (MQ=255)
                                                                                     aTTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGc                                                            <  6:273411/149‑1 (MQ=255)
                                                                                     aTTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGc                                                            >  6:67254/1‑149 (MQ=255)
                                                                                          tCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCACAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCgg                                                                                            <  1:165205/112‑1 (MQ=255)
                                                                                          tCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCACAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCgg                                                                                            >  2:165205/1‑112 (MQ=255)
                                                                                             cccGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCAcc                                                    >  5:248398/1‑149 (MQ=255)
                                                                                               cGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTc                                                  <  8:264284/149‑1 (MQ=255)
                                                                                                     ccATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCAcccc                                            <  1:205186/149‑1 (MQ=255)
                                                                                                     ccATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCAcccc                                            >  2:205186/1‑149 (MQ=255)
                                                                                                      cATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCggg                                                                                           <  1:265108/101‑1 (MQ=255)
                                                                                                      cATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCggg                                                                                           >  2:265108/1‑101 (MQ=255)
                                                                                                      cATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCt                                           >  5:46363/1‑149 (MQ=255)
                                                                                                       aTCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCtt                                          >  8:5665/1‑149 (MQ=255)
                                                                                                        tCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTg                                         <  4:130795/149‑1 (MQ=255)
                                                                                                         ccAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGa                                        >  6:9799/1‑149 (MQ=255)
                                                                                                          cAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGAc                                       <  1:1715/149‑1 (MQ=255)
                                                                                                          cAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGAc                                       <  7:284187/149‑1 (MQ=255)
                                                                                                            gCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCAcccc                                            >  4:35331/1‑142 (MQ=255)
                                                                                                            gCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCAcccc                                            <  3:35331/142‑1 (MQ=255)
                                                                                                               aCGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAg                                  >  3:286918/1‑149 (MQ=255)
                                                                                                               aCGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAg                                  <  4:286918/149‑1 (MQ=255)
                                                                                                                     gcagcGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTgc                            <  6:73213/149‑1 (MQ=255)
                                                                                                                      cagcGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTg                             >  1:5026/1‑147 (MQ=255)
                                                                                                                           gcCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCAc                       <  7:116834/149‑1 (MQ=255)
                                                                                                                                  ccGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCAGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGccc                >  8:334777/1‑149 (MQ=255)
                                                                                                                                     cgccccgcccCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGa             <  5:67254/149‑1 (MQ=255)
                                                                                                                                          cgccccACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAg        <  6:46363/149‑1 (MQ=255)
                                                                                                                                           gccccACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGACCTggctcagg       <  2:166469/149‑8 (MQ=255)
                                                                                                                                                aCCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGgtcgtc  <  6:248398/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                               
CTGTGGGAGTCGCGCGAGAAGGCCTTGGCGGGCACTCTCAGGTGCTCTCGCGCACCATCAGCTCAAAGCCCAGATCCTGCTGCGCATTCGTCACCCGCTTCCCATCCAGCAACGCCAGCAGCGCCTGTGCCGCGCCCCGCCCCACCGCAGCCCGCGGTGTGCGAATGGAGGTCAGCCGTGGCACCATATGCGCCGAGGCCGGGAGATCGTTGAAACCGACCATCGCCACCTGGCGTGGCACCTCCACCCCTTGACGCAAGGCCTGCAGCACCGCGCCCTGAGCCAGGTCGTC  >  NC_002947/3865821‑3866112

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: