Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 197,590 T→C 15.3% D1032D (GAT→GAC PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947197,5900TC15.3% 204.3 / 26.5 85D1032D (GAT→GACPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (32/40);  new base C (5/8);  total (37/48)
Fisher's exact test for biased strand distribution p-value = 7.68e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.62e-01

CCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCAC  >  NC_002947/197442‑197733
                                                                                                                                                    |                                                                                                                                               
cccTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGAt                                                                                                                                                 <  5:143128/149‑1 (MQ=37)
   tGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTc                                                                                                                                              >  4:196023/1‑149 (MQ=37)
     aCGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTcc                                                                                                                                             <  6:209664/148‑1 (MQ=37)
     aCGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTc                                                                                                                                              <  4:143264/147‑1 (MQ=37)
     aCGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCa                                                                                                                                            <  4:210645/149‑1 (MQ=37)
          gAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGAccc                                                                                                                                       <  3:196023/149‑1 (MQ=37)
          gAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGAccc                                                                                                                                       <  6:242578/149‑1 (MQ=37)
             cccGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCgg                                                                                                                                    <  2:8633/149‑1 (MQ=37)
              ccGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCgg                                                                                                                                    <  4:199330/148‑1 (MQ=37)
              ccGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGc                                                                                                                                   <  6:138544/149‑1 (MQ=37)
                 tCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGaa                                                                                                                                <  6:69507/149‑1 (MQ=37)
                        ccGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTc                                                                                                                         >  7:231770/1‑149 (MQ=37)
                        ccGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGCCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGt                                                                                                                          >  4:16473/1‑148 (MQ=37)
                         cGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAAt                                                                                                                               <  5:61625/142‑1 (MQ=37)
                         cGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAAt                                                                                                                               >  6:61625/1‑142 (MQ=37)
                         cGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCt                                                                                                                        <  3:231656/149‑1 (MQ=37)
                                ggcgTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTAcaacaa                                                                                                                 <  2:23769/149‑1 (MQ=37)
                                  cgTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACg                                                                                                               <  1:264973/149‑1 (MQ=37)
                                    tGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGt                                                                                                             <  8:70878/149‑1 (MQ=37)
                                         aCCTACACCGTGATCGTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTTCCGTTTAAGCCGGCAAGACCCATGGCAGCGTCGATTTCCAGACCCCGGCTAATGACGTCTACAACAACGGTTCGa                                                                                                         <  7:125383/148‑1 (MQ=11)
                                          ccTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACt                                                                                                       >  2:210257/1‑149 (MQ=38)
                                          ccTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACt                                                                                                       <  3:16473/149‑1 (MQ=38)
                                           cTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACt                                                                                                       >  8:240225/1‑148 (MQ=38)
                                           cTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCAGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCgcccg                                                                                                      >  5:39308/1‑145 (MQ=255)
                                              cacCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGt                                                                                                     <  5:236392/147‑1 (MQ=39)
                                              cacCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCa                                                                                                   <  3:154868/149‑1 (MQ=39)
                                                ccGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGAc                                                                                                                            <  6:190011/122‑1 (MQ=38)
                                                ccGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGAc                                                                                                                            >  5:190011/1‑122 (MQ=38)
                                                     aTCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCa                                                                                             <  4:182162/148‑1 (MQ=18)
                                                      tCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCAc                                                                                            <  3:179073/148‑1 (MQ=38)
                                                      tCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCAcc                                                                                           <  8:231770/149‑1 (MQ=17)
                                                        cTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACg                                                                                                                           <  8:131843/115‑1 (MQ=35)
                                                        cTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACg                                                                                                                           >  7:131843/1‑115 (MQ=35)
                                                         tGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTc                                                                                                           <  7:226620/130‑1 (MQ=35)
                                                         tGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTc                                                                                                           >  8:226620/1‑130 (MQ=35)
                                                         tGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATc                                                                                        >  3:172738/1‑149 (MQ=37)
                                                         tGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCAt                                                                                         >  8:263321/1‑148 (MQ=17)
                                                            gaaaTCCTGCCCAGCCGCCGGTAACCGTTACCCTGTACAACGGCCAAACCATTACCGTTGAAGCCGGCTAGACCCATGGCAGCGTCGATTTCCAGACCCCTGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCga                                                                                      <  3:230944/146‑1 (MQ=11)
                                                             cAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCgaga                                                                                    <  1:210257/149‑1 (MQ=37)
                                                                  cTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGc                                                                                >  7:209175/1‑148 (MQ=17)
                                                                  cTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGc                                                                                >  5:145906/1‑148 (MQ=17)
                                                                   tGCCCAGACGCCGGAAACCGTGACCCTGTCCAACGGCCAAACCATTACCTTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACt                                                                                                       <  1:259307/124‑1 (MQ=255)
                                                                   tGCCCAGACGCCGGAAACCGTGACACTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACt                                                                                                       >  2:259307/1‑124 (MQ=255)
                                                                     cccAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACCGTCAGCGTCACCATCGAGAACGCcac                                                                             <  7:263321/148‑1 (MQ=17)
                                                                      ccAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAg                                                                           <  4:258575/149‑1 (MQ=9)
                                                                        aGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGc                                                                         <  6:82809/149‑1 (MQ=255)
                                                                                   aCCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACAGTTCGACTGTCAGCGTCGCCATCGAGAACGCCACAGGCGGTAATTTc                                                                >  7:190501/1‑147 (MQ=21)
                                                                                     cGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAg                                                             >  2:84007/1‑148 (MQ=12)
                                                                                     cGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAg                                                             >  6:11502/1‑148 (MQ=32)
                                                                                         aCCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACt                                                        <  4:58205/149‑1 (MQ=12)
                                                                                                  aaCGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGaa                                               >  1:123229/1‑149 (MQ=32)
                                                                                                    cGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCgg                                                                       >  6:77379/1‑123 (MQ=32)
                                                                                                    cGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGATCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCgg                                                                       <  5:77379/123‑1 (MQ=17)
                                                                                                       ccAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGc                                                            >  1:166828/1‑131 (MQ=12)
                                                                                                       ccAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGc                                                            <  2:166828/131‑1 (MQ=12)
                                                                                                        cAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGAc                                         >  2:305247/1‑149 (MQ=34)
                                                                                                         aaaCGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGa                                                              >  4:231317/1‑127 (MQ=12)
                                                                                                         aaaCGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGa                                                              <  3:231317/127‑1 (MQ=12)
                                                                                                         aaaCCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACg                                        >  5:268994/1‑149 (MQ=35)
                                                                                                         aaaCCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACg                                        >  7:74551/1‑149 (MQ=35)
                                                                                                           acgaTTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGAcccc                                      <  1:84007/146‑1 (MQ=255)
                                                                                                             cATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCgg                                    >  3:21556/1‑149 (MQ=37)
                                                                                                              aTTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGc                                   >  2:105409/1‑149 (MQ=37)
                                                                                                                 aCCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTCACCCCGAACCCGACGCCGGCt                                  >  2:25577/1‑147 (MQ=37)
                                                                                                                  ccGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCa                                <  8:74551/148‑1 (MQ=37)
                                                                                                                  ccGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAg                               <  3:285289/149‑1 (MQ=37)
                                                                                                                  ccGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAg                               <  5:11502/149‑1 (MQ=37)
                                                                                                                   cGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGa                              >  8:239612/1‑149 (MQ=37)
                                                                                                                   cGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGa                              >  7:98354/1‑149 (MQ=12)
                                                                                                                    gTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGAc                             >  3:270019/1‑149 (MQ=37)
                                                                                                                     ttGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGACACAGGCGGTAATTTCGAGGACCTGACCCCGAACAGGACCCCGGCTCAGAc                             >  6:188616/1‑148 (MQ=11)
                                                                                                                        aaGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCAc                          >  6:598/1‑148 (MQ=37)
                                                                                                                        aaGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACAGTTCGACTGTCAGCGTCGCCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCAc                          <  8:190501/148‑1 (MQ=255)
                                                                                                                          gCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGAt                       <  6:268994/149‑1 (MQ=37)
                                                                                                                              gCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCAcc                   <  7:239612/149‑1 (MQ=38)
                                                                                                                              gCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATCaac                   <  5:188616/149‑3 (MQ=2)
                                                                                                                               cAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACACCGAACCCGACGCCGGCTCAGACCACGATCAcc                   >  1:123624/1‑148 (MQ=38)
                                                                                                                                    cccAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCgg                                    <  8:23832/126‑1 (MQ=38)
                                                                                                                                     ccAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGa                                                              <  5:206344/99‑1 (MQ=33)
                                                                                                                                     ccAGGGCAGCGTCGATTTCCAGACCCCGGCCAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGa                                                              >  6:206344/1‑99 (MQ=18)
                                                                                                                                       aGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTc          <  7:240225/149‑1 (MQ=39)
                                                                                                                                       aGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTc          <  4:270019/149‑1 (MQ=39)
                                                                                                                                              cGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCg             >  5:252269/1‑139 (MQ=40)
                                                                                                                                              cGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCg             <  6:252269/139‑1 (MQ=40)
                                                                                                                                               gTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACCCCGGCTCAGACCACGATc                      <  2:7957/129‑1 (MQ=17)
                                                                                                                                                tCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGAcgccac  >  2:142580/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
CCCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACCTACACCGTGATCCTGAGCAATCCTGCCCAGACGCCGGTAACCGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAGACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCACCATCGAGAACGCCACAGGCGGTAATTTCGAGCAACTGACCCCGAACCCGACGCCGGCTCAGACCACGATCACCGACTCGGTCGACACCAC  >  NC_002947/197442‑197733

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: