Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,551,118 A→G 7.0% Q76R (CAG→CGG) ‡ PP_3135 → glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,551,1180AG 7.0% 247.6 / 5.3 85Q76R (CAG→CGG) ‡PP_3135glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base A (39/40);  new base G (3/3);  total (42/43)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.79e-01

CGCCGGGCCGACAGCGAGCGGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGC  >  NC_002947/3550970‑3551266
                                                                                                                                                    |                                                                                                                                                    
cgcCGGGCCGACAGCGAGCGGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAg                                                                                                                                                     <  1:88321/150‑1 (MQ=255)
    gggCCGACAGCGAGCGGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAgttgt                                                                                                                                                 <  2:1913295/150‑1 (MQ=255)
     ggCCGACAGCGAGCGGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAgttgtt                                                                                                                                                >  2:926012/1‑150 (MQ=255)
      gCCGACAGCGAGCGGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCTGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGttgttg                                                                                                                                               <  2:1186819/150‑1 (MQ=255)
      gCCGACAGCGAGCGGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACTGGttgttg                                                                                                                                               <  2:1572415/150‑1 (MQ=255)
       ccGACAGCGAGCGGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACTGGTTGTTgg                                                                                                                                              <  2:845815/150‑1 (MQ=255)
        cGACAGCGAGCGGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGgc                                                                                                                                             >  1:99191/1‑150 (MQ=255)
          aCAGCGAGCGGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGgcgc                                                                                                                                           >  1:1965301/1‑150 (MQ=255)
             gcgagcgGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCAcc                                                                                                                                        <  2:1965301/150‑1 (MQ=255)
               gagcgGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTc                                                                                                                                      >  1:2306767/1‑150 (MQ=255)
                          cAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGttgttg                                                                                                                                               <  1:822614/130‑1 (MQ=255)
                          cAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGttgttg                                                                                                                                               >  2:822614/1‑130 (MQ=255)
                          cAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTa                                                                                                                           <  1:1598828/150‑1 (MQ=255)
                          cAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTa                                                                                                                           >  1:516113/1‑150 (MQ=255)
                               aCCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCt                                                                                                                      <  1:1020821/150‑1 (MQ=255)
                                       gtCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTc                                                                                                              >  2:1682048/1‑150 (MQ=255)
                                          gACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAAt                                                                                                           >  2:241582/1‑150 (MQ=255)
                                               tGCCTCGCGCGATTACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGa                                                                                                      <  1:1566192/150‑1 (MQ=255)
                                                gCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGAt                                                                                                     <  2:1386893/150‑1 (MQ=255)
                                                 ccTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATg                                                                                                    <  1:1727193/150‑1 (MQ=255)
                                                      cgcgATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGcc                                                                                                  <  2:372772/147‑1 (MQ=255)
                                                      cgcgATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGcc                                                                                                  >  1:372772/1‑147 (MQ=255)
                                                      cgcgATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTc                                                                                               <  1:1394997/150‑1 (MQ=255)
                                                      cgcgATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTc                                                                                               <  1:2436190/150‑1 (MQ=255)
                                                         gATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGtt                                                                                            <  1:2129999/150‑1 (MQ=255)
                                                            gACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTgcc                                                                                         <  1:1792988/150‑1 (MQ=255)
                                                               tCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTgccgcc                                                                                      >  1:331263/1‑150 (MQ=255)
                                                               tCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTgccgcc                                                                                      <  1:241582/150‑1 (MQ=255)
                                                               tCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTAGTTgccgcc                                                                                      >  2:499283/1‑150 (MQ=255)
                                                                      aGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGc                                                                               >  2:1939485/1‑150 (MQ=255)
                                                                       gcgcATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCt                                                                              <  2:409093/150‑1 (MQ=255)
                                                                         gcATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGc                                                                            >  1:1577716/1‑150 (MQ=255)
                                                                              gCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGc                                                                               >  2:834135/1‑142 (MQ=255)
                                                                              gCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGc                                                                               <  1:834135/142‑1 (MQ=255)
                                                                                ccAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAAc                                                                     >  2:38113/1‑150 (MQ=255)
                                                                                 cAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCt                                                                                                                                       >  2:774959/1‑83 (MQ=255)
                                                                                 cAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCt                                                                                                                                       <  1:774959/83‑1 (MQ=255)
                                                                                   gcagtgCCTCGGGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCt                                                                  >  1:2208832/4‑150 (MQ=255)
                                                                                    gtgtGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACNCTc                                                                 <  2:1509117/150‑1 (MQ=255)
                                                                                     tgtgCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCg                                                                >  1:680802/1‑150 (MQ=255)
                                                                                      gtgCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGTGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTAAGCCCCAATACCGATTTCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTcta                                                               <  1:1879637/150‑3 (MQ=255)
                                                                                      gtgCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTAATGGTTTTCATCAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCAGCCGACTCGCTGCAAAAAACGCTCGa                                                               <  2:62242/150‑1 (MQ=255)
                                                                                         ccccGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGa                                                                                                      <  1:2415703/108‑1 (MQ=255)
                                                                                         ccccGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGa                                                                                                      >  2:2415703/1‑108 (MQ=255)
                                                                                           ccGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACa                                                          >  2:188096/1‑150 (MQ=255)
                                                                                           ccGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACa                                                          >  2:382194/1‑150 (MQ=255)
                                                                                            cGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACAt                                                         <  1:1682048/150‑1 (MQ=255)
                                                                                             gAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATg                                                        <  1:38113/150‑1 (MQ=255)
                                                                                                    tgctgctgAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGccc                                                 >  1:1854113/1‑150 (MQ=255)
                                                                                                       tgctgAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCAcc                                              >  1:889798/1‑150 (MQ=255)
                                                                                                       tgctgAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCAcc                                              >  2:2122840/1‑150 (MQ=255)
                                                                                                        gctgAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCg                                                                                <  1:2035761/115‑1 (MQ=255)
                                                                                                        gctgAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCg                                                                                >  2:2035761/1‑115 (MQ=255)
                                                                                                         ctgAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCg                                            >  2:1826475/1‑150 (MQ=255)
                                                                                                         ctgAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCg                                            <  1:926012/150‑1 (MQ=255)
                                                                                                         ctgAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCAGAAAAACTGAAAAAAACGCTCGACTACATGGACGCCCACCCg                                            <  2:830182/150‑1 (MQ=255)
                                                                                                              caacaaGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGTCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCAcc                                              <  2:1924497/143‑1 (MQ=255)
                                                                                                              caacaaGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCAcc                                              >  1:1924497/1‑143 (MQ=255)
                                                                                                              caacaaGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCAcc                                              <  2:1924735/143‑1 (MQ=255)
                                                                                                              caacaaGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCAcc                                              >  1:1924735/1‑143 (MQ=255)
                                                                                                              caacaaGCAGAACGTGGGTTTCGGCCGGGCCAAAAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTg                                       >  1:564388/1‑150 (MQ=255)
                                                                                                                acaaGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTgcg                                     >  1:240606/1‑150 (MQ=255)
                                                                                                                  aatccgAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACCTGGACGCCCACCCGGAGTgcggc                                   <  2:1181313/145‑1 (MQ=255)
                                                                                                                           gTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTgccgcc                                                                                      <  2:2292167/90‑1 (MQ=255)
                                                                                                                           gTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTgccgcc                                                                                      >  1:2292167/1‑90 (MQ=255)
                                                                                                                           gTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTgggg                          <  1:188096/150‑1 (MQ=255)
                                                                                                                           gTGGGTTGCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTgggg                          >  2:901775/1‑150 (MQ=255)
                                                                                                                            tGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTggggg                         <  2:1870545/150‑1 (MQ=255)
                                                                                                                             gggTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTAAAGGGCAAGTATGCCTTGCTGCTCAATACCTATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGt                        >  1:1226890/1‑150 (MQ=255)
                                                                                                                               gTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTgc                      <  2:516113/150‑1 (MQ=255)
                                                                                                                                     gCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCGCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTgg                >  1:870976/1‑150 (MQ=255)
                                                                                                                                      ccGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTc                                                                 <  2:1515459/100‑1 (MQ=255)
                                                                                                                                      ccGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTc                                                                 >  1:1515459/1‑100 (MQ=255)
                                                                                                                                      ccGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGa               <  1:1715490/150‑1 (MQ=255)
                                                                                                                                       cGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAg              <  2:2208832/150‑1 (MQ=255)
                                                                                                                                           ccAACAACTGGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATAcc                                                                                                        >  1:705770/1‑56 (MQ=38)
                                                                                                                                           ccAACAACTGGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATAcc                                                                                                        <  2:705770/56‑1 (MQ=38)
                                                                                                                                           ccAACAACTGGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTc          >  1:2339422/1‑150 (MQ=255)
                                                                                                                                           ccAACAACTGGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTc          >  1:162487/1‑150 (MQ=255)
                                                                                                                                            caacaaCCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTcg         >  1:1299524/1‑150 (MQ=255)
                                                                                                                                             aacaacCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTa                                >  2:1877535/1‑126 (MQ=255)
                                                                                                                                             aacaacCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTa                                <  1:1877535/126‑1 (MQ=255)
                                                                                                                                                aacCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGAt     >  1:2035400/1‑150 (MQ=255)
                                                                                                                                                  ccAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATgg   >  1:791253/1‑150 (MQ=255)
                                                                                                                                                   cAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGACGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGc  <  2:680802/150‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
CGCCGGGCCGACAGCGAGCGGCTGGCCAGTTACCTGGTGGTCGACAATGCCTCGCGCGATGACTCGGTGCAGCGCATGGCCCAGGTGTGCCCCGAGGCACTGCTGCTGAGCAACAAGCAGAACGTGGGTTTCGGCCGGGCCAACAACCAGTTGTTGGCGCACCTCAAGGGCAAGTATGCCTTGCTGCTCAATACCGATGCCTTCGTTGCCGCCGACTCGCTGCAAAAAACGCTCGACTACATGGACGCCCACCCGGAGTGCGGCGTACTGGGGGTGCGCCTGGAGGGTCGCGATGGC  >  NC_002947/3550970‑3551266

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: