Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,551,264 G→C 7.5% G125R (GGC→CGC)  PP_3135 → glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,551,2640GC 7.5% 244.2 / 7.5 80G125R (GGC→CGC) PP_3135glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base G (37/37);  new base C (2/4);  total (39/41)
Fisher's exact test for biased strand distribution p-value = 6.76e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.49e-01

GTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGTGAAGATGGTCGACGACATGAACTGGGACCACGCCTCGGTGCGCGAGTGTGACTGGTTGCCC  >  NC_002947/3551119‑3551413
                                                                                                                                                 |                                                                                                                                                     
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                                                                                                                                             >  1:2310958/1‑147 (MQ=255)
         gcACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCa                                                                                                                                             <  2:2310958/147‑1 (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                                                                                                                                                  >  1:1313614/1‑121 (MQ=255)
                              cTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCCATCGCGAc                                                                                                                                                  <  2:1313614/121‑1 (MQ=255)
                               ttGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCcaca                                                                                                                    >  2:1841168/1‑150 (MQ=255)
                                   tgctCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGc                                                                                                                >  2:1916936/1‑150 (MQ=255)
                                    gctCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCt                                                                                                               <  2:2304855/150‑1 (MQ=255)
                                    gctCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCt                                                                                                               <  2:47001/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                                                                         >  1:1966410/1‑150 (MQ=255)
                                                                          aaaaaaCGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGttct                                                                         >  2:614465/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:1540848/150‑1 (MQ=255)
                                                                                                       ccGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGTGAAGATGGTCGACGACATg                                            >  2:1595202/1‑150 (MQ=255)
                                                                                                       ccGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGTGAAGATGGTCGACGACATg                                            >  2:274646/1‑150 (MQ=255)
                                                                                                       ccGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGTGAAGATGGTCGACGACATg                                            >  2:1703698/1‑150 (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)
                                                                                                                                                 ggCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGTGAAGATGGTCGACGACATGAACTGGGACCACGCCTCGGTGCGCGAGTGTGACTGGTTGccc  <  2:1743512/150‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                     
GTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGCGACCTGCAGCCCAGTTGTCGCTACTTCCCCACACCGCTCAATGTGTTCGTCGGGCGTACCGGGTTGGGGCGGTTCTTCCCAGGGGTGAAGATGGTCGACGACATGAACTGGGACCACGCCTCGGTGCGCGAGTGTGACTGGTTGCCC  >  NC_002947/3551119‑3551413

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