Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 939,928:1 +C 7.7% coding (5697/18933 nt) PP_0806 ← surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947939,9281.C 7.7% 204.1 / 14.9 65coding (5697/18933 nt)PP_0806surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base . (30/30);  new base C (2/3);  total (32/33)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.62e-01

GGAGGCGTTGCCTGCCGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGCCACCACGTCCA  >  NC_002947/939781‑940076
                                                                                                                                                    |                                                                                                                                                     
ggAGGCGTTGCCTGCCGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGGCGGCGGCCAGCACGCGTACg‑gg                                                                                                                                                     >  1:17334/1‑149 (MQ=255)
    gCGTTGCCTGCCGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTATCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCAC‑CGTACGCGtgctg                                                                                                                                                 <  2:538855/150‑1 (MQ=255)
     cGTTGCCTGCCGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCAC‑CGTACGCGtgctgg                                                                                                                                                <  1:170669/150‑2 (MQ=255)
      gTTGCCTGCCGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGcc                                                                                                                                               <  1:533631/150‑1 (MQ=255)
       ttGCCTGCCGCATCGGTGGCCGTGACTTGACGCACCCCGCCCTCTATCGGTGGCGGGTCGAGGGTGATGCTGACTCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGccc                                                                                                                                              <  1:1096841/150‑1 (MQ=17)
          cctgccGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGt                                                                                                                                           >  1:181209/1‑150 (MQ=255)
                   tCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGgcactg                                                                                                                                  >  1:16953/1‑146 (MQ=255)
                      gTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCg                                                                                                                               <  1:723134/150‑1 (MQ=255)
                        ggCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGt                                                                                                                             <  2:447628/150‑1 (MQ=255)
                           cGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGAc                                                                                                                          <  2:94861/150‑1 (MQ=255)
                           cGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGAc                                                                                                                          >  2:488116/1‑150 (MQ=255)
                              gACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCAc                                                                                                                       <  2:1033851/150‑1 (MQ=255)
                                    aaGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGcc                                                                                                                 >  2:249180/1‑150 (MQ=255)
                                         ccTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑Gtgct                                                                                                                                                  >  2:808595/1‑112 (MQ=39)
                                         ccTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑Gtgct                                                                                                                                                  <  1:808595/112‑1 (MQ=38)
                                           tCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCa                                                                                                          >  1:37349/1‑150 (MQ=38)
                                           tCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCa                                                                                                          >  1:835378/1‑150 (MQ=38)
                                                gTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGccc                                                                                                                                              >  1:19297/1‑109 (MQ=33)
                                                gTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGccc                                                                                                                                              <  2:19297/109‑1 (MQ=33)
                                                         tgGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTtc                                                                                           >  1:368644/1‑150 (MQ=35)
                                                           gCGGGTCGAGGGTGATGCTGAAGCTGCCGGTATCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTtcgg                                                                                         >  2:623836/1‑150 (MQ=25)
                                                              ggTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGAcc                                                                                                                         >  2:841367/1‑116 (MQ=37)
                                                              ggTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGAcc                                                                                                                         <  1:841367/116‑1 (MQ=38)
                                                                tCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTcggtcggt                                                                                    <  2:835378/150‑1 (MQ=35)
                                                                              gAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGt                                                                      <  2:114420/150‑1 (MQ=35)
                                                                               aaGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCAGCGTAGGTG‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTg                                                                     >  2:640868/1‑150 (MQ=11)
                                                                                  cTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGata                                                                  >  1:836872/1‑150 (MQ=35)
                                                                                    gCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATc                                                                <  2:836872/150‑1 (MQ=35)
                                                                                      cGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGTGGTGATATCcg                                                              >  1:482688/1‑150 (MQ=25)
                                                                                       ggTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCcgg                                                             <  1:370983/150‑1 (MQ=35)
                                                                                            cgcCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCAt                                                                                                         <  1:187291/102‑1 (MQ=38)
                                                                                            cgcCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCAt                                                                                                         >  2:187291/1‑102 (MQ=38)
                                                                                            cgcCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCcggtgcgg                                                        >  1:892713/1‑150 (MQ=35)
                                                                                             gcCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGCCCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCcggtgcggt                                                       >  2:411999/1‑150 (MQ=25)
                                                                                                aaCTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGTCGGTGAACAC‑CAGGCACCCTTGTAGACCTAGTATGTCGGTCGGTGCGGCAGGGGCGGTGATATCctgtggtgggat                                                    >  1:406841/1‑143 (MQ=2)
                                                                                                aaCTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGAt                                                    >  1:1063299/1‑150 (MQ=34)
                                                                                                aaCTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGAt                                                    <  1:249180/150‑1 (MQ=34)
                                                                                                   tATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGTCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTtg                                                 >  2:588925/1‑150 (MQ=25)
                                                                                                   tATCACGCTACCCAGCACAGTGGTGCCGGCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCCGCAGGGGCGGTGATATCcggtgcggtggctgg                                                 >  2:901676/1‑145 (MQ=255)
                                                                                                    aTCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTtgg                                                <  2:389026/150‑1 (MQ=34)
                                                                                                    aTCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGAATGGGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTtgg                                                >  1:738895/1‑150 (MQ=11)
                                                                                                     tCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATcc                                                               >  2:483363/1‑134 (MQ=33)
                                                                                                     tCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATcc                                                               <  1:483363/134‑1 (MQ=33)
                                                                                                     tCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGa                                                     <  2:786927/144‑1 (MQ=33)
                                                                                                     tCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGa                                                     >  1:786927/1‑144 (MQ=33)
                                                                                                       aCGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGTGGTGATATCCGGTGCGGTGATTtgggtg                                             <  2:482688/150‑1 (MQ=17)
                                                                                                       aCGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTtgggtg                                             >  1:610628/1‑150 (MQ=32)
                                                                                                       aCGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTtgggtg                                             <  1:360123/150‑1 (MQ=32)
                                                                                                        cGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTtgggtgg                                            <  2:181209/150‑1 (MQ=32)
                                                                                                                  cacaGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGa                                  <  2:633294/150‑1 (MQ=32)
                                                                                                                   acaGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGaa                                 >  2:867857/1‑150 (MQ=32)
                                                                                                                    caGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGaaa                                <  2:115614/150‑1 (MQ=32)
                                                                                                                       tggtgCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAacac                             >  1:317701/1‑150 (MQ=32)
                                                                                                                                    gggcAGCAC‑CGTACGCGTGCTGGCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTtgggt                                              <  2:601217/117‑1 (MQ=11)
                                                                                                                                    gggcAGCAC‑CGTACGCGTGCTGGCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTtgggt                                              >  1:601217/4‑120 (MQ=11)
                                                                                                                                    ggcCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAAt                <  1:488116/150‑1 (MQ=32)
                                                                                                                                     ggcAGCAC‑CGTACGCGTGCTGGCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATg               >  1:256848/3‑150 (MQ=11)
                                                                                                                                     gcCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATg               <  1:867857/150‑1 (MQ=32)
                                                                                                                                     cccAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATg               <  1:36137/149‑1 (MQ=17)
                                                                                                                                       cAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGc             <  2:892713/150‑1 (MQ=32)
                                                                                                                                        aGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGcc            >  2:173631/1‑150 (MQ=32)
                                                                                                                                         gCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTTGGGGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGcca           <  1:623836/150‑1 (MQ=17)
                                                                                                                                                tACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGCCACCACGTc    >  2:85724/1‑150 (MQ=32)
                                                                                                                                                  cG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGGCCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGCCACCACGTAca  <  2:16953/150‑1 (MQ=17)
                                                                                                                                                  cG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAAt                >  2:622259/1‑136 (MQ=33)
                                                                                                                                                  cG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAAt                <  1:622259/136‑1 (MQ=33)
                                                                                                                                                    |                                                                                                                                                     
GGAGGCGTTGCCTGCCGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGCCACCACGTCCA  >  NC_002947/939781‑940076

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