Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_002947 3,614,722 A→C 5.3% T199P (ACC→CCC)  pet → protein Pet18

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0029473,614,7220AC 5.3% 237.6 / 2.4 85T199P (ACC→CCC) petprotein Pet18
Reads supporting (aligned to +/- strand):  ref base A (42/38);  new base C (3/2);  total (45/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 = 9.48e-01

GAAACCGATGCCGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGGCGCAG  >  NC_002947/3614574‑3614865
                                                                                                                                                      |                                                                                                                                                
gAAACCGATGCCGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCa                                                                                                                                                  >  1:92828/1‑149 (MQ=255)
 aaaCCGATGCCGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCAc                                                                                                                                                 >  1:194344/1‑149 (MQ=255)
  aaCCGATGCCGGACAACTTCGTCCATTCCGAGTGGATTACCCTGCATGACAACCCTGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGACGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACAAGCGCAcc                                                                                                                                                >  7:346128/1‑149 (MQ=255)
    ccGATGCCGTACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAAACCCGCATTCAAGGCATTCGTCGCTTTCCTGCAGTCAGAACAGGACCGCGACCCACCTCAGGAGGACGCTTTTGACAGTGACTTCAACCACCGCACCGt                                                                                                                                              >  6:233284/1‑149 (MQ=255)
      gATGCCGGAGTACTTCGTGCATGCCGTGTGGATTACCCTGCATGACAAGGCGGCATTCGAGGCATTCGTCGCTTTCCTGCAGGCCGATCTGGACCTCGTCGGCCCGCAGGAGGCCGCTTTGGGAAGTGACTTCTACCAGCGCAGCGTgg                                                                                                                                            <  2:340009/149‑1 (MQ=255)
      gATGCCGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTgg                                                                                                                                            >  3:96/1‑149 (MQ=255)
           cGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTg                                                                                                                                       >  3:33325/1‑149 (MQ=255)
            ggACTACTTCGTCCATGCCGAGTGGATTAACCTGCATGACAACCCTGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCTAACTGGAACGCGTCGGCCCGCAGGATGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTgg                                                                                                                                      <  8:346128/149‑1 (MQ=255)
                aaCTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCt                                                                                                                                  >  5:342553/1‑149 (MQ=255)
                  cTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCCACGGTGCGCTggtagaag                                                                                                                                <  3:257/149‑7 (MQ=255)
                  cTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCCACGGTGCGCTggtagaag                                                                                                                                <  4:300766/149‑7 (MQ=255)
                  cTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTgg                                                                                                                                >  5:126572/1‑149 (MQ=255)
                         cATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAg                                                                                                                                    >  2:64461/1‑138 (MQ=255)
                         cATGCCGAGTGGAATAACCTTCATGCCAACCCGGCATTCCAGGCATTCGTGGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCATGATGGCGCTTTGTGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAg                                                                                                                                    <  1:64461/139‑1 (MQ=255)
                               gAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGc                                                                                                                  >  8:100336/1‑149 (MQ=255)
                                aGTGGAATACCCTGCATGACAACCCGGCATTCCAGGCAATCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAACTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCt                                                                                                                 >  2:116789/1‑149 (MQ=255)
                                 gTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCtg                                                                                                                <  3:290899/149‑1 (MQ=255)
                                 gTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCtg                                                                                                                <  6:66616/149‑1 (MQ=255)
                                  tGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCtgt                                                                                                               <  7:189400/149‑1 (MQ=255)
                                  tGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCtgt                                                                                                               <  6:193958/149‑1 (MQ=255)
                                   ggATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCtgtg                                                                                                              >  1:52309/1‑149 (MQ=255)
                                     aTTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTGTGGCAAGTGACTTCTACCAGCGCCCCGTGGCGCTGCAGCTGGCGTTCTTCGCTGCTGTTTa                                                                                                            >  6:368013/1‑149 (MQ=255)
                                         cccTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCg                                                                                                        >  1:329677/1‑149 (MQ=255)
                                          ccTGCATGACAACCCGGCATTCCAGGCATTCGTACGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCg                                                                                                        >  1:46021/1‑149 (MQ=255)
                                           cTGCATGTCAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCg                                                                                                        >  5:211090/1‑147 (MQ=255)
                                           cTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGa                                                                                                      <  4:33325/149‑1 (MQ=255)
                                              cATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAg                                                                                                     <  2:52309/147‑1 (MQ=255)
                                                tGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTgg                                                                                                                                            >  4:19058/1‑107 (MQ=255)
                                                tGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTgg                                                                                                                                            <  3:19058/107‑1 (MQ=255)
                                                  aCAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCaa                                                                                                <  2:92828/148‑1 (MQ=255)
                                                    aaCCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGGCGAACTGGACCGCGTGGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCg                                                                                                                      <  8:340455/124‑1 (MQ=255)
                                                    aaCCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCCCCGTGGCGCTGGCGCTGGCTTTCGTCGAGGCTGTGTACCCGGAGGGCAAAt                                                                                              >  8:341300/1‑148 (MQ=255)
                                                    aaCCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCg                                                                                                                      >  7:340455/1‑124 (MQ=255)
                                                    aaCCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATg                                                                                             <  2:46021/149‑1 (MQ=255)
                                                     aCCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGa                                                                                            <  5:33176/149‑1 (MQ=255)
                                                       ccGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGc                                                                                          <  2:339503/149‑1 (MQ=255)
                                                       ccGGCATTCCAGGCACTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAg                                                                                           >  8:222865/1‑148 (MQ=255)
                                                                cAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATc                                                                                 >  3:341641/1‑149 (MQ=255)
                                                                 aGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCt                                                                                >  5:182062/1‑149 (MQ=255)
                                                                   gCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTc                                                                              <  8:317291/149‑1 (MQ=255)
                                                                     aTTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGa                                                                            <  6:126572/149‑1 (MQ=255)
                                                                        cgtcgCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAgcg                                                                         >  7:297498/1‑149 (MQ=255)
                                                                        cgtcgCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAgc                                                                          <  4:140884/148‑1 (MQ=255)
                                                                          tcgCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGcc                                                                       >  7:270401/1‑149 (MQ=255)
                                                                          tcgCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGcc                                                                       >  4:331420/1‑149 (MQ=255)
                                                                          tcgCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGcc                                                                       <  7:100336/149‑1 (MQ=255)
                                                                               tttCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTc                                                                     >  7:138961/1‑147 (MQ=255)
                                                                                 tCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAg                                                                  <  6:211090/148‑1 (MQ=255)
                                                                                  ccTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGATCGCCTCAAgg                                                                 <  6:342553/148‑1 (MQ=255)
                                                                                      cAGGCCGCACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCg                                                              <  7:341300/147‑1 (MQ=255)
                                                                                       aGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTg                                                           <  8:138961/149‑1 (MQ=255)
                                                                                                tGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTgccacgcca                                                   >  4:144653/1‑148 (MQ=255)
                                                                                                tGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTgccacgcca                                                   >  3:351826/1‑148 (MQ=255)
                                                                                                  gACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCGTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACaa                                                <  4:96/149‑1 (MQ=255)
                                                                                                      gcgTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGa                                            <  6:182062/149‑1 (MQ=255)
                                                                                                        gTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAAc                                          <  2:329677/149‑1 (MQ=255)
                                                                                                             cccGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATGTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTa                                     >  2:216482/1‑149 (MQ=255)
                                                                                                              ccGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCACGGCCGCTGCCACGCCACAATGGAACCGCTAc                                    >  8:167549/1‑149 (MQ=255)
                                                                                                              ccGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTAc                                    >  6:142259/1‑149 (MQ=255)
                                                                                                              ccGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTAc                                    <  7:222865/149‑1 (MQ=255)
                                                                                                              ccGCAGGAAGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTAc                                    >  2:280175/1‑149 (MQ=255)
                                                                                                                 cAGGAGGCCGCTTTGTCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTa                                 <  6:133810/149‑1 (MQ=255)
                                                                                                                 cAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTa                                 >  7:20506/1‑149 (MQ=255)
                                                                                                                     aggCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGAcc                             <  4:351826/149‑1 (MQ=255)
                                                                                                                       gCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCAt                           <  2:194344/149‑1 (MQ=255)
                                                                                                                        ccGCTTTGGCATGTGACTTCTACCATCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTAACCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTAAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCAt                           <  5:233284/148‑1 (MQ=255)
                                                                                                                        ccGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCAt                           >  8:140772/1‑148 (MQ=255)
                                                                                                                              tGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTc                     <  4:341641/148‑1 (MQ=255)
                                                                                                                                 gCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTg                  >  5:383183/1‑149 (MQ=255)
                                                                                                                                  cAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGc                 <  1:280175/149‑1 (MQ=255)
                                                                                                                                   aaGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTgcca                                                        >  6:109018/1‑109 (MQ=255)
                                                                                                                                   aaGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTgcca                                                        <  5:109018/109‑1 (MQ=255)
                                                                                                                                   aaGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCg                <  1:190533/149‑1 (MQ=255)
                                                                                                                                   aaGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCg                >  6:238967/1‑149 (MQ=255)
                                                                                                                                     gTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCgg               >  7:29009/1‑148 (MQ=255)
                                                                                                                                     gTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCgg               >  3:59879/1‑148 (MQ=255)
                                                                                                                                      tGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGc              >  2:174889/1‑148 (MQ=255)
                                                                                                                                      tGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGc              <  1:174889/148‑1 (MQ=255)
                                                                                                                                         ccagctccAGCGCCACGGTGCGCTGGAGCTGGCTTTCTACGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAAAGGAACCGCTACGTAGACAATGACTTCGTGCGGCAGAt          >  2:301027/7‑149 (MQ=255)
                                                                                                                                         cTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATg                                                                                       <  4:361008/72‑1 (MQ=255)
                                                                                                                                         cTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATg                                                                                       >  3:361008/1‑72 (MQ=255)
                                                                                                                                          ttCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATg         >  7:311695/1‑149 (MQ=255)
                                                                                                                                          ttCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATg         >  3:54447/1‑149 (MQ=255)
                                                                                                                                           tCTACCAGCGCACCGTGGCGCTGGAGCTGGACTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATg         >  1:36577/1‑149 (MQ=255)
                                                                                                                                             tACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGgc      >  1:362030/1‑149 (MQ=255)
                                                                                                                                                 aGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGGCGCAg  <  2:362030/149‑1 (MQ=255)
                                                                                                                                                      |                                                                                                                                                
GAAACCGATGCCGGACAACTTCGTCCATGCCGAGTGGATTACCCTGCATGACAACCCGGCATTCCAGGCATTCGTCGCTTTCCTGCAGGCCGAACTGGACCGCGTCGGCCCGCAGGAGGCCGCTTTGGCAAGTGACTTCTACCAGCGCACCGTGGCGCTGGAGCTGGCTTTCTTCGATGCTGTGTACCCGGAGGGCAAATGAGCATGGATATCTTCGAGCGCCTCAAGGCCGCTGCCACGCCACAATGGAACCGCTACGTAGACCATGACTTCGTGCGGCAGATGGGCGCAG  >  NC_002947/3614574‑3614865

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