Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,857,961 T→A 5.4% L170L (CTT→CTA)  PP_5134 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,857,9610TA 5.4% 229.7 / 2.6 74L170L (CTT→CTA) PP_5134hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (39/31);  new base A (2/2);  total (41/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.27e-01

GTGATCGATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACCAGCAGGCACAG  >  NC_002947/5857821‑5858108
                                                                                                                                            |                                                                                                                                                   
gtgATCGATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCt                                                                                                                                             >  5:299526/1‑149 (MQ=255)
 tgATCGATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCtt                                                                                                                                            <  6:299526/149‑1 (MQ=255)
      gATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGc                                                                                                                                        <  6:67102/148‑1 (MQ=255)
        tGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCtt                                                                                                                                      >  3:111539/1‑148 (MQ=255)
        tGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCAAGCCGGGCAAGGGCTTc                                                                                                                                     <  2:251339/149‑1 (MQ=255)
             gcgAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCa                                                                                                                                >  7:137647/1‑149 (MQ=255)
              cgAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCAt                                                                                                                               >  7:137209/1‑149 (MQ=255)
                       gCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGt                                                                                                                      >  3:322070/1‑149 (MQ=255)
                         tGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTc                                                                                                                    >  7:226891/1‑149 (MQ=255)
                         tGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTAGCCCGGCTTGGGCTTCACGCATTGGTACGGTTc                                                                                                                    >  4:75821/1‑149 (MQ=255)
                            ccAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGa                                                                                                                  >  7:356462/1‑148 (MQ=255)
                             cAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGa                                                                                                                  >  4:345428/1‑147 (MQ=255)
                                 tACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGt                                                                                                            >  4:40089/1‑149 (MQ=255)
                                   ccGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAg                                                                                                               <  5:278098/144‑1 (MQ=255)
                                   ccGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAg                                                                                                               >  6:278098/1‑144 (MQ=255)
                                    cgcgACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGc                                                                                                         >  6:359032/1‑149 (MQ=255)
                                    cgcgACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGc                                                                                                         <  8:226891/149‑1 (MQ=255)
                                     gcgACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGcc                                                                                                        >  6:72474/1‑149 (MQ=255)
                                      cgacgaACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCa                                                                                                                                    <  6:48647/120‑1 (MQ=255)
                                      cgacgaACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCa                                                                                                                                    <  6:42878/120‑1 (MQ=255)
                                      cgacgaACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCa                                                                                                                                    >  5:48647/1‑120 (MQ=255)
                                      cgacgaACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCa                                                                                                                                    >  5:42878/1‑120 (MQ=255)
                                      cgacgaACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGccc                                                                                                       >  4:267196/1‑149 (MQ=255)
                                       gacgaACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTAGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGcccc                                                                                                      <  3:75821/149‑1 (MQ=255)
                                            aCTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCag                                                                                                 <  8:356462/149‑1 (MQ=255)
                                             cTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCaga                                                                                                <  3:33802/149‑1 (MQ=255)
                                             cTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCa                                                                                                  >  8:50039/1‑147 (MQ=255)
                                               ggAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCagag                                                                                               >  3:184405/1‑148 (MQ=255)
                                                gAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCagag                                                                                               <  5:359032/147‑1 (MQ=255)
                                                  aGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATg                                                                                           <  4:165725/149‑1 (MQ=255)
                                                   gTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATgg                                                                                          <  3:158535/149‑1 (MQ=255)
                                                       tGGCCATCTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCagag                                                                                               <  5:16541/140‑1 (MQ=255)
                                                       tGGCCAGCTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCagag                                                                                               >  6:16541/1‑140 (MQ=255)
                                                           cAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGcc                                                                                  <  8:137209/149‑1 (MQ=255)
                                                              cTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGtt                                                                               <  4:322070/149‑1 (MQ=255)
                                                                 cGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTg                                                                             >  6:87828/1‑148 (MQ=255)
                                                                   gCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCtt                                                                                                                                      <  1:132070/89‑1 (MQ=255)
                                                                   gCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCtt                                                                                                                                      >  2:132070/1‑89 (MQ=255)
                                                                         cGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGt                                                                                         >  8:170465/1‑128 (MQ=255)
                                                                         cGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGt                                                                                         <  7:170465/128‑1 (MQ=255)
                                                                         cGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGc                                                                    >  3:82309/1‑149 (MQ=255)
                                                                          gTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGc                                                                                                                                              <  3:321422/74‑1 (MQ=255)
                                                                          gTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGc                                                                                                                                              >  4:321422/1‑74 (MQ=255)
                                                                            tgctgcCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTg                                                                 <  2:143431/149‑1 (MQ=255)
                                                                                gcCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGgcgc                                                             >  8:97761/1‑149 (MQ=255)
                                                                                   ggTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTg                                                                 <  1:127674/142‑1 (MQ=255)
                                                                                   ggTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTg                                                                 >  2:127674/1‑142 (MQ=255)
                                                                                   ggTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCgctgc                                                          <  3:136076/149‑1 (MQ=255)
                                                                                       cTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCgctgctgct                                                      <  4:111539/149‑1 (MQ=255)
                                                                                        tCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGgcg                                                              >  4:14225/1‑140 (MQ=255)
                                                                                        tCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGgcg                                                              <  3:14225/140‑1 (MQ=255)
                                                                                          aaCTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGt                                                                                         <  5:262669/111‑1 (MQ=255)
                                                                                          aaCTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGt                                                                                         >  6:262669/1‑111 (MQ=255)
                                                                                          aaCTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTgg                                                    >  2:230224/1‑148 (MQ=255)
                                                                                          aaCTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTgg                                                    >  6:178312/1‑148 (MQ=255)
                                                                                                gtcgCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAAcgc                                             >  6:84301/1‑149 (MQ=255)
                                                                                                   gccagcagccagcaCCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAAcgccc                                           <  4:82309/148‑1 (MQ=255)
                                                                                                       gcagccagcaCCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGcccgcccg                                      >  7:313925/1‑149 (MQ=255)
                                                                                                        cagccagcaCCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCgg                                     <  5:72474/149‑1 (MQ=255)
                                                                                                             agcaCCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTg                                >  8:325375/1‑149 (MQ=255)
                                                                                                              gcaCCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTg                                >  1:239641/1‑148 (MQ=255)
                                                                                                                  cgcgAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACg                           >  5:5463/1‑149 (MQ=255)
                                                                                                                   gcgAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACgg                          <  7:325375/149‑1 (MQ=255)
                                                                                                                    cgAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGc                         <  5:87828/149‑1 (MQ=255)
                                                                                                                       gCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGa                      <  3:40089/149‑1 (MQ=255)
                                                                                                                        ccTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGATCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGAt                     >  1:93128/1‑149 (MQ=255)
                                                                                                                         cTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGaaa                                                 <  4:274919/120‑1 (MQ=255)
                                                                                                                         cTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGaaa                                                 >  3:274919/1‑120 (MQ=255)
                                                                                                                                 cGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCgacgac              >  2:34018/1‑147 (MQ=255)
                                                                                                                                     aaGCCCAAGCCGGGCAAGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACcagcag        >  3:162162/1‑149 (MQ=255)
                                                                                                                                      aGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACCAGCAgg       <  3:267196/149‑1 (MQ=255)
                                                                                                                                       gCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACCAGCAGGc      >  2:238569/1‑149 (MQ=255)
                                                                                                                                         ccTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACcagca         >  2:224935/1‑144 (MQ=255)
                                                                                                                                         ccTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACcagca         <  1:224935/144‑1 (MQ=255)
                                                                                                                                           ttGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACCAGCAGGCACAg  >  3:242412/1‑149 (MQ=255)
                                                                                                                                            |                                                                                                                                                   
GTGATCGATGCGCGCGAGCCGGTGCTGGCCAAGTACCGCGACGAACTGGAAGTGCTGGCCAACTGCGGCAAGCCGTTGCTGCCGGTGCTCAACTTCGTCGCCAGCAGCCAGCACCGCGAGCCTGAATGGCGTGAAGCCCTTGCCCGGCTTGGGCTTCACGCATTGGTACGGTTCGACAGCGTGGCCCCGCCAGAGGATGGTGAGCGCCGTTTGTACGAAAGCCTGGCGCTGCTGCTGGAAAACGCCCGCCCGGCGCTGCAACGGCTGATCGACGACCAGCAGGCACAG  >  NC_002947/5857821‑5858108

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: ‑