Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,219,158 T→C 6.8% L232P (CTG→CCG)  PP_1064 → alpha/beta family hydrolase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029471,219,1580TC 6.8% 255.5 / 6.8 88L232P (CTG→CCG) PP_1064alpha/beta family hydrolase
Reads supporting (aligned to +/- strand):  ref base T (43/39);  new base C (3/3);  total (46/42)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.69e-01

GCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAGGCGGCGCAGGCCTACGACCGGGCG  >  NC_002947/1219011‑1219306
                                                                                                                                                    |                                                                                                                                                     
gCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTg                                                                                                                                                      <  5:235746/149‑1 (MQ=255)
gCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTg                                                                                                                                                      <  4:322706/149‑1 (MQ=255)
  gTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGt                                                                                                                                                    >  8:148228/1‑149 (MQ=255)
  gTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGt                                                                                                                                                    <  3:13373/149‑1 (MQ=255)
     cgcCAGTTTCATGTTTGGGGTGGGAGCATGGGCGGCATGATTGCCCAGCACCGGGCGGCGCTGGCGCCCGAGCGTGGGCGCAGCCTGCCGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGGTCCGGCGCCgggagc                                                                                                                                                 <  3:29232/149‑5 (MQ=255)
      gcCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAg                                                                                                                                               <  7:317338/149‑1 (MQ=255)
             tCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTgg                                                                                                                                        >  1:183678/1‑149 (MQ=255)
                 gTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCg                                                                                                                                    >  8:66245/1‑149 (MQ=255)
                  tATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCgg                                                                                                                                   >  6:177995/1‑149 (MQ=255)
                  tATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCgg                                                                                                                                   <  2:213472/149‑1 (MQ=255)
                  tATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCg                                                                                                                                    <  1:51532/148‑1 (MQ=255)
                    ttGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGGTCCGGCGCCggga                                                                                                                                                   >  6:148736/1‑130 (MQ=255)
                    ttGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGGTCCGGCGCCggga                                                                                                                                                   <  5:148736/132‑3 (MQ=255)
                    ttGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGgcg                                                                                                                                 >  4:38852/1‑149 (MQ=255)
                    ttGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCCAGCGTGTGAGCAGCCTGACGCTGGTAATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGgcg                                                                                                                                 >  8:269181/1‑149 (MQ=255)
                       ggTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGgcgcag                                                                                                                              >  7:266105/1‑149 (MQ=255)
                        gggtgATCATGGGCGGCATGATTGCCCAGCTCCTGGCTGCGATGGCGCCCGAGCGTGTGCGCAGCATGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGgcgcagc                                                                                                                             <  4:15198/147‑1 (MQ=255)
                           tgAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTgg                                                                                                                                                     <  7:264307/123‑1 (MQ=255)
                           tgAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTgg                                                                                                                                                     >  8:264307/1‑123 (MQ=255)
                           tgAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGgcgcagcgca                                                                                                                          >  7:246039/1‑149 (MQ=255)
                             aGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCAcc                                                                                                                        >  6:198666/1‑149 (MQ=255)
                              gCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGc                                                                                                                                           >  3:23357/1‑129 (MQ=255)
                              gCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGc                                                                                                                                           <  4:23357/129‑1 (MQ=255)
                              gCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCAcc                                                                                                                        >  5:195988/1‑148 (MQ=255)
                              gCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCAcc                                                                                                                        >  5:50647/1‑148 (MQ=255)
                                    gcggcATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGt                                                                                                                 >  6:131211/1‑149 (MQ=255)
                                        cATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAg                                                                                                             >  8:177656/1‑149 (MQ=255)
                                            aTTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTgg                                                                                                          <  2:183678/148‑1 (MQ=255)
                                                cccAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAAt                                                                                                      <  2:141096/148‑1 (MQ=255)
                                                cccAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATc                                                                                                     <  2:223258/149‑1 (MQ=255)
                                                 ccAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATc                                                                                                     >  2:209784/1‑148 (MQ=255)
                                                 ccAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCg                                                                                                    >  6:267778/1‑149 (MQ=255)
                                                  cAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGa                                                                                                   <  5:339377/149‑1 (MQ=255)
                                                         tGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGcagcag                                                                                            >  5:66985/1‑149 (MQ=255)
                                                          ggcggcGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAgg                                                                                           >  5:344872/1‑149 (MQ=255)
                                                           gcggcgATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGc                                                                                          <  8:246039/149‑1 (MQ=255)
                                                              gcgATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGGGCATGTCCGGTGCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGATGGCCCGGCGCAGCGCACCGACCCGTGAAGTGGCAATCGAGCTGCAGGCCg                                                                                        >  1:246704/1‑149 (MQ=255)
                                                                     cgcCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACctgctg                                                                                >  5:162641/1‑149 (MQ=255)
                                                                      gcCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTgg                                                                               <  5:131211/149‑1 (MQ=255)
                                                                        ccGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCCGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTggcg                                                                             >  8:178873/1‑149 (MQ=255)
                                                                        ccGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTggcg                                                                             >  5:308847/1‑149 (MQ=255)
                                                                        ccGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTggcg                                                                             >  6:91772/1‑149 (MQ=255)
                                                                             cGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCtt                                                                        <  5:267778/149‑1 (MQ=255)
                                                                              gtgtGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTg                                                                       <  6:222067/149‑1 (MQ=255)
                                                                                gtgCGCAGCCTGCCGCGGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTggg                                                                     <  3:297357/149‑1 (MQ=255)
                                                                                gtgCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTggg                                                                     <  3:38852/149‑1 (MQ=255)
                                                                                    gcAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGc                                                                 <  6:50647/149‑1 (MQ=255)
                                                                                       gCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCt                                                              <  6:195988/149‑1 (MQ=255)
                                                                                         cTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGa                                                            <  7:148228/149‑1 (MQ=255)
                                                                                          tGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAg                                                           >  4:273800/1‑149 (MQ=255)
                                                                                            aCGCTGGTCATGTCCAGTTCAGGGGTGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTggc                                                                              <  4:30567/128‑1 (MQ=255)
                                                                                            aCGCTGGTCATGTCCAGTTCAGGGGTGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTggc                                                                              >  3:30567/1‑128 (MQ=255)
                                                                                                   gTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACggg                                                    <  8:266105/148‑1 (MQ=255)
                                                                                                   gTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGga                                                   >  7:278906/1‑149 (MQ=255)
                                                                                                      aTGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGgacga                                                >  5:230159/1‑149 (MQ=255)
                                                                                                      aTGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGgacga                                                <  5:198666/149‑1 (MQ=255)
                                                                                                        gTCCAGTTCAGGGGCGGCGGGGCTCACGGCGCCGGACCCGGCGCTGGTGGAAGTTGATGGCCCGGCGCAGCGCACCGAACCGTCAAGTGGCAATCGAGCAGAAGCCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGCTACGGgacgac                                               >  1:271508/1‑149 (MQ=255)
                                                                                                         tCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACcg                                             >  1:228440/1‑149 (MQ=255)
                                                                                                          ccAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACcgc                                            >  8:110812/1‑149 (MQ=255)
                                                                                                            aGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCGGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGa                                          <  7:66245/149‑1 (MQ=255)
                                                                                                              ttCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAgg                                        >  7:217347/1‑149 (MQ=255)
                                                                                                              ttCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAgg                                        <  7:269181/149‑1 (MQ=255)
                                                                                                                 aGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGCCCTGCTGGCGGCCATGGGCAGCCCTGAGGTAAGGGACGACCGCGAGgtgt                                     >  3:30049/1‑149 (MQ=255)
                                                                                                                 aGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGgtgt                                     <  5:91772/149‑1 (MQ=255)
                                                                                                                   gggCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTtg                                   <  6:344872/149‑1 (MQ=255)
                                                                                                                   gggCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTtg                                   >  8:205137/1‑149 (MQ=255)
                                                                                                                     gcggcgGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTTAGGTACGGGACGACCGCGAGGTGTtgct                                 <  1:209784/149‑1 (MQ=255)
                                                                                                                     gcggcgGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTtgct                                 <  3:273800/149‑1 (MQ=255)
                                                                                                                       ggcggGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTtgctgc                               >  7:244819/1‑149 (MQ=255)
                                                                                                                             gCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGc                                                                 >  7:346730/1‑109 (MQ=255)
                                                                                                                             gCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGc                                                                 <  8:346730/109‑1 (MQ=255)
                                                                                                                             gCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAg                                         <  7:196948/133‑1 (MQ=255)
                                                                                                                             gCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAg                                         >  8:196948/1‑133 (MQ=255)
                                                                                                                                accagcgcCGGGTCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAggcgg                      >  7:39287/5‑149 (MQ=255)
                                                                                                                                 ccagcgcCGGGTCCGGCGCCGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAggcggc                     >  5:39134/4‑149 (MQ=255)
                                                                                                                                 ccGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGCCCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAggcggc                     >  8:116937/1‑149 (MQ=255)
                                                                                                                                 ccGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAggcg                       <  7:110812/147‑1 (MQ=255)
                                                                                                                                 ccGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAggcg                       <  8:217347/147‑1 (MQ=255)
                                                                                                                                  cagcgcCGGGTCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTggg                                                                     >  8:231021/3‑100 (MQ=255)
                                                                                                                                  cGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAGgcggcg                    >  7:166467/1‑149 (MQ=255)
                                                                                                                                    gcgcCGGGTCCGGCGCCGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAGgcggcg                    >  6:60698/1‑147 (MQ=255)
                                                                                                                                     cgcCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACGGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGTTACGGGACGACCGCGAGGTGTTGCTGCACCAGGCGGCGCAg                 <  6:308847/149‑1 (MQ=255)
                                                                                                                                     cgcCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAGGCGGCGCAg                 <  6:230159/149‑1 (MQ=255)
                                                                                                                                      gcCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGCACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCGTGGGCAGCCCTGAGGTACGGGACGCCCGCGAGGTGTTGCTGCACCAGGCGGCGCAgg                >  5:248842/1‑149 (MQ=255)
                                                                                                                                             ccGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAGGCGGCGCAGGCCTACGa         <  8:244819/149‑1 (MQ=255)
                                                                                                                                             ccGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAGGCGGCGCAGGCCTACGa         >  3:122583/1‑149 (MQ=255)
                                                                                                                                                 cgcTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAGGCGGCGCAGGCCTACGACCgg     <  7:116937/149‑1 (MQ=255)
                                                                                                                                                 cgcCGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAGGCGGCGCAGGCCTACGACCgg     <  5:60698/149‑1 (MQ=255)
                                                                                                                                                    tGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAGGCGGCGCAGGCCTACGACCGGGCg  >  8:58019/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                     
GCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAGGCGGCGCAGGCCTACGACCGGGCG  >  NC_002947/1219011‑1219306

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: