Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,117,850 T→C 8.8% H33H (CAT→CAC PP_4504 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,117,8500TC 8.8% 185.9 / 7.8 68H33H (CAT→CACPP_4504hypothetical protein
Reads supporting (aligned to +/- strand):  ref base T (33/29);  new base C (2/4);  total (35/33)
Fisher's exact test for biased strand distribution p-value = 4.21e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.52e-01

CTAGACTGTGCTCATCCGGTAATCACCGGCCTATGAAAGGAGAGACACCATGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGCAAATGCAGCAGCGCCGCGACATTACCGAGCGCTACCTGTCCAAGCTGCCGGCCACCGACCAGAAAGCCATG  >  NC_002947/5117703‑5117997
                                                                                                                                                   |                                                                                                                                                   
cTAGACTGTGCTCATCCGGTAATCACCGGCCTATGAAAGGAGAGACACCATGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCatcatc                                                                                                                                                    <  5:30973/149‑1 (MQ=255)
   gACTGTGCTCATCCGGTAATCACCGGCCTATGAAAGGAGAGACACCATGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATcgcg                                                                                                                                                 >  7:76384/1‑149 (MQ=255)
        tgCTCATCCGGTAATCACCGGCCTATGAAAGGAGAGACACCATGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGAt                                                                                                                                            <  7:202626/149‑1 (MQ=255)
          cTCATCCGGTAATCACCGGCCTATGAAAGGAGAGACACCATGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATg                                                                                                                                           <  5:117554/148‑1 (MQ=255)
             aTCCGGTAATCACCGGCCTATGACAGGAGAGACACCATGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCACCGCGATGATgcggg                                                                                                                                       <  4:240796/149‑3 (MQ=255)
             aTCCGGTAATCACCGGCCTATGAAAGGAGAGACACCATGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGTATCACCGCGATGATgcggg                                                                                                                                       <  8:273754/149‑3 (MQ=255)
                cGGTAATCACCGGCCTATGAAAGGAGAGACACCATGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCttt                                                                                                                                    <  7:64079/149‑1 (MQ=255)
                 ggTAATCACCGGCCTATGAAAGGAGAGACACCATGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCt                                                                                                                                      >  4:14662/1‑146 (MQ=255)
                 ggTAATCACCGGCCTATGAAAGGAGAGACACCATGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCt                                                                                                                                      <  3:14662/146‑1 (MQ=255)
                  gTAATCACCGGCCTATGAAAGGAGAGACACCATGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTc                                                                                                                                   >  2:303772/1‑148 (MQ=255)
                   tAATCACCGGCCTATGAAAGGAGAGACACCATGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATcg                                                                                                                                                   <  1:133323/131‑1 (MQ=255)
                   tAATCACCGGCCTATGAAAGGAGAGACACCATGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATcg                                                                                                                                                   >  2:133323/1‑131 (MQ=255)
                   tAATCACCGGCCTATGAAAGGAGAGACACCATGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCACGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTcca                                                                                                                                 <  7:307249/149‑1 (MQ=255)
                    aaTCACCGGCCTATGAAAGGAGAGACACCATGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTccac                                                                                                                                >  7:31226/1‑149 (MQ=255)
                       cACCGGCCTATGAAAGGAGAGACACCATGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCaa                                                                                                                             >  4:109508/1‑149 (MQ=255)
                                   aaaGGAGAGACACCATGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAg                                                                                                                  <  8:76384/148‑1 (MQ=255)
                                        agagaCACCATGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCcgcg                                                                                                             >  3:162642/1‑148 (MQ=255)
                                         gagaCACCATGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGAt                                                                                                           >  6:63323/1‑149 (MQ=255)
                                           gaCACCATGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGAt                                                                                                           >  3:265129/1‑147 (MQ=255)
                                           gaCACCATGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGAt                                                                                                           >  6:117059/1‑147 (MQ=255)
                                           gaCACCATGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCa                                                                                                         <  4:265129/149‑1 (MQ=255)
                                            acacCATGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAg                                                                                                        <  3:58669/149‑1 (MQ=255)
                                               ccATGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAgcgc                                                                                                     >  2:117215/1‑149 (MQ=255)
                                               ccATGCGCAAGACCCTTATCACCCTGATGTTCGCTACCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAgcgc                                                                                                     >  8:56229/1‑149 (MQ=255)
                                                 aTGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAAc                                                                                                                            >  8:334838/1‑124 (MQ=255)
                                                 aTGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAAc                                                                                                                            <  7:334838/124‑1 (MQ=255)
                                                     gcAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCcag                                                                                                  >  1:353551/1‑146 (MQ=255)
                                                            ccTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGc                                                                                         >  6:248325/1‑148 (MQ=255)
                                                            ccTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGc                                                                                         >  6:242882/1‑148 (MQ=255)
                                                            ccTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGc                                                                                         >  2:155839/1‑148 (MQ=255)
                                                             cTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCaa                                                                                       <  3:109508/149‑1 (MQ=255)
                                                                aTCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGc                                                                                     <  1:28005/148‑1 (MQ=255)
                                                                   gCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGa                                                                                  <  5:78291/148‑1 (MQ=255)
                                                                   gCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGAt                                                                                 >  4:111801/1‑149 (MQ=255)
                                                                    cccTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATg                                                                                >  5:49324/1‑149 (MQ=255)
                                                                      cTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAgcgc                                                                                                     >  5:168524/1‑126 (MQ=255)
                                                                      cTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAgcgc                                                                                                     <  6:168524/126‑1 (MQ=255)
                                                                         aTGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGa                                                                           <  8:31226/149‑1 (MQ=255)
                                                                            ttCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGc                                                                         <  6:49324/148‑1 (MQ=255)
                                                                               gCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGCAAAt                                                                     <  4:81952/149‑1 (MQ=255)
                                                                                cTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCcagc                                                                                                 >  7:342333/1‑120 (MQ=255)
                                                                                cTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCcagc                                                                                                 <  8:342333/120‑1 (MQ=255)
                                                                                         tGCCGACCGTAGCCATGGCCATGCCAGATGGCGGTCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGCGCCGCGATCAGCGCCAGCAGCGCGGCAAGCTGATGGGCGAGCAAATGCagcagc                                                             <  5:315701/147‑1 (MQ=255)
                                                                                         tGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGCAAATGCagcagc                                                             <  3:111801/147‑1 (MQ=255)
                                                                                              aCCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTccacc                                                                                                                               >  7:57085/1‑76 (MQ=255)
                                                                                              aCCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTccacc                                                                                                                               <  8:57085/76‑1 (MQ=255)
                                                                                                  tAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGCAAATGCAGCAGCGCCGCGACAt                                                   >  6:248970/1‑148 (MQ=255)
                                                                                                    gccatggccatgCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGCAAATGCAGCAGCGCCGCGACATTa                                                 <  5:63323/148‑1 (MQ=255)
                                                                                                         ggccatgCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAgc                                                                                                       <  8:45210/89‑1 (MQ=255)
                                                                                                         ggccatgCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAgc                                                                                                       >  7:45210/1‑89 (MQ=255)
                                                                                                             atgCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCca                                                                                                   <  7:201719/89‑1 (MQ=255)
                                                                                                             atgCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCca                                                                                                   >  8:201719/1‑89 (MQ=255)
                                                                                                              tgCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGCAAATGCAGCAGCGCCGCGACATTACCGAGCGCTAc                                      <  4:162642/149‑1 (MQ=255)
                                                                                                                  aGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGCAAATGCAGCAGCGCCGCGACATTACCGAGCGCTACCTGt                                  >  4:158734/1‑149 (MQ=255)
                                                                                                                       gcggcCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGCAAATGCAGCAGCGCCGCGACATTACCGAGCGCTACCTGTCCAAg                             >  5:80022/1‑149 (MQ=255)
                                                                                                                          gcCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGCAAATGCAGCAGCGCCGCGACATTACCGAGCGCTACCTGTCCAAGCTg                          >  8:259138/1‑149 (MQ=255)
                                                                                                                              gcgcCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGCAAATGCAGCAGCGCCGCGACATTACCGAGCGCTACCTGTCCaa                              >  8:70371/1‑141 (MQ=255)
                                                                                                                              gcgcCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGCAAATGCAGCAGCGCCGCGACATTACCGAGCGCTACCTGTCCaa                              <  7:70371/141‑1 (MQ=255)
                                                                                                                              gcgcCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGCAAATGCAGCAGCGCCGCGACATTACCGAGCGCTACCTGTCCAAGCTGCCgg                      <  5:117059/149‑1 (MQ=255)
                                                                                                                                gcCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGCAAATGCAGCAGCGCCGCGACATTACCGAGCGCTACCTGTCCAAGCTGCCGGc                     >  5:347783/1‑148 (MQ=255)
                                                                                                                                  cACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGCAAATGCAGCAGCGCCGCGACATTACCGAGCGCTACCTGTCCAAGCTGCCGGCCAc                  >  2:260679/1‑149 (MQ=255)
                                                                                                                                   acgacgGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGCAAATGCAGCAGCGCCGCGACATTACCGAGCGCTACCTGTCCAAGCTGCCGGCCAcc                 <  1:260679/149‑1 (MQ=255)
                                                                                                                                       aaaggcGCATCACCGCGATGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGCAAATGCAGCAGCGCCGCGACATTAc                                                <  2:131958/109‑1 (MQ=255)
                                                                                                                                       aaaggcGCATCACCGCGATGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGCAAATGCAGCAGCGCCGCGACATTAc                                                >  1:131958/6‑114 (MQ=255)
                                                                                                                                          cccGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGCAAATGCAGCAGCGCCGCGACATTACCGAGCGCTACCTGTCCAAGCTGCCGGCCACCGACCAGa          <  6:347783/149‑1 (MQ=255)
                                                                                                                                           ccGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGCAAATGCAGCAGCGCCGCGACATTACCGAGCGCTACCTGTCCAAGCTGCCGGCCACCGACCAGa          >  4:292775/1‑148 (MQ=255)
                                                                                                                                           ccGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGCAAATGCAGCAGCGCCGCGACATTACCGAGCGCTACCTGTCCAAGCTGCCGGCCACCGACCAGa          >  1:325722/1‑148 (MQ=255)
                                                                                                                                             gCATCACCGCGATGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGCAAATGCAGCAGCGCCGCGACATTACCGAGCGCTACCTGTCCAAGCTGCCGGCCACCGACCAGaaa        <  8:214482/148‑1 (MQ=255)
                                                                                                                                             gCATCACCGCGATGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGCAAATGCAGCAGCGCCGCGACATTACCGAGCGCTACCTGTCCAAGCTGCCGGCCACCGACCAGa          >  7:214482/1‑146 (MQ=255)
                                                                                                                                                  atcGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGCAAATGCAGCAGCGCCGCGACATTACCGAGCGCTACCTGTCCAAGCTGCCGGCCACCGACCAGAAAGCCATg  <  2:325722/149‑1 (MQ=255)
                                                                                                                                                   tcGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGCAAATGCAGCAGCGCCGCGACATTACCGAGCGCTACCTGTCCAAGCTGCCGGCCACCGACCAg           <  1:25739/139‑1 (MQ=255)
                                                                                                                                                   tcGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGCAAATGCAGCAGCGCCGCGACATTACCGAGCGCTACCTGTCCAAGCTGCCGGCCACCGACCAg           >  2:25739/1‑139 (MQ=255)
                                                                                                                                                   |                                                                                                                                                   
CTAGACTGTGCTCATCCGGTAATCACCGGCCTATGAAAGGAGAGACACCATGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGCAAATGCAGCAGCGCCGCGACATTACCGAGCGCTACCTGTCCAAGCTGCCGGCCACCGACCAGAAAGCCATG  >  NC_002947/5117703‑5117997

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: