Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 205,663 C→T 8.0% G3723G (GGC→GGT PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947205,6630CT 8.0% 161.9 / 2.9 61G3723G (GGC→GGTPP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base C (31/25);  new base T (2/3);  total (33/28)
Fisher's exact test for biased strand distribution p-value = 6.53e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.97e-01

TCAGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCG  >  NC_002947/205515‑205810
                                                                                                                                                    |                                                                                                                                                    
tCAGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAggcggc                                                                                                                                                      <  7:166609/149‑1 (MQ=39)
  aGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCga                                                                                                                                                    <  7:305742/149‑1 (MQ=39)
        ggTATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCa                                                                                                                                               <  2:390304/148‑1 (MQ=39)
          tATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTGCCACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCACGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGCGATCACt                                                                                                                                             >  8:168266/1‑148 (MQ=255)
               aCTCGGCGGTGGATGCCAGTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCATGTGACCGACACCACCTACGAAGTGGTGTCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACAcc                                                                                                                                      <  1:236913/149‑1 (MQ=255)
                    gcggTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTaca                                                                                                                                         <  4:278598/142‑1 (MQ=39)
                    gcggTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTaca                                                                                                                                         >  3:278598/1‑142 (MQ=39)
                    gcggTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCAc                                                                                                                                 <  5:284502/149‑1 (MQ=39)
                    gcggTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCAc                                                                                                                                 <  3:305858/149‑1 (MQ=39)
                       gTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCt                                                                                                                              <  2:248297/149‑1 (MQ=39)
                               cACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCTGCGAGATCACTTACACCATCACGCTGACcaac                                                                                                                       >  8:177600/1‑148 (MQ=21)
                                 cTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACcaaca                                                                                                                      >  4:143301/1‑147 (MQ=34)
                                  tGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACcaaca                                                                                                                      >  4:243485/1‑146 (MQ=34)
                                       cGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGt                                                                                                              <  6:277912/149‑1 (MQ=34)
                                        gACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTc                                                                                                             <  2:258898/149‑1 (MQ=34)
                                           cTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGc                                                                                                          >  2:222855/1‑149 (MQ=34)
                                               gagaACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGa                                                                                                       <  2:213263/148‑1 (MQ=35)
                                                agaACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGa                                                                                                                  <  7:138893/136‑1 (MQ=34)
                                                agaACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGa                                                                                                                  >  8:138893/1‑136 (MQ=34)
                                                agaACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGAt                                                                                                      >  6:113884/1‑148 (MQ=37)
                                                 gaACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACAACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATc                                                                                                     >  1:84613/1‑148 (MQ=37)
                                                              ggcggcAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGAc                                                                                                                           >  4:386857/1‑113 (MQ=35)
                                                              ggcggcAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGAc                                                                                                                           <  3:386857/113‑1 (MQ=35)
                                                                 ggcAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCAtt                                                                                    <  7:177600/149‑1 (MQ=37)
                                                                              ttCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCt                                                                        >  2:89322/1‑148 (MQ=37)
                                                                                  aGTGACCGACACAACCGACGAAGTGGTCGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTCCCGAt                                                                                                      >  7:302996/1‑114 (MQ=255)
                                                                                   gTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGa                                                                  >  3:48880/1‑149 (MQ=37)
                                                                                      aCCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACgg                                                               >  2:245704/1‑149 (MQ=34)
                                                                                              caccGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCg                                                       <  1:222855/149‑1 (MQ=34)
                                                                                              caccGACGAAGTGGTGGCCAAGCTGACCGAGACCCCGTCGGTCACCGAAGCCGGCGAGATCACTTACACCATCACGATGACCAACAAAGACGGTCTGCCGAACAACAACCACAGTCCATTGACATTCACGCTGACCGACGGCAAAAcc                                                        >  7:164360/1‑148 (MQ=11)
                                                                                               accGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGt                                                      >  7:67268/1‑149 (MQ=34)
                                                                                                tcaccgaaggggtggcggtcagtttggccaccactTCGGTGACCGAAGGCGGTGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGt                                                      <  7:38310‑M1/113‑1 (MQ=255)
                                                                                                ccGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGtc                                                     <  5:113884/149‑1 (MQ=34)
                                                                                                ccGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTAGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGtc                                                     <  3:7465/149‑1 (MQ=21)
                                                                                                   ccgaaggggtggcggtcagtttggccaccactTCGGTGAACGAAGGCGGTGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGtcat                                                   <  1:7953‑M1/116‑1 (MQ=255)
                                                                                                   acgaAGTGGTGGCCAAGCTGACCGCGACCCAGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTCACCTTCACGCTGAGCGACGGCAAAACCGtcat                                                   >  3:247570/1‑148 (MQ=21)
                                                                                                     cacgggggggcggtcagttgggccaccactTCGGGGACCGAAGGCGGTGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCAc                                                >  5:168585‑M1/31‑149 (MQ=255)
                                                                                                      aaggggtggcggtcagtttggccaccactTCGGTGACCGAAGGCGGTGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGAt                                                                                                      <  7:250450‑M1/65‑1 (MQ=255)
                                                                                                      aaggggtggcggtcagtttggccaccactTCGGTGACCGAAGGCGGTGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGAt                                                                                                      >  8:250450‑M1/30‑94 (MQ=255)
                                                                                                        gtggtgGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGt                                             <  3:188997/149‑1 (MQ=34)
                                                                                                         tggtggCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCg                                                                   <  1:63119/126‑1 (MQ=35)
                                                                                                         tggtggCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCg                                                                   >  2:63119/1‑126 (MQ=35)
                                                                                                           gtggCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGcc                                          >  4:158816/1‑149 (MQ=34)
                                                                                                              gCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGc                                       <  1:89322/149‑1 (MQ=34)
                                                                                                              gCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGc                                       >  2:32739/1‑149 (MQ=34)
                                                                                                                cAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACCACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGcc                                      >  3:174418/1‑148 (MQ=21)
                                                                                                                cAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATAACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGcc                                      >  6:123694/1‑148 (MQ=21)
                                                                                                                  aGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACGGTGCCgttggcc                                   <  8:67268/149‑7 (MQ=255)
                                                                                                                   gctgcacgcgACCCCGTCGGTCACCGAAGTCGGCGAGATAACTTACACCATCACGCTGACCAACAAAGACTGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACg                                  <  4:48880/143‑1 (MQ=11)
                                                                                                                    cTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGTATTTACCTTCACGCTGa                                                                      >  8:199358/1‑112 (MQ=255)
                                                                                                                    cTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGTATTGACCTTCACGCTGa                                                                      <  7:199358/112‑1 (MQ=25)
                                                                                                                        ccGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCAcc                             >  6:40104/1‑149 (MQ=35)
                                                                                                                            gACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTgg                         <  2:84613/149‑1 (MQ=34)
                                                                                                                             aCCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACAGCCATCCCGCTGACCAACAAAGACGGTCTGCCGATCAACACCCACAGTGCATTGCCGTTCACGATGAGCGCCGGCGAAACCGTCATAACCGTGCCGGCCAACGGCCCCGTgg                         >  1:71708/1‑149 (MQ=11)
                                                                                                                               cccGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCa                      >  5:233400/1‑149 (MQ=34)
                                                                                                                                ccGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCAc                     <  1:32739/149‑1 (MQ=34)
                                                                                                                                 cGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACt                    >  3:22568/1‑149 (MQ=34)
                                                                                                                                   TCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCGAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGc                  >  2:71708‑M1/1‑149 (MQ=255)
                                                                                                                                      gTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTAc               <  3:143301/149‑1 (MQ=34)
                                                                                                                                      gTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTAc               >  2:232347/1‑149 (MQ=34)
                                                                                                                                        cACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCg             >  2:233205/1‑149 (MQ=34)
                                                                                                                                           cGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCa          >  3:400196/1‑149 (MQ=34)
                                                                                                                                               ggcggcGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGc      <  6:233400/149‑1 (MQ=34)
                                                                                                                                                gcggcgAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATAAACAACCCCAGTGCATTGACATTCACGCTCAGCGAAGGCAAAACCGTCATCACCGTGCCGGCCAACGCCACCGTGGGCACTGATACCGACACTGca     >  5:90838/1‑148 (MQ=11)
                                                                                                                                                   gcgAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGAGCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTTGGCACTGCTACCGTCACTGCCCCg  <  4:247570/149‑1 (MQ=21)
                                                                                                                                                    |                                                                                                                                                    
TCAGCCTGGGTATCAACTCGGCGGTGGATGCCACTGGTGCGACCTTCGAGAACCTGCAACTGGGCGGCAACGCTTCGGTTCAAGTGACCGACACCACCGACGAAGTGGTGGCCAAGCTGACCGCGACCCCGTCGGTCACCGAAGGCGGCGAGATCACTTACACCATCACGCTGACCAACAAAGACGGTCTGCCGATCAACAACCACAGTGCATTGACCTTCACGCTGAGCGACGGCAAAACCGTCATCACCGTGCCGGCCAACGGCACCGTGGGCACTGCTACCGTCACTGCCCCG  >  NC_002947/205515‑205810

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: