Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,551,261 G→C 8.5% D124H (GAT→CAT)  PP_3135 → glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,551,2610GC 8.5% 253.6 / 9.0 82D124H (GAT→CAT) PP_3135glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (37/38);  new base C (2/5);  total (39/43)
Fisher's exact test for biased strand distribution p-value = 4.36e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.91e-01

AACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGTGAAGATGGTCGACGACATGAACTGGGACCACGCCTCGGTGCGCGAGTGTGA  >  NC_002947/3551114‑3551403
                                                                                                                                                   |                                                                                                                                              
aacCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGAt                                                                                                                                              <  1:143717/150‑1 (MQ=255)
aacCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCCAt                                                                                                                                              <  2:2363101/150‑1 (MQ=255)
     gttgttGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGa                                                                                                                                         >  2:367397/1‑150 (MQ=255)
       tgttgGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCCATCGCGAcc                                                                                                                                       <  1:2179876/150‑1 (MQ=255)
         ttgGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGc                                                                                                                                           <  1:486085/144‑1 (MQ=255)
         ttgGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGc                                                                                                                                           >  2:486085/1‑144 (MQ=255)
           ggCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCCATCGCGAccctcc                                                                                                                                   <  2:2277025/150‑5 (MQ=255)
            gcgcACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAg                                                                                                                                  >  1:1967408/1‑150 (MQ=255)
            gcgcACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCCATCGCGAccctcca                                                                                                                                  <  1:724841/150‑6 (MQ=255)
              gcACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCa                                                                                                                                   <  2:2310958/147‑1 (MQ=255)
              gcACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCa                                                                                                                                   >  1:2310958/1‑147 (MQ=255)
                   tCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGtt                                                                                                                           >  2:1150664/1‑150 (MQ=255)
                   tCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGATGCCGCCGACTCGCTGCAAAAAACGATCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGAGGTCTCGCGATGGCGACCTGCAGCCCAGtt                                                                                                                           >  2:1009277/1‑150 (MQ=255)
                     aaGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGt                                                                                                                         >  2:576306/1‑150 (MQ=255)
                      aGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTc                                                                                                                        <  2:2256022/150‑1 (MQ=255)
                      aGGGCAAGTATGCCGTGCTGCTCAATACCGGTGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTc                                                                                                                        <  1:289519/150‑1 (MQ=255)
                        ggCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGc                                                                                                                      >  1:1798332/1‑150 (MQ=255)
                         gCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCt                                                                                                                     >  1:2137919/1‑150 (MQ=255)
                          cAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTa                                                                                                                    >  1:317852/1‑150 (MQ=255)
                                tGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTccc                                                                                                              >  1:1921633/1‑150 (MQ=255)
                                 gCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTcccc                                                                                                             >  2:1961589/1‑150 (MQ=255)
                                  ccTTGCGGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCca                                                                                                            >  2:348264/1‑150 (MQ=255)
                                   cTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCCATCGCGAc                                                                                                                                        <  2:1313614/121‑1 (MQ=255)
                                   cTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCCATCGCGAc                                                                                                                                        >  1:1313614/1‑121 (MQ=255)
                                    ttGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCcaca                                                                                                          >  2:1841168/1‑150 (MQ=255)
                                        tgctCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGc                                                                                                      >  2:1916936/1‑150 (MQ=255)
                                         gctCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCt                                                                                                     <  2:47001/150‑1 (MQ=255)
                                         gctCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCt                                                                                                     <  2:2304855/150‑1 (MQ=255)
                                          ctCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTc                                                                                                    <  2:1798332/150‑1 (MQ=255)
                                           tCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGc                                                                                                                                 <  2:62083/120‑1 (MQ=255)
                                           tCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACAGGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGc                                                                                                                                 >  1:62083/1‑120 (MQ=255)
                                             aaTACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAAt                                                                                                 <  1:1715364/150‑1 (MQ=255)
                                              aTACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAACACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGACGCGATGGCGACCTGCAGCCAAGTTGTCGCTACTTCCCCACACCGCTCacgg                                                                                                >  2:279526/1‑147 (MQ=255)
                                               tACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAAtgt                                                                                               <  1:1961589/150‑1 (MQ=255)
                                                   gATGCCTCCGTTGCGGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCGAGCCGGCGTGCGGCGTACTGGGGGTCCGCCTGGAGGGTCGCGATGGCTACCTGCAGCCCATTTGTCGCTACTTCCCCACACCGCTCAATGTGTtc                                                                                           <  1:1150664/150‑1 (MQ=255)
                                                    aTGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCGGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTtcg                                                                                          >  1:871893/1‑150 (MQ=255)
                                                       ccTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTtcgtcg                                                                                       <  2:712458/150‑1 (MQ=255)
                                                        cTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCgg                                                                                      <  2:536099/150‑1 (MQ=255)
                                                        cTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCgg                                                                                      <  2:1921633/150‑1 (MQ=255)
                                                         ttCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCggg                                                                                     >  1:285776/1‑150 (MQ=255)
                                                          tCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGc                                                                                    <  1:576306/150‑1 (MQ=255)
                                                                  gccgACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCggg                                                                            >  1:1671425/1‑150 (MQ=255)
                                                                           cTGCAACAAACGCGCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCgg                                                                   <  1:279526/150‑1 (MQ=255)
                                                                           cTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCgg                                                                   <  1:367397/150‑1 (MQ=255)
                                                                            tGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGt                                                                  >  2:1026482/1‑150 (MQ=255)
                                                                            tGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGt                                                                  <  2:1821326/150‑1 (MQ=255)
                                                                               aaaaaaCGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGttct                                                               >  2:614465/1‑150 (MQ=255)
                                                                               aaaaaaCGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGttct                                                               >  1:1966410/1‑150 (MQ=255)
                                                                                  aaaCGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTcc                                                            <  1:597540/150‑1 (MQ=255)
                                                                                    aCGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACATCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGAGGTTCTTCCCa                                                          <  2:1609784/150‑1 (MQ=255)
                                                                                      gCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAgg                                                        <  2:2148843/150‑1 (MQ=255)
                                                                                       cTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAggg                                                       <  2:1966410/150‑1 (MQ=255)
                                                                                        tCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAgggg                                                      >  1:1743512/1‑150 (MQ=255)
                                                                                            cTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCgg                                                                                      >  2:615506/1‑114 (MQ=255)
                                                                                            cTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCgg                                                                                      <  1:615506/114‑1 (MQ=255)
                                                                                            cTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGTGaa                                                  >  2:224420/1‑150 (MQ=255)
                                                                                                aTGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGTGAAGATg                                              <  2:1967408/150‑1 (MQ=255)
                                                                                                   gACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGTGAAGATGGTc                                           <  1:224420/150‑1 (MQ=255)
                                                                                                    aCGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGTGAAGATGGTcg                                          >  1:164742/1‑150 (MQ=255)
                                                                                                      gCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGTGAAGATGGTcgac                                        <  1:1841168/150‑1 (MQ=255)
                                                                                                        ccACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGtt                                                                                            <  1:394068/96‑1 (MQ=255)
                                                                                                        ccACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGtt                                                                                            >  2:394068/1‑96 (MQ=255)
                                                                                                        ccACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGTGAAGATGGTcgacga                                      >  1:1420171/1‑150 (MQ=255)
                                                                                                            ccGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGTGAAGATGGTCGACGACATg                                  >  2:1595202/1‑150 (MQ=255)
                                                                                                            ccGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGTGAAGATGGTCGACGACATg                                  >  2:274646/1‑150 (MQ=255)
                                                                                                            ccGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGTGAAGATGGTCGACGACATg                                  >  2:1703698/1‑150 (MQ=255)
                                                                                                            ccGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGTGAAGATGGTCGACGACATg                                  <  2:1540848/150‑1 (MQ=255)
                                                                                                                  tGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGTGAAGATGGTCGACGACATGAACTgg                            <  1:1703698/150‑1 (MQ=255)
                                                                                                                  tGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGTGAAGATGGTCGACGACATGAACTgg                            <  2:285776/150‑1 (MQ=255)
                                                                                                                              ggggTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGTGAAGATGGTCGACGACATGAACTGGGAc                         <  2:1589937/141‑1 (MQ=255)
                                                                                                                              ggggTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGTGAAGATGGTCGACGACATGAACTGGGAc                         >  1:1589937/1‑141 (MQ=255)
                                                                                                                                ggTGCGGCTGGAGGGTCGCGATGGCGCCCGGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGTGAAGATGGTCGACGACATGAACTGGGACCACGCCTCGGt              <  1:1109186/150‑1 (MQ=255)
                                                                                                                                    cgcCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGTGAAGATGGTCGACGACATGAACTGGGACCACGCCTCGGTgcgc          <  2:2137919/150‑1 (MQ=255)
                                                                                                                                    cgcCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGTGAAGATGGTCGACGACATGAACTGGGACCACGCCTCGGTgcgc          <  2:1700356/150‑1 (MQ=255)
                                                                                                                                     gcCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGTGAAGATGGTCGACGACATGAACTGGGACCGCGCCTCGGTGcgcg         >  2:897317/1‑150 (MQ=255)
                                                                                                                                       cTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGc                                                                                                                      <  2:586669/39‑1 (MQ=255)
                                                                                                                                       cTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGc                                                                                                                      >  1:586669/1‑39 (MQ=255)
                                                                                                                                        tGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGTGAAGATGGTCGACGACATGAACTGGGACCACGCCTCGGTGCGCGAgt      >  1:1657271/1‑150 (MQ=255)
                                                                                                                                         tgcaggTCGCCATCGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGAGAAGATGGTCGACGACATGAACTGGGACCACGCCTCGGTGCGCGAgtg     >  2:1760664/5‑150 (MQ=255)
                                                                                                                                            gggTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGTGAAGATGGTCGACGACATGAACTGGGa                          <  2:524780/126‑1 (MQ=255)
                                                                                                                                            gggTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGTGAAGATGGTCGACGACATGAACTGGGa                          >  1:524780/1‑126 (MQ=255)
                                                                                                                                            gggTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGTGAAGATGGTCGACGACATGAACTGGGACCACGCCTCGGTGCGCGAGTGTGa  <  1:1916936/150‑1 (MQ=255)
                                                                                                                                                   |                                                                                                                                              
AACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGTGAAGATGGTCGACGACATGAACTGGGACCACGCCTCGGTGCGCGAGTGTGA  >  NC_002947/3551114‑3551403

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