Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947196,3900TC34.3% 77.2 / 59.9 70D632D (GAT→GACPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (23/23);  new base C (10/14);  total (33/37)
Fisher's exact test for biased strand distribution p-value = 6.16e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00
Rejected as consensus: Frequency below/above cutoff threshold.

CCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCAC  >  NC_002947/196243‑196533
                                                                                                                                                   |                                                                                                                                               
ccTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGAtt                                                                                                                                                >  1:18269/1‑149 (MQ=34)
ccTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGAtt                                                                                                                                                >  4:184347/1‑149 (MQ=34)
  tGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTc                                                                                                                                              >  3:170207/1‑149 (MQ=255)
    aCGGCGAGCCGGTCGGTCACCGATGTCGGCGTGATCACTTACACCGTGACCCTGAGCAACCATGCCCAGACGCCGGTGACGGTGACCCTGTCCAACGGCCAAACGATTACAGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTc                                                                                                                                              <  1:286734/147‑1 (MQ=11)
    aCGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTcc                                                                                                                                             <  8:152235/148‑1 (MQ=34)
    aCGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTcc                                                                                                                                             <  3:184347/148‑1 (MQ=34)
    aCGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTcc                                                                                                                                             <  6:349529/148‑1 (MQ=255)
    aCGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTc                                                                                                                                              <  6:141501/147‑1 (MQ=255)
         gAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAACACCCAGGGCAGCGTCGACTTCCAGAccc                                                                                                                                       >  4:215004/1‑149 (MQ=11)
            ccggTCGGTCACCGAAGGAGGCTTTATAACTTACACCGTGAGCGTGTGCAATCCTGCCCAGACTCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCATAACCCAGGGCAGCGTCGATTTCCAGACCCCgg                                                                                                                                    <  8:254246/146‑1 (MQ=11)
            cccGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCgg                                                                                                                                    <  2:206899/149‑1 (MQ=255)
             ccGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGCCGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGCTTACCGTTGAAGCCGGCAACACCCAGGGCAGCGTCGACTTCCAGCCCCCGGc                                                                                                                                   >  7:242799/1‑149 (MQ=11)
             ccGTCGGGCACCTACGGCGGCGTGATCACATACACCGTGACCCGGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGc                                                                                                                                   <  4:235377/149‑1 (MQ=11)
               gTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGa                                                                                                                                 <  2:317302/149‑1 (MQ=255)
                  ggTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAAt                                                                                                                               <  6:160880/148‑1 (MQ=255)
                     cACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGc                                                                                                                                   <  2:315707/141‑1 (MQ=34)
                     cACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGc                                                                                                                                   >  1:315707/1‑141 (MQ=34)
                      actgAAGGCGGCGTTATCAATTACACCGTTACCCTTAGCAACACTGCACAGACGTCGGTGACCTTGACCCTTTCCAACTGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCTGCGTCGACTTCCAGACCCCGGCGAATGACg                                                                                                                           <  7:230711/145‑1 (MQ=11)
                      aCCGAAGGCGGCGTGATCACTTTCACCGTGACCATGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACg                                                                                                                           >  6:256033/1‑148 (MQ=11)
                      aCCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGt                                                                                                                          <  6:153715/149‑1 (MQ=34)
                      aCCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGt                                                                                                                          >  5:220033/1‑149 (MQ=34)
                      aCCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACACTGTCCAACGGCCAAACGATTACCTTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTACAAGACCCCGGCAAATAACg                                                                                                                           >  8:303070/1‑148 (MQ=11)
                       ccGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTc                                                                                                                         <  2:250678/149‑1 (MQ=255)
                       ccGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTc                                                                                                                         >  3:185067/1‑149 (MQ=255)
                        cGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCt                                                                                                                        <  2:223240/149‑1 (MQ=34)
                          aaGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTa                                                                                                                       >  2:40345/1‑148 (MQ=34)
                            ggcggcGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTacaa                                                                                                                    <  7:210796/149‑1 (MQ=34)
                             gcggcgTGATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGTCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTacaa                                                                                                                    >  2:257165/1‑148 (MQ=11)
                               ggcgTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTacaacat                                                                                                                 >  4:181235/1‑148 (MQ=21)
                                    gATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGtt                                                                                                            >  7:53333/1‑149 (MQ=34)
                                    gATCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCGAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGtt                                                                                                            <  6:299982/149‑1 (MQ=11)
                                     aTCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTc                                                                                                           >  6:131547/1‑149 (MQ=34)
                                     aTCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGTCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGtt                                                                                                            <  1:257165/148‑1 (MQ=11)
                                      tCACTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATAACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTc                                                                                                           >  8:261715/1‑148 (MQ=11)
                                        aCTTACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGAc                                                                                                        <  2:115323/149‑1 (MQ=255)
                                         cTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACt                                                                                                       <  2:118236/149‑1 (MQ=34)
                                         cTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACt                                                                                                       <  6:220033/149‑1 (MQ=34)
                                          ttACACCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTg                                                                                                      >  3:285655/1‑149 (MQ=255)
                                          ttACACCGTGACCCAGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATCACGTCTACAACAACGGTTCGACTg                                                                                                      >  1:75333/1‑149 (MQ=21)
                                             cacCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCa                                                                                                   <  4:285655/149‑1 (MQ=255)
                                             cacCGTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGACTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCa                                                                                                   <  6:44820/149‑1 (MQ=255)
                                              acCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAg                                                                                                  >  3:314632/1‑149 (MQ=34)
                                                 gTGACCCTGAGCAACCCTGCCCAGACGCCGGTGACCGTGACCCTGTCCAACGGCCAAACGATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGCCTTACAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAgcgt                                                                                               >  5:128278/1‑149 (MQ=11)
                                                      ccTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCa                                                                                          >  3:165874/1‑149 (MQ=37)
                                                          aGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCg                                                                                       <  5:131547/148‑1 (MQ=37)
                                                             aaTCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGaa                                                                                   <  4:165874/149‑1 (MQ=37)
                                                                    cccAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCAcc                                                                            >  2:307733/1‑149 (MQ=38)
                                                                    cccAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAAACATAACAGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTc                                                                                                                                              >  8:226293/1‑83 (MQ=11)
                                                                        gACGCCGGTGACAGTGACCCTGTCTAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACcggcg                                                                        <  5:64489/149‑1 (MQ=255)
                                                                              ggTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACATCGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAAtt                                                                  <  3:181235/149‑1 (MQ=255)
                                                                                   cAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGTCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGa                                                              <  8:313763/148‑1 (MQ=38)
                                                                                   cAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAg                                                             <  1:307733/149‑1 (MQ=38)
                                                                                   cAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAg                                                             <  2:18269/149‑1 (MQ=38)
                                                                                         cccTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACt                                                        >  3:159200/1‑148 (MQ=37)
                                                                                          ccTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGa                                                      <  4:314632/149‑1 (MQ=37)
                                                                                          ccTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGa                                                      >  4:115023/1‑149 (MQ=37)
                                                                                                  tcGGCCAAACCATTTCCGTTGAAGCCGGCAAATCCCAGGGCAGCGTCGATTTCCAGACGCCGGCGAATGACGTCTACAACAACGGTTGGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAAt                                              <  2:75333/148‑1 (MQ=255)
                                                                                                    ggCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGaa                                               <  7:288891/146‑1 (MQ=37)
                                                                                                    ggCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAAAGGTTCGACTGTCAGTGTAACCATCGCGAACTCCCCCGGCGGTAATTTCGAGCAACTGACCCCGaa                                               >  8:288891/1‑146 (MQ=11)
                                                                                                       cAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGACCAACTGAACCAGAATcag                                           >  8:253304/1‑145 (MQ=255)
                                                                                                        aaaCCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTAAGTGTCACCATAGAGAACGCCACCGGAGGTAATTTCGCGCAACTGAACCCGATTCCGaag                                        >  8:320755/1‑147 (MQ=14)
                                                                                                                 ccGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAg                               <  1:40345/149‑1 (MQ=39)
                                                                                                                        aGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGa                        >  3:128239/1‑149 (MQ=39)
                                                                                                                        aGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGa                        >  8:321754/1‑149 (MQ=39)
                                                                                                                            ggCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCTCGAATCCGACGCCGGCTCAGACCACGATCaa                    <  8:53333/149‑1 (MQ=255)
                                                                                                                            ggCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCaa                    >  6:157446/1‑149 (MQ=39)
                                                                                                                             gtaaaaCCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTAcaacaa                                                                                                                 <  8:112635/53‑1 (MQ=14)
                                                                                                                             gtaaaaCCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTAcaacaa                                                                                                                 >  7:112635/3‑55 (MQ=14)
                                                                                                                                 aaCCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTATCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAg                                                             <  3:194596/103‑1 (MQ=38)
                                                                                                                                 aaCCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTATCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAg                                                             >  4:194596/1‑103 (MQ=38)
                                                                                                                                            gCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGc     >  3:324225/1‑148 (MQ=39)
                                                                                                                                              gTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGccac  >  6:154977/1‑149 (MQ=39)
                                                                                                                                                   |                                                                                                                                               
CCTGACGGCGAGCCCGTCGGTCACCGAAGGCGGCGTGATCACTTACACCGTGACCCTGAGCAATCCTGCCCAGACGCCGGTGACAGTGACCCTGTCCAACGGCCAAACCATTACCGTTGAAGCCGGCAAAACCCAGGGCAGCGTCGATTTCCAGACCCCGGCGAATGACGTCTACAACAACGGTTCGACTGTCAGTGTCACCATCGAGAACGCCACCGGCGGTAATTTCGAGCAACTGACCCCGAATCCGACGCCGGCTCAGACCACGATCAACGACTCGGTCGATGCCAC  >  NC_002947/196243‑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: