Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,605,513 C→G 5.9% A310G (GCC→GGC)  PP_4077 → hypothetical protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,605,5130CG 5.9% 216.3 / 4.2 68A310G (GCC→GGC) PP_4077hypothetical protein
Reads supporting (aligned to +/- strand):  ref base C (30/33);  new base G (2/2);  total (32/36)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.51e-01

CCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTCGAATCGTTCGCCCACGA  >  NC_002947/4605368‑4605660
                                                                                                                                                 |                                                                                                                                                   
ccACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGcc                                                                                                                                                    >  7:73969/1‑147 (MQ=255)
    aGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTc                                                                                                                                              <  5:121169/149‑1 (MQ=255)
      gTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGa                                                                                                                                            <  8:357462/149‑1 (MQ=255)
      gTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGa                                                                                                                                            <  8:351833/149‑1 (MQ=255)
       tGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGGCAGCCTCGCTGGCCAGGTCGat                                                                                                                                           <  3:392944/149‑2 (MQ=37)
         ggCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCt                                                                                                                                         <  7:162571/149‑1 (MQ=255)
            aTCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGaa                                                                                                                                      <  2:197238/149‑1 (MQ=255)
            aTCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTACACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCCGCCCAGGTCGAGCTg                                                                                                                                        >  2:200062/1‑147 (MQ=255)
                    ccATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAg                                                                                                                                >  1:329017/1‑147 (MQ=255)
                      aTTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGc                                                                                                                            >  2:212263/1‑149 (MQ=255)
                       tttCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGc                                                                                                                            >  5:45885/1‑148 (MQ=255)
                       tttCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGc                                                                                                                            >  5:411064/1‑148 (MQ=255)
                           aCGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGTCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAgg                                                                                                                        <  8:24931/148‑1 (MQ=255)
                           aCGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGgc                                                                                                                       <  6:45885/149‑1 (MQ=255)
                              cgTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCa                                                                                                                    <  1:325764/149‑1 (MQ=255)
                                tCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGt                                                                                                                  >  8:227316/1‑149 (MQ=255)
                                   aCTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCg                                                                                                               <  7:227316/149‑1 (MQ=255)
                                       ccAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCtt                                                                                                            <  3:181357/148‑1 (MQ=255)
                                        cAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTAc                                                                                                          <  1:359833/149‑1 (MQ=255)
                                             acCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGgcg                                                                                                                      >  4:200308/1‑132 (MQ=255)
                                             acCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGgcg                                                                                                                      <  3:200308/132‑1 (MQ=255)
                                               ccTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCa                                                                                                                                                  >  6:306486/1‑102 (MQ=255)
                                               ccTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCa                                                                                                                                                  <  5:306486/102‑1 (MQ=255)
                                                cTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCTAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCtt                                                                                                  >  3:33015/1‑149 (MQ=255)
                                                 ttCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCtt                                                                                                  >  8:7935/1‑148 (MQ=255)
                                                   cgcgTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCa                                                                                                <  2:329017/148‑1 (MQ=255)
                                                      gTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCCAGCTGAAGCCCAGAAGCGAGACGCAGTCCGCTTACGAAATCTTCAgtg                                                                                             >  4:207729/1‑148 (MQ=255)
                                                                aagaTTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGt                                                                                   >  6:279765/1‑148 (MQ=255)
                                                                    tttCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGtcatc                                                                               <  5:279765/148‑1 (MQ=255)
                                                                             gTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGa                                                                      >  2:298859/1‑148 (MQ=255)
                                                                                 tGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGgcgc                                                                                                                     <  5:356373/97‑1 (MQ=255)
                                                                                 tGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGgcgc                                                                                                                     >  6:356373/1‑97 (MQ=255)
                                                                                  gCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGaacaac                                                                 >  5:189002/1‑148 (MQ=255)
                                                                                   cACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGAt                                                                                        >  2:369876/1‑124 (MQ=255)
                                                                                   cACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGAt                                                                                        <  1:369876/124‑1 (MQ=255)
                                                                                   cACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGaacaac                                                                 >  3:353360/1‑147 (MQ=255)
                                                                                   cACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGaacaac                                                                 <  4:353360/147‑1 (MQ=255)
                                                                                   cACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGaacaa                                                                  <  6:125302/146‑1 (MQ=255)
                                                                                   cACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCa                                                               <  6:411064/149‑1 (MQ=255)
                                                                                    aCCCCGTCTGTCATCTTGTTGAATCACAGTGCCGATCCCATCGACCTGGCCATCGTTGCTGACCATGTCGAGCTGAAGCCCAGAAGCGATGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAAccc                                                               <  3:207729/148‑2 (MQ=255)
                                                                                      cccGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGaa                                                                                     >  6:300861/1‑124 (MQ=255)
                                                                                      cccGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACt                                                             >  6:185505/1‑148 (MQ=255)
                                                                                      cccGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACt                                                             >  3:403644/1‑148 (MQ=255)
                                                                                              tCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAg                                                     <  4:403644/148‑1 (MQ=255)
                                                                                              tCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCGGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGt                                                    <  1:200062/149‑1 (MQ=255)
                                                                                                aaCCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCtt                                                                                                  <  1:197129/101‑1 (MQ=255)
                                                                                                aaCCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCtt                                                                                                  >  2:197129/1‑101 (MQ=255)
                                                                                                aaCCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTAc                                                  >  8:265601/1‑149 (MQ=255)
                                                                                                       tCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCTAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGa                                            <  4:33015/148‑1 (MQ=255)
                                                                                                          aGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTATCGGTGAGCAc                                        >  2:186455/1‑149 (MQ=255)
                                                                                                              caGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGc                                    <  7:7935/149‑1 (MQ=255)
                                                                                                                   ccGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGa                                 <  2:249525/147‑1 (MQ=255)
                                                                                                                   ccGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGa                                 >  1:249525/1‑147 (MQ=255)
                                                                                                                   ccGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGAc                                <  1:212263/148‑1 (MQ=255)
                                                                                                                   ccGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACt                               >  6:331592/1‑149 (MQ=255)
                                                                                                                          cATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGcc                        <  1:298859/149‑1 (MQ=255)
                                                                                                                           aTCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGccc                       >  4:35358/1‑149 (MQ=255)
                                                                                                                            tCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACATGCGGACTTTTCGCCCt                      >  5:344991/1‑149 (MQ=255)
                                                                                                                              gACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCttt                    <  8:73969/149‑1 (MQ=255)
                                                                                                                               aCCTGGGCAGCCTCGCTGGCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGtcatc                                                                               <  6:24922/89‑1 (MQ=37)
                                                                                                                               aCCTGGGCAGCCTCGCTGGCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGtcatc                                                                               >  5:24922/1‑89 (MQ=37)
                                                                                                                               aCCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTc                   <  6:189002/149‑1 (MQ=255)
                                                                                                                               aCCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTc                   >  2:323977/1‑149 (MQ=255)
                                                                                                                                ccTGGGCAGCCTCGCTGGCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTc                   >  4:28314/1‑148 (MQ=38)
                                                                                                                                    gCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTCGAAt               <  5:185505/148‑1 (MQ=255)
                                                                                                                                      cAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACtt                              <  8:227750/131‑1 (MQ=255)
                                                                                                                                      cAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTCGAATCGt            >  6:394271/1‑149 (MQ=255)
                                                                                                                                       aGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTCGAATCGtt           >  2:7158/1‑149 (MQ=255)
                                                                                                                                        gCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTCGAATCGTTc          <  7:265601/149‑1 (MQ=255)
                                                                                                                                                 cccAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTAAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTCGAATCTTTCGCCCACGa  >  8:165096/1‑148 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
CCACAGGTGGGCATCGAGTTCCATTTCACGCGTCAACTGCCAGACACCCTTCGCGTTGGCACAGAAGATTTCAGCCTGTTCTGCACCCCGGCTGTCAACCTGTTCAAGCACAGTGCCGAGCCCATCGACCTGGCCAGCGAGGCTGCCCAGGTCGAGCTGAAGCCCAGAAGCGAGGCGCAGTCCGCTTACGAAATCTTCAGTGTCGATGAAGTCATCAGCACCCGAACAACCACTGACGGAAGTACCGGTGAGCACCTGCGGACTTTTCGCCCTTTCGAATCGTTCGCCCACGA  >  NC_002947/4605368‑4605660

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