Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 5,117,855 G→A 8.8% G35D (GGT→GAT)  PP_4504 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029475,117,8550GA 8.8% 190.8 / 6.5 68G35D (GGT→GAT) PP_4504hypothetical protein
Reads supporting (aligned to +/- strand):  ref base G (32/30);  new base A (2/4);  total (34/34)
Fisher's exact test for biased strand distribution p-value = 6.73e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.22e-01

TGCTCATCCGGTAATCACCGGCCTATGAAAGGAGAGACACCATGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGCAAATGCAGCAGCGCCGCGACATTACCGAGCGCTACCTGTCCAAGCTGCCGGCCACCGACCAGAAAGCCATGAAGGAT  >  NC_002947/5117711‑5118003
                                                                                                                                                |                                                                                                                                                    
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                                                                                                                                            <  3:14662/146‑1 (MQ=255)
         ggTAATCACCGGCCTATGAAAGGAGAGACACCATGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCt                                                                                                                                            >  4:14662/1‑146 (MQ=255)
          gTAATCACCGGCCTATGAAAGGAGAGACACCATGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTc                                                                                                                                         >  2:303772/1‑148 (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                                                                                                                 >  6:117059/1‑147 (MQ=255)
                                   gaCACCATGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGAt                                                                                                                 >  3:265129/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                                                                                                                                     <  8:57085/76‑1 (MQ=255)
                                                                                      aCCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTccacc                                                                                                                                     >  7:57085/1‑76 (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                                                                                                         >  8:201719/1‑89 (MQ=255)
                                                                                                     atgCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCca                                                                                                         <  7:201719/89‑1 (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                >  1:325722/1‑148 (MQ=255)
                                                                                                                                   ccGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGCAAATGCAGCAGCGCCGCGACATTACCGAGCGCTACCTGTCCAAGCTGCCGGCCACCGACCAGa                >  4:292775/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                 >  2:25739/1‑139 (MQ=255)
                                                                                                                                           tcGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGCAAATGCAGCAGCGCCGCGACATTACCGAGCGCTACCTGTCCAAGCTGCCGGCCACCGACCAg                 <  1:25739/139‑1 (MQ=255)
                                                                                                                                                gTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGCAAATGCAGCAGCGCCGCGACATTACCGAGCGCTACCTGTCCAAGCTGCCGGCCACCGACCAGAAAGCCATGAAGGAt  <  3:292775/149‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                    
TGCTCATCCGGTAATCACCGGCCTATGAAAGGAGAGACACCATGCGCAAGACCCTTATCGCCCTGATGTTCGCTGCCGCTTTGCCGACCGTAGCCATGGCCATGCCAGATGGCGGCCCGCGCCACGACGGCCCGCATCATCGCGGTGATGCGCCTTTCCACCAACTGGACCTGAGCCGCGATCAGCGCCAGCAGATCGGCAAGCTGATGGGCGAGCAAATGCAGCAGCGCCGCGACATTACCGAGCGCTACCTGTCCAAGCTGCCGGCCACCGACCAGAAAGCCATGAAGGAT  >  NC_002947/5117711‑5118003

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: