Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 196,390 T→C 27.1% D632D (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_002947196,3900TC27.1% 90.4 / 42.0 68D632D (GAT→GACPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (21/27);  new base C (8/10);  total (31/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 6.64e-01

CTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCAC  >  NC_002947/196244‑196533
                                                                                                                                                  |                                                                                                                                               
cTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGAttt                                                                                                                                               <  1:86736/149‑1 (MQ=34)
   aCGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCa                                                                                                                                            <  4:321586/149‑1 (MQ=34)
   aCGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTcc                                                                                                                                             <  2:282860/148‑1 (MQ=255)
     ggCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGa                                                                                                                                          <  5:289354/149‑1 (MQ=34)
     ggCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCATGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGa                                                                                                                                          <  6:276716/149‑1 (MQ=21)
     ggCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGa                                                                                                                                          >  1:31844/1‑149 (MQ=255)
          gCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCg                                                                                                                                     <  5:307103/149‑1 (MQ=255)
           cccGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCgg                                                                                                                                    <  7:65567/149‑1 (MQ=255)
             cGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGc                                                                                                                                   >  2:217845/1‑148 (MQ=34)
               tCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGaa                                                                                                                                <  3:35028/149‑1 (MQ=255)
                 ggTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTCAAGCCGGAAACACCCAAGGCACCGGCGAATTCAAAACCCCCGCCAAt                                                                                                                               >  8:47663/1‑148 (MQ=11)
                    cACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACg                                                                                                                           <  7:73016/149‑1 (MQ=255)
                     aCCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACg                                                                                                                           >  6:217410/1‑148 (MQ=34)
                     aCCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGt                                                                                                                          >  1:206239/1‑149 (MQ=34)
                      ccGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTc                                                                                                                         <  2:31844/149‑1 (MQ=255)
                       cGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGCTTTCCAGACCCCGGCGAATGACGTc                                                                                                                         >  7:111970/1‑148 (MQ=21)
                       cGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTc                                                                                                                         >  3:6931/1‑148 (MQ=255)
                       cGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCAGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTc                                                                                                                         >  2:6359/1‑148 (MQ=11)
                        gAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTa                                                                                                                       <  1:258702/149‑1 (MQ=34)
                         aaGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTac                                                                                                                      <  5:217410/149‑1 (MQ=34)
                                   gATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCATGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCACACCCCCGGGCAGCGTCGAGTTCCAGCCCCCGGCGAATGACGTCGACAACACCGGtg                                                                                                            >  8:251981/1‑148 (MQ=11)
                                     tCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCTGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTc                                                                                                           <  8:177019/148‑1 (MQ=21)
                                     tCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGAACCTGTCCAACGGCCAAACCATTACCGTTGAAGCAGGCAAAACCCAGGGCAGCGACGATATCCAGACCCCGGCGAAAGACGGCTACAACAACGGATCg                                                                                                          >  6:174293/1‑149 (MQ=11)
                                       aCTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGa                                                                                                         <  2:206239/148‑1 (MQ=34)
                                       aCTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGa                                                                                                         <  6:8005/148‑1 (MQ=255)
                                        atctcacCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACt                                                                                                       <  2:233384/145‑1 (MQ=11)
                                         ttACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCAACGGCAGCGTCGATTTCAAGACAACGGAGAATGACGTATACAACAACGGTACGACt                                                                                                       >  1:298586/1‑148 (MQ=11)
                                             acCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAg                                                                                                  <  7:268376/149‑1 (MQ=34)
                                             acCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCGAGTACAACGGCCAAACGATTTCCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTGCCGACCCCGGCGAATGACGTCTACAACAACCGTTCGACt                                                                                                       >  1:233384/1‑144 (MQ=11)
                                             acCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAATCCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGt                                                                                                     <  8:117794/146‑1 (MQ=11)
                                             acCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAGCCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGt                                                                                                     >  7:117794/1‑146 (MQ=11)
                                             acCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCat                                                                                                  >  4:46404/1‑148 (MQ=255)
                                               cGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGAGACCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAgcg                                                                                                >  1:274998/1‑147 (MQ=11)
                                                    cccTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGCGTCAc                                                                                            >  5:194107/1‑148 (MQ=11)
                                                     ccTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGTGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCa                                                                                          <  6:67064/149‑1 (MQ=37)
                                                       tGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATc                                                                                        >  8:184906/1‑149 (MQ=37)
                                                         aGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCga                                                                                      <  4:266767/149‑1 (MQ=37)
                                                          gCAATCCTGCCCAGACTCCTGTGACAGTGACCCTGTCCAACGGCCTAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGGGTCGATTTCCAGACCCCGGCGAATGACGTCTAGAACAACGGTTCGACTGTCAGTGTCACCATCga                                                                                      <  5:174293/148‑1 (MQ=11)
                                                          gCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCTACGGCCAAACCATTACCGTTTAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCga                                                                                      <  7:151342/148‑1 (MQ=255)
                                                           cAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGTTTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCgaga                                                                                    <  8:111970/149‑1 (MQ=37)
                                                             aTCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAAc                                                                                  >  7:199470/1‑149 (MQ=37)
                                                               ccTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGc                                                                                >  5:317885/1‑149 (MQ=37)
                                                                  gCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGATTGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCAc                                                                             >  8:172915/1‑149 (MQ=37)
                                                                     cAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACcgg                                                                          <  6:317885/149‑1 (MQ=39)
                                                                                 aCAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGa                                                              <  7:184906/149‑1 (MQ=39)
                                                                                     tGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGt                                                                      <  2:35451/137‑1 (MQ=12)
                                                                                     tGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCGGGTAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGt                                                                      >  1:35451/1‑137 (MQ=12)
                                                                                     tGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGa                                                              >  3:138138/1‑145 (MQ=37)
                                                                                     tGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGa                                                              <  4:138138/145‑1 (MQ=37)
                                                                                         ccTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCTCCGGCGGTAATTTCGAGCAACTGa                                                      >  3:8938/1‑149 (MQ=37)
                                                                                                   ggCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATcc                                            <  1:217845/149‑1 (MQ=39)
                                                                                                      cAAACCATTACCTTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACTTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGAc                                         <  2:298586/149‑1 (MQ=255)
                                                                                                      cAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCg                                                                                       >  7:168920/1‑103 (MQ=35)
                                                                                                      cAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCg                                                                                       <  8:168920/103‑1 (MQ=35)
                                                                                                               aCCGTTGAAGCCGGCAAAACCCAGGGCACCGGCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCATTGTCACAATCGAGAACGCCACCGGCGGAAATTTCCAGAAACTGACCCCGAAACCGCCACCGGCt                                  >  6:142976/1‑147 (MQ=11)
                                                                                                                ccGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCa                                <  8:199470/148‑1 (MQ=39)
                                                                                                                ccGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTcc                                                                                                                                             <  2:123352/39‑1 (MQ=255)
                                                                                                                   ttGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCa                                >  4:213183/1‑145 (MQ=39)
                                                                                                                   ttGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCa                                <  3:213183/145‑1 (MQ=39)
                                                                                                                   ttGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGAc                             >  5:28417/1‑148 (MQ=39)
                                                                                                                     gAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTAcaacaa                                                                                                                 <  4:267185/62‑1 (MQ=35)
                                                                                                                     gAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTAcaacaa                                                                                                                 >  3:267185/1‑62 (MQ=35)
                                                                                                                      aaGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCAc                          <  7:172915/148‑1 (MQ=39)
                                                                                                                            gCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCTCCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAAc                   <  4:8938/149‑1 (MQ=255)
                                                                                                                                         cAGCGTCGATTTCCAGCCCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATg      <  7:251981/149‑1 (MQ=255)
                                                                                                                                           gCGTCGATTTCCAGACCCCGGCGAATGACGTCAACAACAACGGTTCGACTGTCAGTGTCAGCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTACGACCGCGATCAACCACTAGGTCCATgac    >  4:61025/1‑147 (MQ=14)
                                                                                                                                             gTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGccac  >  4:151685/1‑149 (MQ=39)
                                                                                                                                                 aTTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACg                         <  2:310367/122‑1 (MQ=38)
                                                                                                                                                 aTTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACg                         >  1:310367/1‑122 (MQ=38)
                                                                                                                                                  |                                                                                                                                               
CTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCAC  >  NC_002947/196244‑196533

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: