Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,219,156 C→T 6.7% D229D (GAC→GAT) ‡ 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,1560CT 6.7% 213.8 / 5.0 75D229D (GAC→GAT) ‡PP_1064alpha/beta family hydrolase
Reads supporting (aligned to +/- strand):  ref base C (33/37);  new base T (2/3);  total (35/40)
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.67e-01

GACGCGTTGGGCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAGGCGGCGCAGGCC  >  NC_002947/1219007‑1219300
                                                                                                                                                     |                                                                                                                                                
gACGCGTTGGGCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGAc                                                                                                                                                  >  2:516299/1‑150 (MQ=255)
gACGAGTTGGGCGTGCGCCAGTTTCATGTATTGGGTGTGAGCCTCGGCGGAATCACTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGAc                                                                                                                                                  <  2:38826/150‑1 (MQ=255)
    cgTTGGGCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCgg                                                                                                                                              <  2:989155/150‑1 (MQ=255)
       tGGGCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGgcgc                                                                                                                                           <  1:171481/150‑1 (MQ=255)
          gCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGGTCCGGCGCCgg                                                                                                                                        <  2:100940/150‑1 (MQ=255)
          gCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTgg                                                                                                                                        <  2:143402/150‑1 (MQ=255)
               cgcCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAg                                                                                                                                   <  2:1215877/150‑1 (MQ=255)
                gcCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTg                                                                                                                                      <  2:1386612/146‑1 (MQ=255)
                gcCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGt                                                                                                                                  >  1:890775/1‑150 (MQ=255)
                gcCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGGTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTg                                                                                                                                      >  1:1386612/1‑146 (MQ=255)
                 ccAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGtt                                                                                                                                 >  2:698794/1‑150 (MQ=255)
                  cAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCAGCGCCGGGTCCGgcgccg                                                                                                                                         <  1:857712/141‑3 (MQ=255)
                  cAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCAGCGCCGGGTCCGgcgccg                                                                                                                                         >  2:857712/1‑139 (MQ=255)
                   aGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGc                                                                                                                               <  1:516299/150‑1 (MQ=255)
                    gTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCt                                                                                                                              <  1:698794/150‑1 (MQ=255)
                    gTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCt                                                                                                                              <  1:534942/150‑1 (MQ=255)
                     tttCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTg                                                                                                                             <  2:1424240/150‑1 (MQ=255)
                     tttCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTg                                                                                                                             >  2:715433/1‑150 (MQ=255)
                     tttCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCAAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTg                                                                                                                             <  1:747420/150‑1 (MQ=255)
                       tCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGc                                                                                                                           >  1:1202839/1‑150 (MQ=255)
                       tCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGc                                                                                                                           >  1:991984/1‑150 (MQ=255)
                         aTGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGccc                                                                                                                         >  2:1193510/1‑150 (MQ=255)
                         aTGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGccc                                                                                                                         >  2:898875/1‑150 (MQ=255)
                         aTGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGccc                                                                                                                         >  2:223304/1‑150 (MQ=255)
                          tGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGACTCCCGGCGCCGGACCCGGCGCTg                                                                                                                                         <  1:648676/133‑1 (MQ=255)
                          tGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGACTCCCGGCGCCGGACCCGGCGCTg                                                                                                                                         >  2:648676/1‑133 (MQ=255)
                            tATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTgcgccggg                                                                                                                      <  2:890775/150‑8 (MQ=255)
                             aTTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGTCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCGCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGgcg                                                                                                                     <  2:624157/150‑1 (MQ=255)
                             aTTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTg                                                                                                                             <  2:234993/142‑1 (MQ=255)
                             aTTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTg                                                                                                                             >  1:234993/1‑142 (MQ=255)
                              ttGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTg                                                                                                                             >  1:418158/1‑141 (MQ=255)
                              ttGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTg                                                                                                                             <  2:418158/141‑1 (MQ=255)
                                 ggTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGgcgcagc                                                                                                                 <  1:261248/150‑1 (MQ=255)
                                  gtgtGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGgcgcagcg                                                                                                                <  1:1193510/150‑1 (MQ=255)
                                  gtgtGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGCGCGCAGACCGTCCCGGCGCCGGACCCAGCGGTTGGGCAGTGGCGGGCCCGgtgcagcg                                                                                                                >  2:514859/1‑150 (MQ=255)
                                   tgtgAGCATGGGCGGCATGATTGCCAAGAACCTGACGGGGAGCGCGGCCGAGCGTGTGAGACCCCTGAAGCTGGCACTGTCCGGCCCCTTGGTCGCGGGTCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTCGCTGGCCCGaggcagcgc                                                                                                               >  1:324942/1‑150 (MQ=255)
                                              gcggcATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTg                                                                                                    >  2:613166/1‑150 (MQ=255)
                                                       ttGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCaa                                                                                           >  1:261651/1‑150 (MQ=255)
                                                            cAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGag                                                                                      >  2:1053593/1‑150 (MQ=255)
                                                                   tGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAgg                                                                               >  1:1112764/1‑150 (MQ=255)
                                                                   tGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAgg                                                                               >  2:39581/1‑150 (MQ=255)
                                                                     gcggcgATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGcc                                                                             <  1:898875/150‑1 (MQ=255)
                                                                                  ccGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTg                                                                                                    >  2:1165682/1‑114 (MQ=255)
                                                                                  ccGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTg                                                                                                    <  1:1165682/114‑1 (MQ=255)
                                                                                    gAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGcc                                                              >  2:144146/1‑150 (MQ=255)
                                                                                     aGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGTCGGCctgctgctcgatt                                                             <  1:39581/150‑8 (MQ=255)
                                                                                             cgcAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGc                                                     <  1:362726/150‑1 (MQ=255)
                                                                                                tgCCTGACGCTGGTACTGTCCTGTTCCGGGGTGGCGGTGCGTCCGGCGTGTGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCt                                                  <  2:1084384/149‑1 (MQ=255)
                                                                                                 gCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTg                                                 <  2:1028525/150‑1 (MQ=255)
                                                                                                  ccTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGTGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGa                                                >  1:1084384/1‑150 (MQ=255)
                                                                                                  ccTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGa                                                <  2:991984/150‑1 (MQ=255)
                                                                                                      aCGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTa                                            <  2:709576/150‑1 (MQ=255)
                                                                                                       cGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCTAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTAc                                           <  1:352582/150‑1 (MQ=255)
                                                                                                          tGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACggg                                        >  2:1152531/1‑150 (MQ=255)
                                                                                                            gTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGgac                                      <  1:1053593/150‑1 (MQ=255)
                                                                                                             tCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGAGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGgacg                                     >  1:1045832/1‑150 (MQ=255)
                                                                                                               aTGTCCAGCTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGga                                       >  2:607456/1‑146 (MQ=255)
                                                                                                               aTGTCCAGCTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGga                                       <  1:607456/146‑1 (MQ=255)
                                                                                                                tGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTggcggc                                                               <  2:1382017/121‑1 (MQ=255)
                                                                                                                tGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTggcggc                                                               >  1:1382017/1‑121 (MQ=255)
                                                                                                                 gTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACcg                                 >  1:859647/1‑150 (MQ=255)
                                                                                                                  tCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGagca                                                                                    <  1:746905/98‑1 (MQ=255)
                                                                                                                  tCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGagca                                                                                    >  2:746905/1‑98 (MQ=255)
                                                                                                                    cAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACTTTGAAGTGGCAATCGAGCAGCAGGCCGAACTGAAGGCAGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGa                              >  1:318270/1‑150 (MQ=255)
                                                                                                                     aGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAg                             <  2:324942/150‑1 (MQ=255)
                                                                                                                     aGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGAACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAg                             <  1:715433/150‑1 (MQ=255)
                                                                                                                       ttCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAg                             <  2:891774/148‑1 (MQ=255)
                                                                                                                       ttCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAg                             >  1:891774/1‑148 (MQ=255)
                                                                                                                           ggggCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTtg                       <  2:1112764/150‑1 (MQ=255)
                                                                                                                            gggCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTtgc                      <  1:613166/150‑1 (MQ=255)
                                                                                                                             ggcggcGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTtgct                     <  2:261651/150‑1 (MQ=255)
                                                                                                                                       cTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAggcg           >  2:1160305/1‑150 (MQ=255)
                                                                                                                                         accagcgcCGGGTCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCCGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTggg                                                         >  2:1034419/5‑102 (MQ=255)
                                                                                                                                         accagcgcCGGGTCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCCCCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTggg                                                         <  1:1034419/98‑1 (MQ=255)
                                                                                                                                                ccGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAGGCCGCGCAGGcc  >  2:812730/1‑150 (MQ=255)
                                                                                                                                                   gACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATc                                                                                         <  1:983219/60‑1 (MQ=255)
                                                                                                                                                   gACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATc                                                                                         >  2:983219/1‑60 (MQ=255)
                                                                                                                                                     |                                                                                                                                                
GACGCGTTGGGCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAGGCGGCGCAGGCC  >  NC_002947/1219007‑1219300

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