Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 217,045 T→A 6.5% A7517A (GCT→GCA PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947217,0450TA 6.5% 210.1 / 2.7 71A7517A (GCT→GCAPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (31/34);  new base A (3/2);  total (35/36)
Fisher's exact test for biased strand distribution p-value = 6.69e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

TAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTGACACTGGACAAGACCCCGGTCAGCACCGCCGTGACCGACGAAC  >  NC_002947/216900‑217190
                                                                                                                                                  |                                                                                                                                                  
tAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTAcc                                                                                                                                                 <  3:249461/149‑1 (MQ=35)
   gatttgatcacagggtcgttggtaccgacgtagacgttgtccggggcagtgacggtagcagtgcCCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGt                                                                                                                                               >  7:195026‑M1/65‑148 (MQ=255)
      aaGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACt                                                                                                                                           >  6:319352/1‑149 (MQ=34)
          cggcgAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGccc                                                                                                                                       <  3:96907/149‑1 (MQ=34)
              gagaTCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGa                                                                                                                                   <  7:11288/149‑1 (MQ=34)
                  tCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAAc                                                                                                                               <  6:39519/149‑1 (MQ=34)
                    aCCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGt                                                                                                                             <  5:8518/149‑1 (MQ=34)
                     ccTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGt                                                                                                                             >  6:218919/1‑148 (MQ=34)
                      cTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTc                                                                                                                            >  4:277320/1‑148 (MQ=34)
                      cTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCt                                                                                                                           <  8:61357/149‑1 (MQ=34)
                      cTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCt                                                                                                                           <  5:319352/149‑1 (MQ=34)
                           ccATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTc                                                                                                                      >  2:20152/1‑149 (MQ=34)
                            cATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCg                                                                                                                     <  2:23884/149‑1 (MQ=34)
                             aTCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCgg                                                                                                                    <  3:231856/149‑1 (MQ=34)
                                    tGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAAc                                                                                                             >  6:283439/1‑149 (MQ=34)
                                         aacaaAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGAccc                                                                                                        <  4:173228/149‑1 (MQ=34)
                                          acaaAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCa                                                                                                               <  2:71443/141‑1 (MQ=34)
                                          acaaAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCa                                                                                                               >  1:71443/1‑141 (MQ=34)
                                            aaaGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCtgt                                                                                                     >  4:134709/1‑149 (MQ=34)
                                                aTGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATc                                                                                                 <  5:268173/149‑1 (MQ=33)
                                                    tCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAat                                                                                            >  6:86304/1‑149 (MQ=32)
                                                    tCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAat                                                                                            <  1:20152/149‑1 (MQ=32)
                                                       acggtagcagtgcccacggtgccgttggccggcacggtgatgacggttttgccgtcgctcagcgtgaaggtcaatgcactgtggGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAatc                                                                                           >  7:195026‑M2/85‑148 (MQ=255)
                                                          aaTCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCTGGACAACGTCTACGTCGGTAc                                                                                                                 <  3:177796/123‑1 (MQ=11)
                                                          aaTCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCTGGACAACGTCTACGTCGGTAc                                                                                                                 >  4:177796/1‑123 (MQ=11)
                                                                  gcccacggtgccgttggccggcacggtgatgacggttttgccgtcgctcagcgtgaaggtcaatgcactgtggGCACTGCCACCGTCACTGCCCCGGACAATGTCTACGTCGGTACCAACGACCCTGTCGTGATGTCGATCGCAACCGt                                                                               >  3:242584‑M2/74‑149 (MQ=255)
                                                                   cccacggtgccgttggcaggcacggtgatgacggttttgccgtcgctcagcgtgaaggtcaatgcactgtggGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGtt                                                                              >  1:230268‑M2/73‑149 (MQ=255)
                                                                   CCACAGTGCATTGACCTTCACGCTGAGCGAGCGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCg                                                                                                                                     <  1:28030‑M1/95‑1 (MQ=255)
                                                                   CCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCg                                                                                                                                     >  2:28030‑M1/1‑94 (MQ=255)
                                                                   CCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGtt                                                                              <  6:33313‑M1/148‑1 (MQ=255)
                                                                   CCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTg                                                                             <  2:216814‑M1/149‑1 (MQ=255)
                                                                   CCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTg                                                                             >  5:21000‑M1/1‑149 (MQ=255)
                                                                    ccacggtgccgttggccggcacggtgatgacggttttgccgtcgctcagcgtgaaggtcaatgcactgtggGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGGTCGCAACCGTTg                                                                             >  8:282247‑M2/72‑149 (MQ=255)
                                                                    CACAGTGCATTGAACTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCGAACGGCTCCGTGGGCACTGCTACCGTCACTGCCCCAGACAACTTCTACGTCGGTAGCAACGACCCTGTGATCAAATCGATCGCAACCGTTgg                                                                            >  4:129148‑M1/1‑149 (MQ=255)
                                                                     cacggtgccgttggccggcacggtgatgacggttttgccgtcgctcagcgtgaaggtcaatgcactgtggGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTg                                                                             >  8:148202‑M2/71‑148 (MQ=255)
                                                                     cacggtgccgttggccggcacggtgatgacggttttgccgtcgctcagcgtgaaggtcaatgcactgtggGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCCCAACCGTTg                                                                             >  5:68932‑M2/71‑148 (MQ=255)
                                                                      acgttgccgttggccggcacggtgatgacggttttgcagtcgatcagcgttaaggtcaatgcactgtggTCACTGCTACCGTCACTGCCCCGGACAAGGTCTACGTCGGTACCAACGACCCTGTGATCAAatcgat                                                                                        <  3:195062‑M2/67‑1 (MQ=255)
                                                                      acggtgccgttggccggcacggtgatgacggttttgccgtcgctcagcgtgaaggtcaatgcactgtggGCACTGCTACCGTCACTGACCCGGACAACGTCGACGTCGGTACCAACGACCCTGTGATCAAatcgat                                                                                        >  4:195062‑M2/70‑136 (MQ=255)
                                                                      CAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCg                                                                          >  5:283439‑M1/1‑149 (MQ=255)
                                                                        ggtgccgttggccggcacggtgatgacggttctgccgtcgctcagcgtgaaggtcaatgcactgtggGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGCACCAACGACCCTGTCGTGATGTCGATCGCAACCGTTggagg                                                                         >  4:220513‑M2/68‑149 (MQ=255)
                                                                        gTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATAGATCGCAACCGTTggcggc                                                                        >  2:46485/1‑149 (MQ=18)
                                                                           cATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGgcgc                                                                      <  5:86304/148‑1 (MQ=35)
                                                                               ttggccggcacggtgatgacggttttgccgtcgctcagcgtgaaggtcaatgcactgtggGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGAt                                                                  <  8:195026‑M2/89‑1 (MQ=255)
                                                                               gACCTTCACGCTGAGCGACTGCAAAACCGTCATCACCGTGCCGGCCTACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGAt                                                                  <  7:276483/148‑1 (MQ=255)
                                                                               gACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATg                                                                 >  6:105549/1‑149 (MQ=38)
                                                                                aCCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGt                                                                <  5:235609/149‑1 (MQ=39)
                                                                                aCCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGt                                                                >  6:303017/1‑149 (MQ=39)
                                                                                       cacggtgatgacggttttgccgtcgctcagcgtgaaggtcaatgcactgtggGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGATGTCGGCaa                                                          <  7:148202‑M2/97‑1 (MQ=255)
                                                                                       cGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCATCGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATTCGATCGCAACCGTTGGCGGCGCGGATGTTGGCaa                                                          <  1:46485/149‑1 (MQ=255)
                                                                                         cTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGtt                                                       <  5:303017/149‑1 (MQ=255)
                                                                                          tGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTc                                                      >  1:149843/1‑149 (MQ=255)
                                                                                            aGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTc                                                      <  4:209891/147‑1 (MQ=255)
                                                                                               tgacggttttgccgtcgctcagcgtgaaggtcaatgcactgtggGCACTGCAACCGTGACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGAAGGTGCGGATGTCGGCAAGTTCGAAc                                                  <  2:230268‑M2/105‑1 (MQ=255)
                                                                                                aCGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGtt                                                               <  3:45323/134‑1 (MQ=40)
                                                                                                aCGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGtt                                                               >  4:45323/1‑134 (MQ=40)
                                                                                                  cGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGc                                               <  1:160159/149‑1 (MQ=255)
                                                                                                   ggCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGa                                                    >  3:294414/1‑143 (MQ=255)
                                                                                                   ggCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGa                                                    <  4:294414/143‑1 (MQ=255)
                                                                                                     cAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTGa                                            >  2:102033/1‑149 (MQ=255)
                                                                                                        aaCCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTGacac                                         >  8:139922/1‑149 (MQ=255)
                                                                                                             tcatcaCCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTGACACTGGAc                                    <  2:149843/149‑1 (MQ=255)
                                                                                                                          ggCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTGACCCTGGACAAGACCCCGCTcc                       >  4:52886/1‑148 (MQ=255)
                                                                                                                           gCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTGACACTGGACAAGAccc                             >  3:269284/1‑142 (MQ=255)
                                                                                                                           gCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTGACACTGGACAAGAccc                             <  4:269284/142‑1 (MQ=255)
                                                                                                                             cAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGaa                                                   >  3:194340/1‑118 (MQ=255)
                                                                                                                             cAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGaa                                                   <  4:194340/118‑1 (MQ=255)
                                                                                                                              aaCGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTgg                                                             >  3:275516/1‑107 (MQ=255)
                                                                                                                              aaCGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTgg                                                             <  4:275516/107‑1 (MQ=255)
                                                                                                                                     ccGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTGACACTGGACAAGCCCCCGGTCAGCACCGCCGTg            <  7:139922/149‑1 (MQ=255)
                                                                                                                                               tGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTGACACTGGACAAGACCCCGGTCAGCACCGCCGTCACCGACGAAc  >  8:109178/1‑149 (MQ=255)
                                                                                                                                                 cTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTGACACTGGACAAGACCCCGGTCAGCACCGCCGTGACCGACGAAc  <  3:277320/147‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                  
TAACCGAAGGCGGCGAGATCACCTACACCATCACCCTGACCAACAAAGATGGTCTGCCAATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCGGACAACGTCTACGTCGGTACCAACGACCCTGTGATCAAATCGATCGCAACCGTTGGCGGCGCGGATGTTGGCAAGTTCGAACAGCTGACACTGGACAAGACCCCGGTCAGCACCGCCGTGACCGACGAAC  >  NC_002947/216900‑217190

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: