Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 4,360,333 G→C 8.3% A278G (GCC→GGC)  PP_3835 ← two‑component system sensor histidine kinase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029474,360,3330GC 8.3% 403.2 / 24.4 133A278G (GCC→GGC) PP_3835two‑component system sensor histidine kinase
Reads supporting (aligned to +/- strand):  ref base G (63/58);  new base C (6/5);  total (69/64)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.45e-01

Only 100 of 134 total aligned reads displayed.

GCTCGAAGATGCGCGGCAGTATCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTCAAGCAGGCCACGCGTAGCGTTAAC  >  NC_002947/4360184‑4360481
                                                                                                                                                     |                                                                                                                                                     
gCTCGAAGATGCGCGGCAGTATCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCGTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTTATGCGGCCATGCCCGTCGATgg                                                                                                                                                       <  1:1144051/150‑1 (MQ=255)
  tCGAAGATGCGCGGCAGTATCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCa                                                                                                                                                     <  1:872291/150‑1 (MQ=255)
  tCGAAGATGCGCGGCAGTATCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTGATGCCAt                                                                                                                                                    >  1:1974180/1‑150 (MQ=255)
    gAAGATGCGCGGCAGTATCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGGTGGCATc                                                                                                                                                   >  2:287377/1‑150 (MQ=255)
    gAAGATGCGCGGCAGTATCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATc                                                                                                                                                   >  1:2737356/1‑150 (MQ=255)
    gAAGATGCGCGGCAGTATCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATc                                                                                                                                                   >  1:2139978/1‑150 (MQ=255)
    gAAGATGCGCGGCAGTATCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATc                                                                                                                                                   >  2:1222876/1‑150 (MQ=255)
          gcgcGGCAGTATCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTGATGCCATCGACGGGc                                                                                                                                            <  2:1974180/150‑1 (MQ=255)
            gcgGCAGTATCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCAt                                                                                                                                          <  2:2681257/150‑1 (MQ=255)
             cgGCAGTATCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCAtt                                                                                                                                         <  1:92733/150‑1 (MQ=255)
             cgGCAGTATCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCAtt                                                                                                                                         <  2:2051830/150‑1 (MQ=255)
                    aTCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGaacaac                                                                                                                                  >  1:841539/1‑150 (MQ=255)
                     tctcTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGGTGGCATCACGGGCATTGAacaaca                                                                                                                                 <  1:287377/150‑1 (MQ=255)
                     tctcTGTGGGTATGCCTTCGCCGCTGTCTTCAACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAacaaca                                                                                                                                 >  2:1431798/1‑150 (MQ=255)
                       tctGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCa                                                                                                                                           <  1:1295878/138‑1 (MQ=255)
                       tctGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCa                                                                                                                                           >  2:1295878/1‑138 (MQ=255)
                       tctGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAgg                                                                                                                               >  2:131460/1‑150 (MQ=255)
                         tgtgGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGCTGGCATCACGGGCATTGAACAACAGGtt                                                                                                                             >  1:2955043/1‑150 (MQ=255)
                         tgtgGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGtt                                                                                                                             >  2:365527/1‑150 (MQ=255)
                           tgGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGa                                                                                                                           <  2:802579/150‑1 (MQ=255)
                           tgGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGa                                                                                                                           <  1:365527/150‑1 (MQ=255)
                           tgGGTATGCCTTCGCCGCTGTCTTCAACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGa                                                                                                                           <  1:1431798/150‑1 (MQ=255)
                             ggTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGc                                                                                                                         >  2:2503609/1‑150 (MQ=255)
                             ggTATGCCTGCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGc                                                                                                                         <  1:131460/150‑1 (MQ=255)
                               tATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTGATGCCATCGACggg                                                                                                                                             >  2:2509696/1‑128 (MQ=255)
                               tATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTGATGCCATCGACggg                                                                                                                                             <  1:2509696/128‑1 (MQ=255)
                                   ccTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACAcc                                                                                                                   >  1:578588/1‑150 (MQ=255)
                                   ccTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACAcc                                                                                                                   >  2:1131119/1‑150 (MQ=255)
                                     ttCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCAt                                                                                                                 >  1:152810/1‑150 (MQ=255)
                                          cgcTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCa                                                                                                            <  1:2503609/150‑1 (MQ=255)
                                             tGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAac                                                                                                                      <  1:645095/137‑1 (MQ=255)
                                                 ttCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAg                                                                                                                          <  1:2787299/129‑1 (MQ=255)
                                                 ttCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAg                                                                                                                          >  2:2787299/1‑129 (MQ=255)
                                                 ttCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGccc                                                                                                     >  2:2847847/1‑150 (MQ=255)
                                                     aCCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCa                                                                                                 >  2:532646/1‑150 (MQ=255)
                                                        tCGAAGCGCACGACGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATcc                                                                                              >  1:612331/1‑150 (MQ=255)
                                                             gcgcACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCAACTGCAGTTGCCCGGCATCCACCTc                                                                                         <  2:1832990/150‑1 (MQ=255)
                                                               gcACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCa                                                                                                            >  2:1983161/1‑129 (MQ=255)
                                                                 aCGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATc                                                                                     >  1:310020/1‑150 (MQ=255)
                                                                 aCGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATc                                                                                     >  2:2991171/1‑150 (MQ=255)
                                                                   gAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTGATGCCATCGACGGGCAtgg                                                                                                                                        <  2:1701527/97‑3 (MQ=38)
                                                                    aaGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAAc                                                                                  <  2:612331/150‑1 (MQ=255)
                                                                    aaGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAAc                                                                                  <  1:2991171/150‑1 (MQ=255)
                                                                     aGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACg                                                                                 >  2:1583757/1‑150 (MQ=255)
                                                                       cTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGTCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGc                                                                               >  1:773330/1‑150 (MQ=255)
                                                                            cttcagcccttcaACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTc                                                                          >  1:894898/1‑150 (MQ=255)
                                                                             ttcagcccttcaACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGa                                                                                                                                       >  1:2233645/1‑88 (MQ=255)
                                                                                agcccttcaACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGtt                                                                                                                             <  1:1001084/95‑1 (MQ=255)
                                                                                agcccttcaACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGtt                                                                                                                             >  2:1001084/1‑95 (MQ=255)
                                                                                    cttcaACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGACATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTgca                                                                  >  1:1658612/1‑150 (MQ=255)
                                                                                           ttagCCGGACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGCCTATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGCCCGGTTGCAGCAacac                                                           >  2:1375722/4‑150 (MQ=255)
                                                                                                 gAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCaa                                                     <  1:1222876/150‑1 (MQ=255)
                                                                                                  aacaacTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTgca                                                                  <  1:1394647/136‑1 (MQ=255)
                                                                                                   acaacTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGt                                                   <  2:152810/150‑1 (MQ=255)
                                                                                                       cTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCt                                                         >  2:1717298/1‑140 (MQ=255)
                                                                                                       cTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCt                                                         <  1:1717298/140‑1 (MQ=255)
                                                                                                        tCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAg                                                                             <  2:112498/119‑1 (MQ=255)
                                                                                                        tCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAg                                                                             >  1:112498/1‑119 (MQ=255)
                                                                                                             cacaTCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATc                                              >  2:971235/1‑145 (MQ=255)
                                                                                                             cacaTCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCt                                         >  2:466537/1‑150 (MQ=255)
                                                                                                               catcgacgggcatggccgcatcagcctgactgtcgatgGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTg                                       >  1:2252104‑M2/39‑150 (MQ=255)
                                                                                                               CATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTg                                       >  1:2143761‑M1/1‑150 (MQ=255)
                                                                                                                ATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGt                                                                      <  1:182880‑M1/118‑1 (MQ=255)
                                                                                                                ATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGt                                                                      >  2:182880‑M1/1‑118 (MQ=255)
                                                                                                                 TCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACTGGTATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGc                                     <  2:2227485‑M1/150‑1 (MQ=255)
                                                                                                                 TCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACAcc                                                                                                                   <  2:408247‑M1/72‑1 (MQ=255)
                                                                                                                 TCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACAcc                                                                                                                   >  1:408247‑M1/1‑72 (MQ=255)
                                                                                                                   GACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTggcga                                   <  1:2803441‑M1/150‑3 (MQ=255)
                                                                                                                      aGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGc                                <  2:619203/150‑1 (MQ=255)
                                                                                                                        tCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCa                                                                                                            <  1:1790703/72‑1 (MQ=255)
                                                                                                                        tCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCCGTTGTGTGGCATCCCCCGGCAGGGATGTCATGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCAt                              <  1:1910846/150‑1 (MQ=255)
                                                                                                                        tCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCACTCGAGCTTggcatggcat                              >  2:2521711/1‑150 (MQ=255)
                                                                                                                        tCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATc                                              <  2:1266923/134‑1 (MQ=255)
                                                                                                                        tCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCAt                              <  2:452653/150‑1 (MQ=255)
                                                                                                                        tCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCAt                              <  1:1307569/150‑1 (MQ=255)
                                                                                                                           ggCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGtt                           <  2:310020/150‑1 (MQ=255)
                                                                                                                            gCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTc                          >  2:1364385/1‑150 (MQ=255)
                                                                                                                            gCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTc                          <  1:1109368/150‑1 (MQ=255)
                                                                                                                              tGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAAAAACTTCAAGTCAATCGCGCTTGGCAATGCATGTTCaa                        <  2:894898/150‑1 (MQ=255)
                                                                                                                                 tGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTCAAGCa                     >  1:2207816/1‑150 (MQ=255)
                                                                                                                                   cGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACGCCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTCAAGCAgg                   <  2:1819644/150‑1 (MQ=255)
                                                                                                                                    ggCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTCAAGCAGGc                  >  2:2209574/1‑150 (MQ=255)
                                                                                                                                    ggCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTCAAGCAGGc                  <  2:2737356/150‑1 (MQ=255)
                                                                                                                                     gCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTgca                                                                  <  2:2159197/101‑1 (MQ=255)
                                                                                                                                     gCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTgca                                                                  >  1:2159197/1‑101 (MQ=255)
                                                                                                                                      ccATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTCAAGCa                     <  2:1696456/145‑1 (MQ=255)
                                                                                                                                      ccATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTCAAGCAGGCCa                >  1:2562218/1‑150 (MQ=255)
                                                                                                                                      ccATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTAAAGTCAATCGAGCTTGGCAATGCATGTTCAAGCa                     >  1:1696456/1‑145 (MQ=255)
                                                                                                                                       aaTGCCCGTGATGCCATCGACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCaa                                                     <  1:1478914/111‑1 (MQ=255)
                                                                                                                                       aaTGCCCGTGATGCCATCGACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCaa                                                     >  2:1478914/2‑112 (MQ=255)
                                                                                                                                       aaTGCCCGTGATGCCATCGACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGCCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTCAAGCAGGCCAc               >  2:1261641/2‑150 (MQ=255)
                                                                                                                                        aTGCCCGTGATGCCATCGACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTgc                                                                   <  1:592243/97‑1 (MQ=38)
                                                                                                                                        aTGCCCGTGATGCCATCGACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTgc                                                                   >  2:592243/1‑97 (MQ=38)
                                                                                                                                              gtcgatgGCATCACGGGCATTGAACAACAGGTTGCCCCATACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTCAAGCAGGCCACGCGTAGc        <  2:1492595‑M2/143‑1 (MQ=255)
                                                                                                                                              gtcgatgGCATCACGGGCATTGAACAACAGGTTGAGCAGCACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTCAAGCAGGCCACGCGTAGc        >  1:1779119‑M2/8‑150 (MQ=255)
                                                                                                                                              gtcgatgGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTCAAGCAGGCCACGCGTAGc        <  2:773330‑M2/143‑1 (MQ=255)
                                                                                                                                                 gatgGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGa                                                                                       <  1:1229659‑M2/64‑1 (MQ=255)
                                                                                                                                                    gGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGCGCATCATTAAACATGTTCAAGCAGGCCACGCGTCGCGTTAAc  >  2:1022192‑M2/2‑150 (MQ=255)
                                                                                                                                                    gGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTCAAGCAGGCCACGCGTAGCGTTAAc  <  2:839905‑M2/149‑1 (MQ=255)
                                                                                                                                                    gGCATCACGGACATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTCAAGCAGGCCACGCGTAGCGTTAAc  <  2:1658612‑M2/149‑1 (MQ=255)
                                                                                                                                                     |                                                                                                                                                     
GCTCGAAGATGCGCGGCAGTATCTCTGTGGGTATGCCTTCGCCGCTGTCTTCCACCTCGAAGCGCACGAAGCTGCCCTTCAGCCCTTCAACTTCGCCGAACAACTCGACCACATCGACAGTCAGGCTGATGCGGCCATGCCCGTCGATGGCATCACGGGCATTGAACAACAGGTTGAGCAACACCATCTGCAGTTGCCCGGCATCCACCTCGATCAACGGCAGGTCCGGTTGCAGCAACACCTTCAAGTCAATCGAGCTTGGCAATGCATGTTCAAGCAGGCCACGCGTAGCGTTAAC  >  NC_002947/4360184‑4360481

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: