Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,066,374 G→C 5.4% Q155H (CAG→CAC PP_2678 → hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,066,3740GC 5.4% 293.5 / 4.9 93Q155H (CAG→CACPP_2678hydrolase
Reads supporting (aligned to +/- strand):  ref base G (44/44);  new base C (2/3);  total (46/47)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.85e-01

CAGGCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAACGCGCCCTGACCACGC  >  NC_002947/3066225‑3066519
                                                                                                                                                     |                                                                                                                                                 
cAGGCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAg                                                                                                                                                   >  1:2353401/1‑150 (MQ=255)
 aGGCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGc                                                                                                                                                  >  1:756745/1‑150 (MQ=255)
 aGGCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGc                                                                                                                                                  >  2:594370/1‑150 (MQ=255)
 aGGCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGc                                                                                                                                                  >  2:2073264/1‑150 (MQ=255)
  ggCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGcc                                                                                                                                                 >  1:607784/1‑150 (MQ=255)
  ggCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGcc                                                                                                                                                 <  2:1319052/150‑1 (MQ=255)
      tttGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGc                                                                                                                                             <  1:925382/150‑1 (MQ=255)
      tttGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGc                                                                                                                                             <  2:787066/150‑1 (MQ=255)
          ccGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGc                                                                                                                                         <  1:232114/150‑1 (MQ=255)
            gATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCACGCCCGGCTGCAGCACCCgggtcggcag                                                                                                                                       <  2:454218/150‑8 (MQ=255)
               gtgCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAg                                                                                                                                    <  1:780560/150‑1 (MQ=255)
                tgCCGATCGGTGCCCTGGCCGGCACCCGCGACATGCTGCGCGAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAgg                                                                                                                                   <  2:529019/150‑1 (MQ=255)
                    gATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtggg                                                                                                                               >  1:342908/1‑150 (MQ=255)
                      tCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtgggtg                                                                                                                             <  2:607784/150‑1 (MQ=255)
                       cGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGGTGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGgtgggtgg                                                                                                                            <  1:1983073/150‑1 (MQ=255)
                       cGGTGCCCTGGCCGGCACCCGCCACCGGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGgtgggtgg                                                                                                                            <  1:733086/150‑1 (MQ=255)
                         gTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGCTGCGCGGTACCGAGGTGGTGCTGCCGACCCGGCTGCAGCACCCgg                                                                                                                                               <  2:2203150/129‑1 (MQ=255)
                         gTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGCTGCAGCACCCgg                                                                                                                                               >  1:2203150/1‑129 (MQ=255)
                             ccTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCg                                                                                                                      >  2:2309873/1‑150 (MQ=255)
                               tGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtgggt                                                                                                                              <  1:1678501/140‑1 (MQ=255)
                               tGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtgggt                                                                                                                              >  2:1678501/1‑140 (MQ=255)
                               tGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAg                                                                                                                                    >  2:139797/1‑134 (MQ=255)
                               tGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAg                                                                                                                                    <  1:139797/134‑1 (MQ=255)
                                 gCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGCTGCAGCACCCgggtc                                                                                                                                            >  2:593231/1‑122 (MQ=255)
                                 gCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGCTGCAGCACCCgggtc                                                                                                                                            <  1:593231/124‑3 (MQ=255)
                                   aagCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAg                                                                                                                >  2:1902570/3‑150 (MQ=255)
                                      cACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATTCTCGGCACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGGGGTGGGTGGGCACCTGCTGAGGtt                                                                                                             >  1:1489512/1‑150 (MQ=255)
                                       aCCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTg                                                                                                            >  1:2389573/1‑150 (MQ=255)
                                         ccGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGACGGGGACTGGATGCTCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTgct                                                                                                          >  1:1499851/1‑150 (MQ=255)
                                          cgcgACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtg                                                                                                                                 >  1:2075293/1‑126 (MQ=255)
                                          cgcgACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtg                                                                                                                                 <  2:2075293/126‑1 (MQ=255)
                                           gcgACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTgctgg                                                                                                        <  2:1155675/150‑1 (MQ=255)
                                           gcgACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTgctgg                                                                                                        <  2:952946/150‑1 (MQ=255)
                                             gACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTgctgggg                                                                                                      >  2:1115427/1‑150 (MQ=255)
                                                ctgctgCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTCGGTGGGCACCGCCTGAGGTTgctggggctg                                                                                                   >  2:663418/1‑150 (MQ=255)
                                                 tgctgcGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGGCTGGTTGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGGGGGCACCGCCTGAGGTTgctggggctgg                                                                                                  >  1:1623594/1‑150 (MQ=255)
                                                 tgctgcGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTgctggggctgg                                                                                                  >  1:755742/1‑150 (MQ=255)
                                                   ctgcGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTgggc                                                                                                <  1:2309873/150‑1 (MQ=255)
                                                    tgcGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTgggcg                                                                                               <  2:1884751/150‑1 (MQ=255)
                                                     gcgcCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTgggcgg                                                                                              <  1:716962/150‑1 (MQ=255)
                                                        ccAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCa                                                                                           >  2:789985/1‑150 (MQ=255)
                                                           gCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCAt                                                                                                                                       >  1:2318363/1‑103 (MQ=255)
                                                           gCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCAt                                                                                                                                       <  2:2318363/103‑1 (MQ=255)
                                                           gCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATAc                                                                                        >  2:2185495/1‑150 (MQ=255)
                                                               ggCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGt                                                                                    <  2:2389573/150‑1 (MQ=255)
                                                                  gATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtgggtg                                                                                                                             <  1:1760891/106‑1 (MQ=255)
                                                                  gATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtgggtg                                                                                                                             >  2:1760891/1‑106 (MQ=255)
                                                                  gATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGcc                                                                                 >  2:2228099/1‑150 (MQ=255)
                                                                  gATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGAAGCTGCCGACCCGGGAGCTGCAGCCGGGCGTGAAGGAGGTGGGTGGGGACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATCGGGGTAcc                                                                                 >  2:1166293/1‑150 (MQ=255)
                                                                    tGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGa                                                                               <  2:1443213/150‑1 (MQ=255)
                                                                        aTGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTg                                                                           >  1:315153/1‑150 (MQ=255)
                                                                         tGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTgg                                                                          >  1:339742/1‑150 (MQ=255)
                                                                          ggCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGc                                                                         <  2:315153/150‑1 (MQ=255)
                                                                             cGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCa                                                                                                                                        <  1:1666049/84‑1 (MQ=255)
                                                                             cGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCa                                                                                                                                        >  2:1666049/1‑84 (MQ=255)
                                                                             cGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGAt                                                                      >  2:1412448/1‑150 (MQ=255)
                                                                                 aaCATGTATCGCCTGGTCGGCGACTGGATGTGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTc                                                                  <  2:2198294/150‑1 (MQ=255)
                                                                                    aTGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGAt                                                               >  2:847321/1‑150 (MQ=255)
                                                                                    aTGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGAt                                                               <  1:1412448/150‑1 (MQ=255)
                                                                                    aTGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGAt                                                               >  1:230515/1‑150 (MQ=255)
                                                                                      gTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGa                                                             >  1:1751072/1‑150 (MQ=255)
                                                                                        aTCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGaaa                                                           <  2:1903475/150‑1 (MQ=255)
                                                                                         tCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtgggt                                                                                                                              <  1:1436933/82‑1 (MQ=255)
                                                                                         tCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAggtgggt                                                                                                                              >  2:1436933/1‑82 (MQ=255)
                                                                                         tCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTg                                                                                                                  >  1:1014651/1‑94 (MQ=255)
                                                                                         tCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTg                                                                                                                  <  2:1014651/94‑1 (MQ=255)
                                                                                           gCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGAGCTGGCGATAATCGAGTAAAAGa                                                        >  2:1268737/1‑150 (MQ=255)
                                                                                              tGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTTCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTCAGACGTGTCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCg                                                     <  1:1268737/150‑1 (MQ=255)
                                                                                                     cGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGAGCGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCt                                              <  2:342908/150‑1 (MQ=255)
                                                                                                      gACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGa                                                                                                                 >  2:1726434/1‑82 (MQ=255)
                                                                                                      gACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGa                                                                                                                 <  1:1726434/82‑1 (MQ=255)
                                                                                                       aCTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGa                                                                >  2:650200/1‑130 (MQ=255)
                                                                                                       aCTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGa                                                                <  1:650200/130‑1 (MQ=255)
                                                                                                         tGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTc                                          <  2:1751072/150‑1 (MQ=255)
                                                                                                            aTGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGcc                                       <  1:594370/150‑1 (MQ=255)
                                                                                                             tGCGCGGTACCGAGGTGGTGCTGGCGACCCGGGTGCTGCAGGCGGGCGTGCATGAGGTGGGTGGGCACCGCCGGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCg                                      >  2:1853667/1‑150 (MQ=255)
                                                                                                              gcgcGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCAt                                                                                                                                       <  1:1732908/52‑1 (MQ=255)
                                                                                                              gcgcGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCAt                                                                                                                                       >  2:1732908/1‑52 (MQ=255)
                                                                                                                     aCCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTg                              >  2:586242/1‑150 (MQ=255)
                                                                                                                             ggtgCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGACGGGGCTGGGCGGGCTGACGGGGTCCGTCCTGGCGATTTTCTATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTa                      <  2:2078868/150‑1 (MQ=255)
                                                                                                                             ggtcCTGCCGACCCGGGTGCTGCAGGCGGGCGTGCATGTGTGGGGTGGGCACCGCCTGAGGTTGCTGGGGCTAGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTa                      <  1:1737510/150‑1 (MQ=255)
                                                                                                                              gtgCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCTCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTAc                     <  1:1341134/150‑1 (MQ=255)
                                                                                                                              gtgCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTAc                     <  1:1115427/150‑1 (MQ=255)
                                                                                                                              gtgCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTAc                     <  1:586242/150‑1 (MQ=255)
                                                                                                                                 ctgcCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGCCGGTTTTCGATGAAAAGACCGGGGTGCTTTGCGCCGGCGACCTCGTGCTCTACCaa                  >  2:542769/1‑150 (MQ=255)
                                                                                                                                  tgcCGACCCGGGTGCTGCAGGCGGGCGTGCATGACGTGGGTGCGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGAGCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAAc                 <  2:1082686/150‑1 (MQ=255)
                                                                                                                                  tgcCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAAc                 >  2:2030359/1‑150 (MQ=255)
                                                                                                                                         ccGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGCCTGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAACGCGCCCt          >  1:1039658/1‑150 (MQ=255)
                                                                                                                                         ccGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAACGCGCCCt          >  2:2274310/1‑150 (MQ=255)
                                                                                                                                         ccGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAACGCGCCCt          <  1:2030359/150‑1 (MQ=255)
                                                                                                                                               gctgcAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAACGCGCCCTGACCAc    <  1:847321/150‑1 (MQ=255)
                                                                                                                                                ctgcAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGc                                               >  1:1369348/1‑106 (MQ=255)
                                                                                                                                                ctgcAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGc                                               <  2:1369348/106‑1 (MQ=255)
                                                                                                                                                 tgcAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAACGCGCCCTGACCACGc  >  2:335000/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                 
CAGGCCTTTGCCGATGTGCCGATCGGTGCCCTGGCCGGCACCCGCGACCTGCTGCGCCAGCAAGGCGATGCCATGGCCGAGAACATGTATCGCCTGGTCGGCGACTGGATGCGCGGTACCGAGGTGGTGCTGCCGACCCGGGTGCTGCAGCCGGGCGTGCATGAGGTGGGTGGGCACCGCCTGAGGTTGCTGGGGCTGGGCGGGCATACGGGTGCCGACCTGGCGATTTTCGATGAAAAGACCGGGGTGCTGTTCGCCGGCGACCTGGTGTTCTACCAACGCGCCCTGACCACGC  >  NC_002947/3066225‑3066519

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: