Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 202,062 T→A 7.2% L2523Q (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_002947202,0620TA 7.2% 281.3 / 10.5 97L2523Q (CTG→CAG) PP_0168putative surface adhesion protein
Reads supporting (aligned to +/- strand):  ref base T (43/46);  new base A (4/3);  total (48/49)
Fisher's exact test for biased strand distribution p-value = 7.12e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.74e-01

GGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTC  >  NC_002947/201913‑202211
                                                                                                                                                     |                                                                                                                                                     
ggTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCt                                                                                                                                                       <  2:1138956/150‑1 (MQ=255)
  tCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTgg                                                                                                                                                     <  2:1157572/150‑1 (MQ=255)
    aGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTc                                                                                                                                                   <  2:877330/150‑1 (MQ=255)
    aGCACTACCGCCCATGACGAGCCAGGCACTCCAGGCAACGACGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTc                                                                                                                                                   <  1:1364224/150‑1 (MQ=255)
       aCTAACGTCACTGACGAGCCAGGAACCCGAGGCAACGAAGGCGACATCGTCAAGGTCACCATAACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTg                                                                                                                                                <  2:638190/150‑1 (MQ=14)
        cTACCGTCACTGACGAGCCACACACTCCAGGCAACGAAGGCGCCATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCTTGa                                                                                                                                               <  1:900337/150‑1 (MQ=255)
                      gAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCaacaac                                                                                                                                 >  1:1163137/1‑150 (MQ=39)
                         ccAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCTAGAACGTCAAGACAAAATTCACCGTGCACGGCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCt                                                                                                                              >  1:441632/1‑150 (MQ=14)
                            ggCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAg                                                                                                                           <  2:474198/150‑1 (MQ=40)
                               aCTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTc                                                                                                                        >  1:752537/1‑150 (MQ=255)
                                 tCCAGGCAACTAAGGCGACATCCTCAAGGTCACCATCACGGCCGACCAGACTTCATGAGCCCAGAACGTCAAGCCAACGTTCACCGTGCACGTCAACCAGCCGCTTGCTCACGACCTGGTCGTGACCCTGAGTAACAACGCTCAGGTCAc                                                                                                                      <  2:750852/150‑1 (MQ=14)
                                   cAGGTAACGAAGGCGACATCGTCAAGGTCACCATCCCGGCCGACCAGACTTCGGTAGCCGAGAACGTAAAGCCTCCCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTTACCCTGAGCAACAACGATCAGGTCACCa                                                                                                                    >  1:767313/1‑150 (MQ=11)
                                   cAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCa                                                                                                                    >  1:96998/1‑150 (MQ=40)
                                    aGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCAt                                                                                                                   >  1:104379/1‑150 (MQ=40)
                                      gCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCa                                                                                                                 <  2:1163137/150‑1 (MQ=255)
                                       cAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCaa                                                                                                                <  1:88664/150‑1 (MQ=255)
                                          cGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGCCACCATCAAGGc                                                                                                             >  2:187710/1‑150 (MQ=255)
                                              ggCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGt                                                                                                         <  1:1904259/150‑1 (MQ=255)
                                                cGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTga                                                                                                       <  2:927072/150‑1 (MQ=255)
                                                  aCATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTgaga                                                                                                     <  2:31607/150‑1 (MQ=255)
                                                       gTCAAGGTCATTACCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGt                                                                                                <  1:1267527/150‑1 (MQ=255)
                                                       gTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGt                                                                                                <  1:1210274/150‑1 (MQ=255)
                                                         cAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTgc                                                                                              <  2:918652/150‑1 (MQ=255)
                                                          aaGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTgcg                                                                                             >  1:1173266/1‑150 (MQ=255)
                                                                aCCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTac                                                                                       >  1:1274687/1‑150 (MQ=255)
                                                                 ccATCACGGCCGCCCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTaca                                                                                      >  2:79694/1‑150 (MQ=255)
                                                                  cATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTacac                                                                                     <  1:1459760/150‑1 (MQ=255)
                                                                  cATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTacac                                                                                     >  1:635310/1‑150 (MQ=255)
                                                                    tCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCCGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCCAGCAAGCCCGGGTCACCATCAAGGCCCGTGAGCCCAGTGCGCCGTACAccc                                                                                   <  2:273229/150‑1 (MQ=14)
                                                                        ggCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGTTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCacga                                                                               <  2:65396/150‑1 (MQ=255)
                                                                           cGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGAcgc                                                                            <  2:104379/150‑1 (MQ=255)
                                                                             aCCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACgcgc                                                                          >  2:1283894/1‑150 (MQ=255)
                                                                             aCCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACgcgc                                                                          >  2:31175/1‑150 (MQ=255)
                                                                               cAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTTGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCCCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCaa                                                                        >  1:452588/1‑150 (MQ=255)
                                                                                  aCTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCa                                                                                                                            <  2:803059/95‑1 (MQ=255)
                                                                                  aCTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCa                                                                                                                            >  1:803059/1‑95 (MQ=255)
                                                                                  aCTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGc                                                                     >  1:1629889/1‑150 (MQ=255)
                                                                                  aCTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGACCATCAACAAGCAGCTCGATAACGAAATAGTGGCGAGTAGGAACAGCGTGTTTCAGGTGGCGTTGAAGTCCGTTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGc                                                                     >  2:1674853/1‑150 (MQ=9)
                                                                                    ttCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGa                                                                   <  2:1629889/150‑1 (MQ=255)
                                                                                       ggTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTc                                                                                                                             >  1:781944/1‑89 (MQ=40)
                                                                                       ggTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTc                                                                                                                             <  2:781944/89‑1 (MQ=255)
                                                                                        gTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGAc                                                               <  2:752537/150‑1 (MQ=255)
                                                                                                    gTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTgaga                                                                                                     >  1:1678837/1‑100 (MQ=255)
                                                                                                    gTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTgaga                                                                                                     <  2:1678837/100‑1 (MQ=255)
                                                                                                    gTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGAc                                                   >  2:60754/1‑150 (MQ=255)
                                                                                                    gTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGAc                                                   >  2:1856616/1‑150 (MQ=255)
                                                                                                    gTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGAc                                                   >  2:1292736/1‑150 (MQ=255)
                                                                                                     tCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACg                                                  >  1:311126/1‑150 (MQ=255)
                                                                                                       aaGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCt                                                >  1:385090/1‑150 (MQ=255)
                                                                                                        aGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGATCTGGTCGTGAGCCTGAGCAACAACGCTCAGGTCACCATCCAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTg                                               <  2:96998/150‑1 (MQ=255)
                                                                                                         gCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGcaa                                                                                                                                     <  1:79733/63‑1 (MQ=37)
                                                                                                         gCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGcaa                                                                                                                                     >  2:79733/1‑63 (MQ=35)
                                                                                                         gCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTgg                                              >  2:654635/1‑150 (MQ=255)
                                                                                                          ccAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGc                                             <  1:1403466/150‑1 (MQ=255)
                                                                                                            aaCCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACACCGCTCAGGTCACCATCAAGGCCGGTGAGCAAAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGccg                                           >  1:520428/1‑149 (MQ=255)
                                                                                                             aCCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAg                                          >  2:1280095/1‑150 (MQ=255)
                                                                                                               cTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCa                                                       >  2:1089389/1‑135 (MQ=255)
                                                                                                               cTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCa                                                       <  1:1089389/135‑1 (MQ=255)
                                                                                                               cTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGAt                                        >  1:659704/1‑150 (MQ=255)
                                                                                                                ttCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATc                                       >  2:228215/1‑150 (MQ=255)
                                                                                                                ttCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATc                                       <  2:1173266/150‑1 (MQ=255)
                                                                                                                ttCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATc                                       >  1:143494/1‑150 (MQ=255)
                                                                                                                   aCCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAg                                                                                                               >  1:343782/1‑75 (MQ=35)
                                                                                                                   aCCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAg                                                                                                               <  2:343782/75‑1 (MQ=35)
                                                                                                                    ccGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGt                                                                                                >  1:695573/1‑89 (MQ=255)
                                                                                                                    ccGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGt                                                                                                <  2:695573/89‑1 (MQ=255)
                                                                                                                     cGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCt                                  <  1:970514/150‑1 (MQ=255)
                                                                                                                       tGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCAACATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTgg                                <  2:239240/150‑1 (MQ=255)
                                                                                                                          aCGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCa                             >  1:1018452/1‑150 (MQ=255)
                                                                                                                            gTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACAccc                                                                                   <  2:1543578/94‑1 (MQ=255)
                                                                                                                            gTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACAccc                                                                                   >  1:1543578/1‑94 (MQ=255)
                                                                                                                            gTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGc                                                 >  1:109968/1‑128 (MQ=255)
                                                                                                                            gTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGc                                                 <  2:109968/128‑1 (MQ=255)
                                                                                                                            gTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTg                                 >  1:302989/1‑144 (MQ=255)
                                                                                                                            gTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTg                                 <  2:302989/144‑1 (MQ=255)
                                                                                                                            gTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATc                           >  1:1580566/1‑150 (MQ=255)
                                                                                                                              cAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCAc                         <  1:1283894/150‑1 (MQ=255)
                                                                                                                                   aGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCg                    <  2:413390/150‑1 (MQ=255)
                                                                                                                                     ccGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCt                  >  2:64067/1‑150 (MQ=255)
                                                                                                                                        cTGGCTCACGACCTGGTCGTGACCCTGAGCCACCACGCTCAGGTCACCATCAAGGACAGCGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCtgtg               <  1:1289274/150‑1 (MQ=17)
                                                                                                                                        cTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCtgtg               <  1:31175/150‑1 (MQ=255)
                                                                                                                                        cTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCtgtg               <  1:1280095/150‑1 (MQ=255)
                                                                                                                                         tGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTgg              <  1:79694/150‑1 (MQ=255)
                                                                                                                                          gggtCACGACCAGGTCGTGAGCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGa                                         >  1:1326154/4‑122 (MQ=255)
                                                                                                                                          gggtCACGACCAGGTCGTGAGCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGa                                         <  2:1326154/119‑1 (MQ=255)
                                                                                                                                          gggtCACGACCAGGTCGTGAGCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGa             >  1:937817/4‑150 (MQ=255)
                                                                                                                                           gCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGAt            <  1:60754/150‑1 (MQ=255)
                                                                                                                                            cTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAg                                                      >  1:998542/1‑107 (MQ=255)
                                                                                                                                            cTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAg                                                      <  2:998542/107‑1 (MQ=255)
                                                                                                                                             tCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGAt            <  2:1048360/148‑1 (MQ=255)
                                                                                                                                             tCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGAt            >  1:1048360/1‑148 (MQ=255)
                                                                                                                                             tCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGt          >  1:17938/1‑150 (MQ=255)
                                                                                                                                             tCACGACCGGTTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGt          >  2:75501/1‑150 (MQ=255)
                                                                                                                                             tCACGACCAGGTCGTGAGCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGAc                                                   >  1:95120/1‑109 (MQ=255)
                                                                                                                                             tCACGACCAGGTCGTGAGCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGAc                                                   <  2:95120/109‑1 (MQ=255)
                                                                                                                                               aCGACCAGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGAt                                        >  1:112937/1‑118 (MQ=255)
                                                                                                                                               aCGACCAGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGAt                                        <  2:112937/118‑1 (MQ=255)
                                                                                                                                                 gACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGAc      <  2:452588/150‑1 (MQ=255)
                                                                                                                                                     tGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTc  <  1:64067/150‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GGTCAGCACTACCGTCACTGACGAGCCAGGCACTCCAGGCAACGAAGGCGACATCGTCAAGGTCACCATCACGGCCGACCAGACTTCGGTAGCCGAGAACGTCAAGCCAACCTTCACCGTGCACGTCAACCAGCCGCTGGCTCACGACCTGGTCGTGACCCTGAGCAACAACGCTCAGGTCACCATCAAGGCCGGTGAGACCAGTGCGCCGTACACCCACGACGCGCAAGGCGATGACGTCTACCAGGACGCTGGCCAGATCAGCCTGGGCATCACGTCCGCTGTGGATGTAGACGGTC  >  NC_002947/201913‑202211

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