Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 2,099,995 G→C 6.9% intergenic (‑63/‑142) PP_1875 ← / → PP_1876 two‑component system sensor histidine kinase/response regulator/AAA family ATP/GTP‑binding protein

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029472,099,9950GC 6.9% 212.0 / 4.0 68intergenic (‑63/‑142)PP_1875/PP_1876two‑component system sensor histidine kinase/response regulator/AAA family ATP/GTP‑binding protein
Reads supporting (aligned to +/- strand):  ref base G (32/31);  new base C (3/2);  total (35/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 = 1.00e+00

TTGTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTC  >  NC_002947/2099856‑2100133
                                                                                                                                           |                                                                                                                                           
ttGTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCG                                                                                                                                   >  2:1329554‑M1/1‑150 (MQ=255)
 tGTATGACAGGCGATTTGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGC                                                                                                                                  <  1:533745‑M1/150‑1 (MQ=255)
 tGTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGC                                                                                                                                  <  1:533543‑M1/150‑1 (MQ=255)
         aGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGC                                                                                                                          <  2:698838‑M1/150‑1 (MQ=255)
          ggCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGTGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCG                                                                                                                         <  1:1359400‑M1/150‑1 (MQ=255)
             gATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGC                                                                                                                      >  2:190364‑M1/1‑150 (MQ=255)
                          gATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGTGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGGCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCca                                                                                                         >  1:686332/1‑150 (MQ=255)
                           atatCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCcac                                                                                                        >  2:1117796/1‑150 (MQ=255)
                            tctccATGATTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGGTTGCCCTgtcccacc                                                                                                       >  2:1421561‑M1/3‑142 (MQ=255)
                              tCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTGCGCGCCCCTCCGCGGGCTTGCCCGctc                                                                                                            <  1:191217‑M1/143‑4 (MQ=255)
                              tCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTGCGCGCCCCTCCGCGGGCTTGCCCGctc                                                                                                            >  2:191217‑M1/1‑140 (MQ=255)
                               ccATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGt                                                                                                    >  2:737841/1‑150 (MQ=255)
                                 aTGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCCAGCGGGGGCCGCTGCGCGCCCCTCCGCGGGCTTGCCCGctcccacagggct                                                                                                  >  2:749495‑M1/1‑137 (MQ=255)
                                 aTGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTAc                                                                                                  >  2:1271205/1‑150 (MQ=255)
                                        ccGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAg                                                                                           <  1:1407579/150‑1 (MQ=255)
                                         cGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGc                                                                                          >  2:112787/1‑150 (MQ=255)
                                         cGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCAGGTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGc                                                                                          <  2:1183440/150‑1 (MQ=255)
                                                 tttCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCACCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCCGGCt                                                                                 >  2:1027578/1‑150 (MQ=255)
                                                   tCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggct                                                                                >  2:786521‑M1/1‑119 (MQ=255)
                                                     ccGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAAc                                                                             <  1:112787/150‑1 (MQ=255)
                                                      cGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaag                                                                             >  1:270429‑M1/1‑116 (MQ=255)
                                                        ttCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGCCCCCGCTTGGCGa                                                                                                                                  <  2:449868‑M1/95‑2 (MQ=255)
                                                        ttCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGCACCCGCTTGGCGa                                                                                                                                  >  1:449868‑M1/1‑94 (MQ=255)
                                                         tCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCaa                                                                         <  2:943148/150‑1 (MQ=255)
                                                         tCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCaa                                                                         <  1:190364/150‑1 (MQ=255)
                                                            ttGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGcc                                                                      <  2:266196/150‑1 (MQ=255)
                                                                gctgcgcAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGaa                                                                              <  1:479701/138‑1 (MQ=255)
                                                                gctgcgcAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGaa                                                                              >  2:479701/1‑138 (MQ=255)
                                                                  tgcgcAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgg                                                                 <  1:786521‑M1/150‑47 (MQ=255)
                                                                  tgcgcAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCtgtg                                                                 >  2:359702/1‑149 (MQ=255)
                                                                  tgcgcAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCtgtg                                                                 <  1:359702/149‑1 (MQ=255)
                                                                  tgcgcAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTgg                                                                <  1:1271205/150‑1 (MQ=255)
                                                                  tgcgcAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGGCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTgg                                                                <  2:686332/150‑1 (MQ=255)
                                                                    ggcACGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGa                                                              <  2:631941/149‑1 (MQ=255)
                                                                     gcAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTGCGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcattcgggggctaagctgtacctgtgggaa                                                              <  1:749495‑M1/150‑50 (MQ=255)
                                                                               ccGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGCGcgg                                                          <  2:230122‑M2/143‑4 (MQ=255)
                                                                               ccGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgg                                                          >  1:230122‑M2/1‑140 (MQ=255)
                                                                                  aaTAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgggcaagcccgc                                                >  1:1117796‑M2/1‑137 (MQ=255)
                                                                                   aTAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGCCCCCGCTTGGC                                                                                                                                    <  1:619648‑M1/66‑1 (MQ=255)
                                                                                   aTAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGCCCCCGCTTGGC                                                                                                                                    >  2:619648‑M1/1‑66 (MQ=255)
                                                                                    tAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCCGTTGATCCCCCAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgggcaagcccgcgg                                              <  1:959326‑M2/150‑16 (MQ=255)
                                                                                     aGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaacgggcaagcccgcgga                                              <  2:270429‑M1/150‑66 (MQ=255)
                                                                                     aGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTCTGCGGCTAACCGTTTGCCGCTTGGCGCCCCTCCGCTTGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgggcaagcccgcgga                                             >  1:1027578‑M2/1‑134 (MQ=255)
                                                                                       ccccGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgggcaagcccgcggagg                                           >  1:913923‑M2/1‑132 (MQ=255)
                                                                                         ccGGCTGCTGGTGGGTCTGTCTTTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTa                                                                                            <  1:318311/100‑1 (MQ=255)
                                                                                         ccGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaacggg                                                          <  2:1071867‑M1/134‑54 (MQ=255)
                                                                                         ccGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTa                                                                                            >  2:318311/1‑100 (MQ=255)
                                                                                         ccGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgggcaagcccgcggagggg                                         <  1:7367‑M2/150‑21 (MQ=255)
                                                                                         ccGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCACTCCGCGGGCTTGCCCGctcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaacggg                                                          >  1:1071867‑M1/1‑81 (MQ=255)
                                                                                          cGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGCCTGCTTGGGGCCGTTCAGCTGGGGTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgggcaagcccgcggaggggc                                        >  1:737841‑M2/1‑129 (MQ=255)
                                                                                                 tggtggGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGc                                 >  1:1211794/1‑150 (MQ=255)
                                                                                                   gtggGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCCTGCCCGTTCCCCCAGGTCCTGCTTCGCCCCGTTGGTGACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCgg                               >  2:1261929/1‑150 (MQ=37)
                                                                                                     gggTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGcc                             <  1:634776/150‑1 (MQ=255)
                                                                                                           gtctgtTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCaaa                       >  2:1308991/1‑150 (MQ=255)
                                                                                                              tgtTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTGCAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGGAAGCCCGCGGTGTGGCGCGCAGCGGCCCCCaaaaaa                    <  1:1261929/150‑1 (MQ=255)
                                                                                                              tgtTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTAAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCaaaaaa                    >  1:666888/1‑150 (MQ=255)
                                                                                                               gtTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgg                                                          >  1:743414‑M2/1‑108 (MQ=255)
                                                                                                               gtTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgg                                                          <  2:743414‑M2/111‑4 (MQ=255)
                                                                                                                ttCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCaaaaaaat                  >  1:691681/1‑148 (MQ=255)
                                                                                                                    tttCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTgaga              >  1:1010409/1‑150 (MQ=255)
                                                                                                                        aGTTGATCGCCAAGTGGGGGCCGCTGCGCGCCCCTCGGCGGGCTTGCCCGcgcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaacgggcaagcccgcggag                                             >  2:713907‑M1/1‑50 (MQ=255)
                                                                                                                        aGTTGATCGCCAAGTGGGGGCCGCTGCGCGCCCCTCCGCGGGCTTGCCCGcgcccacagggcttgcgttcagcctggggctaagctgtacctgtgggaacgggcaagcccgcggag                                             <  1:713907‑M1/116‑67 (MQ=255)
                                                                                                                        aGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGt          <  2:1211794/150‑1 (MQ=255)
                                                                                                                              tCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGCGcgggcaagcccgcggagggg                                         <  2:902768‑M2/113‑21 (MQ=255)
                                                                                                                              tCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgggcaagcccgcggagggg                                         >  1:902768‑M2/1‑93 (MQ=255)
                                                                                                                                gCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTGCCAAAGGTGCAGGTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGGCCGCGGGGGGGCGGGCGGCGGCCGCGAAAAAACTGAGCTTGTTGGGAgtc  >  1:413717/1‑150 (MQ=255)
                                                                                                                                gCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgggcaagcccgcggagg                                           >  2:1186342‑M2/1‑91 (MQ=255)
                                                                                                                                gCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgggcaagcccgcggagg                                           <  1:1186342‑M2/109‑19 (MQ=255)
                                                                                                                                        ggggCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgg                                                          >  1:284076‑M2/1‑83 (MQ=255)
                                                                                                                                        ggggCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGcgg                                                          <  2:284076‑M2/86‑4 (MQ=255)
                                                                                                                                        ggggCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTga                <  2:1105456/128‑1 (MQ=255)
                                                                                                                                        ggggCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTga                >  1:1105456/1‑128 (MQ=255)
                                                                                                                                           gCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCtgt                                                                  <  2:553201/75‑1 (MQ=255)
                                                                                                                                           gCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCtgt                                                                  >  1:553201/1‑75 (MQ=255)
                                                                                                                                           |                                                                                                                                           
TTGTATGACAGGCGATTGGTGAACTTGATATCCATGAGTCCCGAGCCGGTTTCCCGTTCCTTGCGCTGCGCAAGCATAGCCGAATAGCCCCGGCTGCTGGTGGGTCTGTCTGTTCATTTCAGTTGATCGCCAAGCGGGGGCCGCTTGGCGCCCCTCCGCGGGCTTGCCCGTTCCCACAGGTACAGCTTAGCCCCAGGCTGAACGCAAGCCCTGTGGGAGCGGGCAAGCCCGCGGAGGGGCGCGCAGCGGCCCCCAAAAAACTGAGATTGTTAGGAGTC  >  NC_002947/2099856‑2100133

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: