Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 910,542 C→T 64.8% E170K (GAA→AAA)  poxB ← pyruvate dehydrogenase, thiamine triphosphate‑binding, FAD‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_000913910,5420CT64.8% 106.4 / 104.1 108E170K (GAA→AAA) poxBpyruvate dehydrogenase, thiamine triphosphate‑binding, FAD‑binding
Reads supporting (aligned to +/- strand):  ref base C (19/19);  new base T (34/36);  total (53/55)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

Only 100 of 108 total aligned reads displayed.

CTCTTTATGCGCCCCCGCGCAGCCGCTGCCACACATCAGGGCGATATTGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTCTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCG  >  NC_000913/910396‑910686
                                                                                                                                                  |                                                                                                                                                
ctctTTATGCGCCCCCGCGCAGCCGCTGCCACACATCAGGGCGATATTGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTCTg                                                                                                                                                <  1:247405/149‑1 (MQ=255)
 tctTTATGCGCCCCCGCGCAGCCGCTGCCACACATCAGGGCGATATTGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTgg                                                                                                                                               >  5:331251/1‑149 (MQ=255)
   tttATGCGCCCCCGCGCAGCCGCTGCCACACATCAGGGCGATATTGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGgcg                                                                                                                                             >  5:365085/1‑149 (MQ=255)
    ttATGCGCCCCCGCGCAGCCGCTGCCACACATCAGGGCGATATTGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGgcgc                                                                                                                                            <  6:103466/149‑1 (MQ=255)
       tGCGCCCCCGCGCAGCCGCTGCCACACATCAGGGCGATATTGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAgg                                                                                                                                         >  2:328853/1‑149 (MQ=255)
       tGCGCCCCCGCGCAGCCGCTGCCACACATCAGGGCGATATTGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAgg                                                                                                                                         <  7:96845/149‑1 (MQ=255)
       tGCGCCCCCGCGCAGCCGCTGCCACACATCAGGGCGATATTGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAg                                                                                                                                          <  6:365085/148‑1 (MQ=255)
          gcCCCCGCGCAGCCGCTGCCACACATCAGGGCGATATTGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGttt                                                                                                                                      >  8:371757/1‑149 (MQ=255)
          gcCCCCGCGCAGCCGCTGCCACACATCAGGGCGATATTGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGttt                                                                                                                                      >  4:145778/1‑149 (MQ=255)
           cccccTCGCAGCAGCTGCCACACATCAGGGCGATCTTGCTGGAATAACGCAGCAGTTGCGCCAGTTTGGGGAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGGGCATGGTTGCCCCTTTTGGCGCAGGtttt                                                                                                                                     <  2:65695/149‑1 (MQ=255)
           cccccGCGCAGCCGCTGCCACACATCAGGGCGATATTGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGtttt                                                                                                                                     <  7:221758/149‑1 (MQ=255)
           ccacctcgcAGCAGCTGCCACCAATCAGGGCGATATTGCTGTAATAACGCAGAAGTTGCGACAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATCCCCGTGCATGTGTGCCCCTTTTGGCGCAGTtttt                                                                                                                                     <  2:59496/143‑1 (MQ=255)
            ccccGCGCAGCCGCTGCCACACATCAGGGCGATATTGCTGGAATACCGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTCTGGCGCAGGtttt                                                                                                                                     >  3:272483/1‑148 (MQ=255)
              ccGCGCAGCCGCTGCCACACATCAGGGCGATATTGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAAc                                                                                                                                  <  5:67929/149‑1 (MQ=255)
              ccGCGCAGCCGCTGCCACACATCAGGGCGATATTGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAAc                                                                                                                                  <  8:282588/149‑1 (MQ=255)
                 cgcAGCCGCTGCCACACATCAGGGCGATATTGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGcc                                                                                                                               <  6:93308/149‑1 (MQ=255)
                    aGCCGCTGCCACACATCAGGGCGATATTGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTCTGGCGCAGGTTTTAAc                                                                                                                                  <  6:310162/143‑1 (MQ=255)
                    aGCCGCTGCCACACATCAGGGCGATATTGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTCTGGCGCAGGTTTTAAc                                                                                                                                  >  5:310162/1‑143 (MQ=255)
                           gCCACACATCAGGGCGATATTGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTg                                                                                                                     >  6:130522/1‑149 (MQ=255)
                            ccAGACATCAGGGCGATATTGCTGGATAAACGCTGCAGGTGCGCGAGTTTGCGTAACTCTTCTTCTTCCTGCGGCACGACTTGTTGTGGCGCTTGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACTTCGCCt                                                                                                                      <  2:87368/147‑1 (MQ=255)
                            ccACACATCAGGGCGATATTGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCt                                                                                                                      <  6:331251/147‑1 (MQ=255)
                            ccACACATCAGGGCGATATTGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTgg                                                                                                                    >  1:115422/1‑149 (MQ=255)
                              acacaTCAGGGCGATATTGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGt                                                                                                                   >  3:345853/1‑148 (MQ=255)
                              acacaTCAGGGCGATATTGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTCTGGCGCAGGTTTTAACGCCACGTCGCCTGGt                                                                                                                   >  5:71808/1‑148 (MQ=255)
                                  aTCAGGCCGATATTGATGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTATTCTTCTTCCGGCGTCACGACTGGTTGTTGCGCATGATACCAGTGCATGGTTGCCCCTTCTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAacac                                                                                                              >  1:412982/1‑149 (MQ=255)
                                     aGGGCGATATTGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTCTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGAc                                                                                                           <  6:71808/149‑1 (MQ=255)
                                      gggCGATATTGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTCTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGAc                                                                                                           >  7:398787/1‑148 (MQ=255)
                                       ggCGATATTGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACaa                                                                                                         >  3:258826/1‑149 (MQ=255)
                                         cGATATTGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAAcc                                                                                                       >  6:179175/1‑149 (MQ=255)
                                            tatTGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGa                                                                                                     >  7:292270/1‑148 (MQ=255)
                                             atTGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTgg                                                                                                                                               >  6:339368/1‑105 (MQ=255)
                                             atTGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTgg                                                                                                                                               <  5:339368/105‑1 (MQ=255)
                                              ttGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGaaa                                                                                                   <  6:366595/148‑1 (MQ=255)
                                              ttGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGaaa                                                                                                   >  5:366595/1‑148 (MQ=255)
                                              ttGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAAc                                                                                                  <  1:328853/149‑1 (MQ=255)
                                                 cTGGAATAAAGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTTGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCCTGTTTTACCGCCACGTCGCCTGGTAACACGACAACCGAAACGc                                                                                                >  6:243874/1‑148 (MQ=255)
                                                    gAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCa                                                                                              <  7:395148/147‑1 (MQ=255)
                                                    gAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTCTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACg                                                                                            <  5:161205/149‑1 (MQ=255)
                                                     aaTAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGt                                                                                          >  4:82835/1‑149 (MQ=255)
                                                     aaTAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTCTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACgg                                                                                           <  1:375766/149‑1 (MQ=255)
                                                         aCGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTCTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTaa                                                                                       <  7:91344/149‑1 (MQ=255)
                                                          cGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTCTGGCGCAGGTTTTAAc                                                                                                                                  <  5:367157/105‑1 (MQ=255)
                                                          cGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTCTGGCGCAGGTTTTAAc                                                                                                                                  >  6:367157/1‑105 (MQ=255)
                                                            cagcagTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGTTACCAGTGCATGGTTGCCCCTTCTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCAc                                                                                             >  8:150334/1‑140 (MQ=255)
                                                            cagcagTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGTTACCAGTGCATGGTTGCCCCTTCTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCAc                                                                                             <  7:150334/140‑1 (MQ=255)
                                                            cagcagTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCa                                                                                    <  5:179175/149‑1 (MQ=255)
                                                              gcagTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCCTCAAGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCATGtt                                                                                                                                       >  1:301481/1‑96 (MQ=255)
                                                                  ttGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTCTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCtt                                                                              <  5:74444/149‑1 (MQ=255)
                                                                    gcgcCAGTTTGCGTAACTCTTCTTCTTCCGGCGGCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAGCGCCCCGTGGCCTGGTAACACGACAACCGACACGCCACGGTTAAGCACCGCTATg                                                                            >  2:303977/1‑149 (MQ=255)
                                                                      gcCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTCTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTgcg                                                                          >  5:56394/1‑149 (MQ=255)
                                                                       ccAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTGGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTgcgc                                                                         <  8:370759/149‑1 (MQ=255)
                                                                        cAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTCTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCa                                                                        <  4:272483/149‑1 (MQ=255)
                                                                         aGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATg                                                                      <  3:82835/149‑1 (MQ=255)
                                                                           tttGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTg                                                                                                                                                <  8:287227/74‑1 (MQ=255)
                                                                           tttGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTg                                                                                                                                                >  7:287227/1‑74 (MQ=255)
                                                                               cGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAAt                                                                 >  3:149913/1‑149 (MQ=255)
                                                                                gTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATc                                                                <  2:115422/149‑1 (MQ=255)
                                                                                   aCTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGc                                                              <  5:168862/148‑1 (MQ=255)
                                                                                    ctctTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGc                                                              >  3:267434/1‑147 (MQ=255)
                                                                                     tcttcttctTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTCTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAg                                                           >  2:328148/1‑149 (MQ=255)
                                                                                       ttcttcttcCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTCTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTgc                                                                           >  6:231395/1‑131 (MQ=255)
                                                                                         cttcttcCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAAc                                                                                                                                  <  4:43578/74‑1 (MQ=255)
                                                                                         cttcttcCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAAc                                                                                                                                  >  3:43578/1‑74 (MQ=255)
                                                                                         cttcttcCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTCTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCAcc                                                                                  <  2:290971/122‑1 (MQ=255)
                                                                                         cttcttcCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTCTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCAcc                                                                                  >  1:290971/1‑122 (MQ=255)
                                                                                            cttcCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTtgt                                                    >  4:156996/1‑149 (MQ=255)
                                                                                             ttcCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATg                                                                      >  5:22548/1‑130 (MQ=255)
                                                                                             ttcCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATg                                                                      <  6:22548/130‑1 (MQ=255)
                                                                                               ccGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTtgtg                                                   <  7:371757/147‑1 (MQ=255)
                                                                                               ccGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTtgt                                                    <  8:292270/146‑1 (MQ=255)
                                                                                                 ggCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATg                                                                      >  8:279143/1‑126 (MQ=255)
                                                                                                 ggCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATg                                                                      <  7:279143/126‑1 (MQ=255)
                                                                                                 ggCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTCTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTtgtg                                                   >  1:152987/1‑145 (MQ=255)
                                                                                                 ggCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTCTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTtgtg                                                   <  2:152987/145‑1 (MQ=255)
                                                                                                 ggCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTCTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGa                                                >  2:205093/1‑148 (MQ=255)
                                                                                                       aCGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTc                                         >  4:402388/1‑149 (MQ=255)
                                                                                                         gACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTCTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGa                                                                                                            <  2:55148/80‑1 (MQ=255)
                                                                                                         gACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTCTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGa                                                                                                            >  1:55148/1‑80 (MQ=255)
                                                                                                         gACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTCTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCg                                       >  6:12918/1‑149 (MQ=255)
                                                                                                            tGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCg                                       <  7:15417/146‑1 (MQ=255)
                                                                                                            tGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCg                                       >  8:15417/1‑146 (MQ=255)
                                                                                                              gTTGTGGCGGATGATACCAGTTCATGGTTGCCCCTTCTGTCGCAGGTTTTAACGCCACGTCGCCTTGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATATCCAGTACTTTTGGGATCTGCTCCGGGCt                                   <  2:412982/148‑1 (MQ=255)
                                                                                                              gTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGc                                    <  3:145778/147‑1 (MQ=255)
                                                                                                                  tgGCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGaaa                              <  1:125291/149‑1 (MQ=255)
                                                                                                                   ggCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAAc                             <  4:303495/149‑1 (MQ=255)
                                                                                                                   ggCGCATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAAc                             >  2:264778/1‑149 (MQ=255)
                                                                                                                    gcgcATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAAcc                            <  4:149913/149‑1 (MQ=255)
                                                                                                                      gcATGATACCAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAg                          >  3:48527/1‑149 (MQ=255)
                                                                                                                      gcATGATACCAGTGCATGGTTGCCCCTTCTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAg                          >  7:305441/1‑149 (MQ=255)
                                                                                                                             aCCAGTGCATGGTTGCCCCTTCTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCaa                   >  3:200121/1‑149 (MQ=255)
                                                                                                                               cAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACCCGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCaa                   >  6:336541/1‑147 (MQ=255)
                                                                                                                               cAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCt                        >  7:237547/1‑142 (MQ=255)
                                                                                                                               cAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATa                 <  4:267434/149‑1 (MQ=255)
                                                                                                                               cAGTGCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATa                 <  6:223916/149‑1 (MQ=255)
                                                                                                                                   gCATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTATGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGa             <  7:83758/149‑1 (MQ=255)
                                                                                                                                    cATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGc                         >  4:396190/1‑136 (MQ=255)
                                                                                                                                    cATGGTTGCCCCTTTTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGc                         <  3:396190/136‑1 (MQ=255)
                                                                                                                                    cATGGTTGCCCCTTCTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGAc            <  4:200121/149‑1 (MQ=255)
                                                                                                                                         ttGCCCCTTCTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACAt       <  8:305441/149‑1 (MQ=255)
                                                                                                                                              ccTTCTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTcgcg  <  5:12918/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
CTCTTTATGCGCCCCCGCGCAGCCGCTGCCACACATCAGGGCGATATTGCTGGAATAACGCAGCAGTTGCGCCAGTTTGCGTAACTCTTCTTCTTCCGGCGTCACGACTGGTTGTGGCGCATGATACCAGTGCATGGTTGCCCCTTCTGGCGCAGGTTTTAACGCCACGTCGCCTGGTAACACGACAACCGAAACGCCACGGTTAAGCACCGCTTTGCGCATGGCAATCGCCAGTACTTGTGGGATCTGCTCCGGGCTGGAAACCAGCTCGCAATAGTGACTACATTCGCG  >  NC_000913/910396‑910686

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