Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 939,922 Δ1 bp 7.9% coding (5703/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,9220G. 7.9% 176.6 / 11.5 63coding (5703/18933 nt)PP_0806surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base G (30/28);  new base . (2/3);  total (32/31)
Fisher's exact test for biased strand distribution p-value = 6.72e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.40e-01

TGGCGGTGGAGGCGTTGCCTGCCGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGCCA  >  NC_002947/939774‑940067
                                                                                                                                                    |                                                                                                                                                   
tGGCGGTGGAGGCGTTGCCTGCCGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCccgcccgc                                                                                                                                                    >  2:439548/1‑150 (MQ=255)
    ggTGGAGGCGTTGCCTGCCGCATCGATGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTAc                                                                                                                                                >  2:370983/1‑150 (MQ=255)
       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                                                                                                                                         <  1:808595/112‑1 (MQ=38)
                                                ccTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑Gtgct                                                                                                                                         >  2:808595/1‑112 (MQ=39)
                                                  tCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCa                                                                                                 >  1:37349/1‑150 (MQ=38)
                                                  tCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCa                                                                                                 >  1:835378/1‑150 (MQ=38)
                                                       gTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGccc                                                                                                                                     <  2:19297/109‑1 (MQ=33)
                                                       gTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGccc                                                                                                                                     >  1:19297/1‑109 (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:249180/150‑1 (MQ=34)
                                                                                                       aaCTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGAt                                           >  1:1063299/1‑150 (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                                            >  1:786927/1‑144 (MQ=33)
                                                                                                            tCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGa                                            <  2:786927/144‑1 (MQ=33)
                                                                                                              aCGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGTGGTGATATCCGGTGCGGTGATTtgggtg                                    <  2:482688/150‑1 (MQ=17)
                                                                                                              aCGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTtgggtg                                    <  1:360123/150‑1 (MQ=32)
                                                                                                              aCGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTtgggtg                                    >  1:610628/1‑150 (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                                     >  1:601217/4‑120 (MQ=11)
                                                                                                                                           gggcAGCAC‑CGTACGCGTGCTGGCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTtgggt                                     <  2:601217/117‑1 (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)
                                                                                                                                                    |                                                                                                                                                   
TGGCGGTGGAGGCGTTGCCTGCCGCATCGGTGGCCGTGACTTGAAGCACCTCGCCGTCTATCTGTGGCGGGTCGAGGGTGATGCTGAAGCTGCCGGTAGCGCCAACTATCACGCTACCCAGCACAGTGGTGCCGTCGGCGGCCAGCACGCGTACG‑GTGCTGCCCGGTTCGGCGCGGCCGGTGACCACGCTGCCGTCGGCATTGATC‑ACCAGTTCGGTCGGTGCGGCAGGGGCGGTGATATCCGGTGCGGTGATTTGGGTGGGCAGCGAGGAAACACCGCCGTCATCAATGGCCA  >  NC_002947/939774‑940067

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: ‑