Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 1,219,156 C→T 8.6% 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 8.6% 197.2 / 6.8 70D229D (GAC→GAT) ‡PP_1064alpha/beta family hydrolase
Reads supporting (aligned to +/- strand):  ref base C (33/31);  new base T (3/3);  total (36/34)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.58e-01

GTTGGGCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAGGCGGCGCAGGCCT  >  NC_002947/1219012‑1219301
                                                                                                                                                |                                                                                                                                                 
gTTGGGCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGgc                                                                                                                                              >  2:910641/1‑150 (MQ=255)
     gCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGGTCCGGCGCCgg                                                                                                                                         <  1:641405/150‑1 (MQ=255)
       gTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTg                                                                                                                                       <  1:828130/150‑1 (MQ=255)
            ccAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGtt                                                                                                                                  >  2:628293/1‑150 (MQ=255)
                tttCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGTACCCGGCGCTGGTGCAGTTGCTg                                                                                                                              <  1:119224/150‑1 (MQ=255)
                tttCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTg                                                                                                                              >  1:682730/1‑150 (MQ=255)
                tttCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTg                                                                                                                              >  2:243106/1‑150 (MQ=255)
                     tGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGGTCCGgcgccggga                                                                                                                                       >  2:659225/1‑131 (MQ=255)
                     tGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGGTCCGGCGCCggga                                                                                                                                       <  1:659225/136‑3 (MQ=255)
                      gTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGGTGCTGGCCCgg                                                                                                                        >  1:848578/1‑150 (MQ=255)
                         ttGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCt                                                                                                                               <  2:668323/140‑1 (MQ=255)
                         ttGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCt                                                                                                                               >  1:668323/1‑140 (MQ=255)
                          tGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGgcgca                                                                                                                    <  1:1182287/150‑1 (MQ=255)
                              tgtgAGCATGGGCGGCATGATTGCCCANCACCTGGCGGCGNTGGCGCNCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCGCGACCCGGCGCTGGTGCAGTTGCTGGCCCGgcgcagcgc                                                                                                                >  1:969057/1‑150 (MQ=255)
                                      tGGGCGGCATGATTGCCAAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAAcc                                                                                                        >  1:447856/1‑150 (MQ=255)
                                        ggcggcATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGt                                                                                                      <  2:848578/150‑1 (MQ=255)
                                             cATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGt                                                                                                 >  2:87928/1‑150 (MQ=255)
                                             cATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGCGAAAAGACTGGCCTCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGt                                                                                                 >  2:635075/1‑150 (MQ=255)
                                              aTGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTg                                                                                                <  1:164812/150‑1 (MQ=255)
                                                  ttGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCaa                                                                                            >  2:255158/1‑150 (MQ=255)
                                                        aGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGTCTTCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGagc                                                                                      <  2:1002589/150‑1 (MQ=255)
                                                        aGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGt                                                                                                 >  2:825711/1‑139 (MQ=255)
                                                        aGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGt                                                                                                 <  1:825711/139‑1 (MQ=255)
                                                          cACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGagcag                                                                                    >  1:648918/1‑150 (MQ=255)
                                                           aCCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGagcagc                                                                                   <  2:447856/150‑1 (MQ=255)
                                                              tGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAgg                                                                                >  1:477644/1‑150 (MQ=255)
                                                                  ggcgATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGa                                                                            <  1:316594/150‑1 (MQ=255)
                                                                  ggcgATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGa                                                                            >  2:1088607/1‑150 (MQ=255)
                                                                   gcgATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGAc                                                                           <  1:1061724/150‑1 (MQ=255)
                                                                     gATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACct                                                                         >  1:233556/1‑150 (MQ=255)
                                                                         gcgcCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACctgctg                                                                     <  2:1129382/150‑1 (MQ=255)
                                                                            cccGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTggcg                                                                  <  1:628293/150‑1 (MQ=255)
                                                                              cGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGAGGCGGGGCTCCCGGCGCCGGACCCGGCGGTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGGAGGCCGACATGCTggcggc                                                                >  2:870030/1‑150 (MQ=255)
                                                                                  cGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTg                                                            <  1:910641/150‑1 (MQ=255)
                                                                                   gtgtGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCa                                                                                             >  2:170234/1‑116 (MQ=255)
                                                                                   gtgtGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCa                                                                                             <  1:170234/116‑1 (MQ=255)
                                                                                    tgtgCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTggg                                                          >  2:912506/1‑150 (MQ=255)
                                                                                             ccTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGa                                                 <  2:477644/150‑1 (MQ=255)
                                                                                             ccTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGa                                                 >  1:294052/1‑150 (MQ=255)
                                                                                               tGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAgg                                               >  2:964380/1‑150 (MQ=255)
                                                                                               tGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGATCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTg                                                            >  2:1051077/1‑137 (MQ=255)
                                                                                               tGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGATCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTg                                                            <  1:1051077/137‑1 (MQ=255)
                                                                                                 aCGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTa                                             <  1:160011/150‑1 (MQ=255)
                                                                                                     tGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACggg                                         >  2:112578/1‑150 (MQ=255)
                                                                                                      ggTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGga                                        <  1:964380/150‑1 (MQ=255)
                                                                                                        tCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGagcag                                                                                    <  2:194902/104‑1 (MQ=255)
                                                                                                        tCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGagcag                                                                                    >  1:194902/1‑104 (MQ=255)
                                                                                                        tCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGgacg                                      >  2:531751/1‑150 (MQ=255)
                                                                                                            gTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGATCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACcg                                  >  1:463321/1‑150 (MQ=255)
                                                                                                            gTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACcg                                  >  1:34120/1‑150 (MQ=255)
                                                                                                             tCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGCACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGGCGACcgc                                 >  1:282393/1‑150 (MQ=255)
                                                                                                               cAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGa                               <  2:165832/150‑1 (MQ=255)
                                                                                                                  ttCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGgt                            <  1:255158/150‑1 (MQ=255)
                                                                                                                   tCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGgtg                           <  2:682730/150‑1 (MQ=255)
                                                                                                                      ggggCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTtg                        >  1:752130/1‑150 (MQ=255)
                                                                                                                         gcggcgGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTtgctg                     >  2:333439/1‑150 (MQ=255)
                                                                                                                            gcggGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCAc                  <  1:87928/150‑1 (MQ=255)
                                                                                                                                 gCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGCGTTGCTGCACCAggc             >  2:596586/1‑150 (MQ=255)
                                                                                                                                  cTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAggcg            <  2:294052/150‑1 (MQ=255)
                                                                                                                                    accagcgcCGGGTCCGGCGCCGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAgg                                                                                >  2:88528/5‑80 (MQ=38)
                                                                                                                                    accagcgcCGGGTCCGGCGCCGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAgg                                                                                <  1:88528/76‑1 (MQ=38)
                                                                                                                                     ccGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCCGgcggcg         >  2:708724/1‑150 (MQ=255)
                                                                                                                                     ccGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAGgcggcg         >  1:1130806/1‑150 (MQ=255)
                                                                                                                                      cGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACg                                           <  1:922118/115‑1 (MQ=255)
                                                                                                                                      cGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACg                                           >  2:922118/1‑115 (MQ=255)
                                                                                                                                        gcgcCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAGGCGGCGCAg      <  2:1155825/150‑1 (MQ=255)
                                                                                                                                        gcgcCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAGGCGGCGCAg      <  1:531751/150‑1 (MQ=255)
                                                                                                                                          gcCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAGGCGGCGCAGGc    <  2:752130/150‑1 (MQ=255)
                                                                                                                                          gcCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCGGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAGGCGGCGCAGGc    <  1:1088607/150‑1 (MQ=255)
                                                                                                                                            cGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAGGCGGCGCAGGCCt  <  2:969057/150‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                                 
GTTGGGCGTGCGCCAGTTTCATGTATTGGGTGTGAGCATGGGCGGCATGATTGCCCAGCACCTGGCGGCGATGGCGCCCGAGCGTGTGCGCAGCCTGACGCTGGTCATGTCCAGTTCAGGGGCGGCGGGGCTCCCGGCGCCGGACCCGGCGCTGGTGCAGTTGCTGGCCCGGCGCAGCGCACCGAACCGTGAAGTGGCAATCGAGCAGCAGGCCGACCTGCTGGCGGCCTTGGGCAGCCCTGAGGTACGGGACGACCGCGAGGTGTTGCTGCACCAGGCGGCGCAGGCCT  >  NC_002947/1219012‑1219301

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: