Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 214,014 T→A 9.7% L6507Q (CTG→CAG)  PP_0168 → putative surface adhesion protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_002947214,0140TA 9.7% 209.2 / 12.7 72L6507Q (CTG→CAG) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (29/36);  new base A (3/4);  total (32/40)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.21e-01

GGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAATTCGGCGGTAGACGCCACCGGC  >  NC_002947/213865‑214162
                                                                                                                                                     |                                                                                                                                                    
ggTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCt                                                                                                                                                      <  2:989416/150‑1 (MQ=255)
     gCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGCCCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCg                                                                                                                                                 <  2:1001133/150‑1 (MQ=255)
       aCTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCAGGTCGTg                                                                                                                                               >  1:332383/1‑150 (MQ=255)
        cTTCGGTCGCCGACGAGCCGGGTAGCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGa                                                                                                                                              <  2:1590660/150‑1 (MQ=255)
        cTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGTCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCAAGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGa                                                                                                                                              <  1:55066/150‑1 (MQ=255)
        cTTCGGTCACCGAAGCACCGGGTACCACACCGAACCAAGACTACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCAGGTCGTgg                                                                                                                                              <  2:1468835/150‑2 (MQ=14)
         ttCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGAc                                                                                                                                             >  2:514639/1‑150 (MQ=255)
         ttCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCAGGTCGTggg                                                                                                                                             <  2:332383/150‑3 (MQ=255)
              tCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGa                                                                                                                                        <  1:2074356/150‑1 (MQ=255)
               cACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAg                                                                                                                                       <  1:987460/150‑1 (MQ=255)
                aCCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGc                                                                                                                                      <  1:1280596/150‑1 (MQ=255)
                  cgacgaGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGcaa                                                                                                                                    >  2:1816175/1‑150 (MQ=255)
                  cgacgaGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGcaa                                                                                                                                    <  1:514639/150‑1 (MQ=255)
                   gacgaGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGcaac                                                                                                                                   <  1:80962/150‑1 (MQ=255)
                         ccGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGcc                                                                                                                             <  1:1064054/150‑1 (MQ=255)
                           aggTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCa                                                                                                                           >  2:1651025/2‑150 (MQ=255)
                            ggTACCCCAGGCAACGAAGGCGCCCCGCTCAATGTCACCATCACCTCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAg                                                                                                                          <  1:1134510/150‑1 (MQ=255)
                            ggTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAg                                                                                                                          <  2:597502/150‑1 (MQ=255)
                            ggTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAg                                                                                                                          <  1:2304768/150‑1 (MQ=255)
                                   cAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCAGGTCGTgg                                                                                                                                              >  1:1705181/1‑122 (MQ=255)
                                   cAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCAGGTCGTgg                                                                                                                                              <  2:1705181/123‑2 (MQ=255)
                                       cAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACTTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCaa                                                                                                               <  2:1449334/150‑1 (MQ=255)
                                            aaGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTCGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCg                                                                                                          <  1:916149/150‑1 (MQ=255)
                                              ggCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGc                                                                                                        <  1:1816175/150‑1 (MQ=255)
                                                 gACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGGCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGaa                                                                                                     >  1:1031195/1‑150 (MQ=38)
                                                  aCCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGaaa                                                                                                    <  1:1435442/150‑1 (MQ=38)
                                                       gTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAgc                                                                                               >  2:383205/1‑150 (MQ=35)
                                                        tCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAg                                                                                                                                       <  1:2357680/109‑1 (MQ=38)
                                                        tCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAg                                                                                                                                       >  2:2357680/1‑109 (MQ=38)
                                                          aaGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGa                                                                                                                                              <  1:1093933/100‑1 (MQ=38)
                                                          aaGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGa                                                                                                                                              >  2:1093933/1‑100 (MQ=39)
                                                          aaGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAgcgcg                                                                                            >  2:297367/1‑150 (MQ=35)
                                                          aaGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCAGGTCGTgg                                                                                                                                              <  2:1580477/100‑2 (MQ=255)
                                                          aaGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCAGGTCGTgg                                                                                                                                              >  1:1580477/1‑99 (MQ=255)
                                                                 ccATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGCCCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCAGCATCAAGGCCGGCGAAACCAGCGCGCCGTaca                                                                                     <  1:973702/150‑1 (MQ=255)
                                                                  cATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGCCCTGGTCGTGACCCTGAGCAACAACGccc                                                                                                                            <  1:2245679/110‑1 (MQ=255)
                                                                  cATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGccc                                                                                                                            >  2:2245679/1‑110 (MQ=38)
                                                                     cACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCa                                                                                 <  2:257499/150‑1 (MQ=38)
                                                                        cgccgACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCacga                                                                              >  1:2071957/1‑150 (MQ=40)
                                                                         gccgACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAg                                                                                                                          <  2:2360112/105‑1 (MQ=38)
                                                                         gccgACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAg                                                                                                                          >  1:2360112/1‑105 (MQ=38)
                                                                                aGACCTCGGTGGCAGAGGACGTCAAACCGACCTCCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAg                                                                      <  1:2057225/150‑1 (MQ=255)
                                                                                  aCCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGc                                                                    >  2:2353660/1‑150 (MQ=255)
                                                                                    cTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGa                                                                  >  1:1574678/1‑150 (MQ=40)
                                                                                     tCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGAt                                                                 <  1:1651025/150‑1 (MQ=39)
                                                                                     tCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGAt                                                                 <  1:297367/150‑1 (MQ=38)
                                                                                           gCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTc                                                           <  1:2383233/150‑1 (MQ=38)
                                                                                             cGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCgg                                                                                                         >  1:1654229/1‑102 (MQ=38)
                                                                                             cGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCgc                                                                                                         <  2:1654229/102‑2 (MQ=38)
                                                                                                        aaCCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTg                                              <  1:1381757/150‑1 (MQ=255)
                                                                                                              ccTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGACaaa                                        >  2:667511/1‑148 (MQ=255)
                                                                                                               cTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGAt                                       >  2:364847/1‑150 (MQ=255)
                                                                                                                ttCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCt                                                          <  1:1258890/130‑1 (MQ=40)
                                                                                                                ttCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCt                                                          >  2:1258890/1‑130 (MQ=40)
                                                                                                                ttCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATc                                      >  2:417779/1‑150 (MQ=255)
                                                                                                                    cTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGcc                                  >  1:738988/1‑150 (MQ=255)
                                                                                                                       tgCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTgg                               <  1:667511/150‑1 (MQ=40)
                                                                                                                       tgCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTgg                               <  1:2353660/150‑1 (MQ=255)
                                                                                                                            gTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCa                                                                                 >  2:2005992/1‑95 (MQ=34)
                                                                                                                            gTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCa                                                                                 <  1:2005992/95‑1 (MQ=33)
                                                                                                                                        cTGGCCCACGACCTGGTCGTGACCCTGAGCGACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCac                                                                                >  1:604355/1‑84 (MQ=14)
                                                                                                                                        cTGGCCCACGACCTGGTCGTGACCCTGAGCGACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCac                                                                                <  2:604355/84‑1 (MQ=14)
                                                                                                                                          ggCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAATTCGGCGGTAGa            >  1:2129511/1‑150 (MQ=255)
                                                                                                                                             ccACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCAt                           >  1:1697383/1‑132 (MQ=40)
                                                                                                                                             ccACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCAt                           <  2:1697383/132‑1 (MQ=40)
                                                                                                                                             ccACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAATTCGGCGGTAGACGc         >  2:815694/1‑150 (MQ=255)
                                                                                                                                               aCGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCa                                     >  1:620944/1‑120 (MQ=40)
                                                                                                                                               aCGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCa                                     <  2:620944/120‑1 (MQ=40)
                                                                                                                                                 gACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAATTCGGCGGTAGACGCCAcc     >  2:483798/1‑150 (MQ=255)
                                                                                                                                                 gACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAATTCGGCGGTAGACGCCAcc     >  1:659971/1‑150 (MQ=255)
                                                                                                                                                    cTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAATTCGGCGGTAGACGCCACCGgc  <  1:483798/150‑1 (MQ=255)
                                                                                                                                                    cTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAATTCGGCGGTAGACGCCACCGgc  >  2:425210/1‑150 (MQ=255)
                                                                                                                                                    cTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAATTCGGCGGTAGACGCCACCGgc  >  2:810577/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                    
GGTCAGCACTTCGGTCACCGACGAGCCGGGTACCCCAGGCAACGAAGGCGACCCGGTCAAGGTCACCATCACCGCCGACCAGACCTCGGTGGCCGAGAACGTCAAACCGACCTTCACTGTGCACGTCAACCAGCCGCTGGCCCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGCGAAACCAGCGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACAACGATGCTGGCCAGATCAGCCTGGGCATCAATTCGGCGGTAGACGCCACCGGC  >  NC_002947/213865‑214162

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