Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,129,512 T→C 5.1% T47A (ACC→GCC)  wcaE ← putative glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,129,5120TC 5.1% 239.1 / 2.8 79T47A (ACC→GCC) wcaEputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base T (33/42);  new base C (2/2);  total (35/44)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 2.30e-01

AAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCGACAGTGATTATGCTAAGCAACATGTTCTTATT  >  NC_000913/2129367‑2129659
                                                                                                                                                 |                                                                                                                                                   
aaaaTATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGcc                                                                                                                                                  <  8:152719/149‑1 (MQ=255)
      tCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGtt                                                                                                                                            <  1:330345/149‑1 (MQ=255)
        gCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTgg                                                                                                                                          <  5:254707/149‑1 (MQ=255)
        gCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTg                                                                                                                                           <  3:217644/148‑1 (MQ=255)
         cccGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTgg                                                                                                                                          >  8:231804/1‑148 (MQ=255)
         cccGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTgg                                                                                                                                          >  4:253146/1‑148 (MQ=255)
         cccGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCAGTCGTTgg                                                                                                                                          >  6:314830/1‑148 (MQ=255)
          ccGACTTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGaa                                                                                                                                        <  2:5234/149‑1 (MQ=255)
          ccGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGaa                                                                                                                                        >  4:189437/1‑149 (MQ=255)
              aTTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGc                                                                                                                                    <  7:231804/149‑1 (MQ=255)
                  agaAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCAt                                                                                                                                 <  5:190925/148‑1 (MQ=255)
                    aaaCAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGa                                                                                                                              >  8:123067/1‑149 (MQ=255)
                    aaaCAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGa                                                                                                                              <  7:15510/149‑1 (MQ=255)
                     aacaacGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGac                                                                                                                             <  4:190243/149‑1 (MQ=255)
                      acaacGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGac                                                                                                                             >  7:101076/1‑148 (MQ=255)
                             gAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATc                                                                                                                      >  4:169362/1‑148 (MQ=255)
                              aaCTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCa                                                                                                                    <  4:346940/149‑1 (MQ=255)
                              aaCTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCa                                                                                                                    <  3:253146/149‑1 (MQ=255)
                               aCTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCAt                                                                                                                   <  6:211561/149‑1 (MQ=255)
                                   gCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTc                                                                                                                 <  1:307887/147‑1 (MQ=255)
                                         tgCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTg                                                                                                          <  5:314830/148‑1 (MQ=255)
                                           ccATCGCAATGCGTTTGTTCATGGCGTCGTAGATACCGTTATGTGGGTCGCTTACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGata                                                                                                       <  8:173414/149‑1 (MQ=255)
                                                gCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTc                                                                                                  <  3:169362/149‑1 (MQ=255)
                                                gCAATGCCTTTGTTCAGGGCGCCGTAGATACCGGTATCTGGCTCGCTGACAAAGCGTAGGTTAACGATACCATGGAGCTTTTCCAGATACGCACGGGTGCCGTCGTGGGAACCGCCATCGACAACAATCCCTGCGACGCTGATATCtt                                                                                                   <  2:48686/148‑1 (MQ=255)
                                                 cAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTc                                                                                                  >  2:41251/1‑148 (MQ=255)
                                                  aaTGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTCACAAAGCGTAGGTTAAAGATAACATTGAGATTTTCCAGATACTCACGGGTGCCGACGTTGGAACCGCCATCGACAACAATCCAATCGAAGCTGATATCTTCCa                                                                                                >  1:338134/1‑149 (MQ=255)
                                                    tGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCAcc                                                                                              <  4:119590/149‑1 (MQ=255)
                                                        tttGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTgcg                                                                                          >  6:79531/1‑149 (MQ=255)
                                                           gTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTgcgc                                                                                         >  4:297538/1‑147 (MQ=255)
                                                               aTGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGTCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGAt                                                                                    <  3:189437/148‑1 (MQ=255)
                                                                tGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATgc                                                                                  <  3:275978/149‑1 (MQ=255)
                                                                 ggCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATgcg                                                                                 >  8:39116/1‑149 (MQ=255)
                                                                    gtcgtAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCa                                                                              <  6:332450/149‑1 (MQ=255)
                                                                       gtAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATgcgc                                                                                <  8:149027/144‑1 (MQ=255)
                                                                       gtAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATgcgc                                                                                >  7:149027/1‑144 (MQ=255)
                                                                          gATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAgg                                                                        >  7:268225/1‑149 (MQ=255)
                                                                             aCCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGCGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCAt                                                                     >  2:58603/1‑149 (MQ=255)
                                                                             aCCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGCGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCAt                                                                     >  2:64387/1‑149 (MQ=255)
                                                                              ccGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCAt                                                                     <  3:297538/148‑1 (MQ=255)
                                                                               cGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCAtgt                                                                   >  8:58869/1‑149 (MQ=255)
                                                                                   aTCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGttt                                                               >  5:189768/1‑149 (MQ=255)
                                                                                    tCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGttt                                                               <  7:123067/148‑1 (MQ=255)
                                                                                      tGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGACGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGa                                                            >  2:272040/1‑149 (MQ=255)
                                                                                        gCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGCCAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGtttt                                                              <  5:79531/145‑1 (MQ=255)
                                                                                        gCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGa                                                            <  8:268225/147‑1 (MQ=255)
                                                                                        gCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGAc                                                           >  6:230542/1‑148 (MQ=255)
                                                                                          tCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTAt                                                        >  4:293531/1‑149 (MQ=255)
                                                                                            gCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGcc                                                                                                                                                  >  6:61652/1‑57 (MQ=255)
                                                                                            gCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGcc                                                                                                                                                  <  5:61652/57‑1 (MQ=255)
                                                                                              tGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCt                                                    >  2:353420/1‑149 (MQ=255)
                                                                                                   aaGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCAc                                                                                               >  7:403257/1‑101 (MQ=255)
                                                                                                   aaGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCAc                                                                                               <  8:403257/101‑1 (MQ=255)
                                                                                                   aaGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAg                                               <  6:189768/149‑1 (MQ=255)
                                                                                                       gTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTa                                           <  7:39116/149‑1 (MQ=255)
                                                                                                          ggTTAAAGATACCATTGAGATTTTCCAGATACTGACTGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACg                                         <  3:206667/148‑1 (MQ=255)
                                                                                                          ggTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACg                                         >  4:2286/1‑148 (MQ=255)
                                                                                                              aaaGATACCATTGAGATTTTCCAGATACTCACGGGCGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAAc                                     <  1:64387/148‑1 (MQ=255)
                                                                                                              aaaGATACCATTGAGATTTTCCAGATACTCACGGGCGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAAc                                     <  1:58603/148‑1 (MQ=255)
                                                                                                               aaGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAAcgc                                   >  5:42049/1‑149 (MQ=255)
                                                                                                                 gATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCGa                                 >  6:283842/1‑149 (MQ=255)
                                                                                                                    aCCATTGAGATTTTCCAGATACTCACGGGTGCGGGCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCGACa                               <  2:338134/148‑1 (MQ=255)
                                                                                                                     ccATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTAGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCGACAGt                             >  5:163260/1‑149 (MQ=255)
                                                                                                                      cATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGc                                                                            <  8:133751/101‑1 (MQ=255)
                                                                                                                      cATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGc                                                                            >  7:133751/1‑101 (MQ=255)
                                                                                                                       aTTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCGACa                               <  1:300516/145‑1 (MQ=255)
                                                                                                                       aTTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCGACa                               >  2:300516/1‑145 (MQ=255)
                                                                                                                       aTTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCGACAGTGa                           <  8:101076/149‑1 (MQ=255)
                                                                                                                         tGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCGACAGTGAt                          >  2:59343/1‑148 (MQ=255)
                                                                                                                           tgaTTTTACAGATACTCACGGGTGGCGTCGTTGGAACCGCCATCGACTACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGGGTTTTGACTATCCCTTCGAGGTTACGAAACGCGACAGTGATTa                        <  5:219000/147‑1 (MQ=255)
                                                                                                                           agaTTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCGACAGTGATTa                        <  1:41251/148‑1 (MQ=255)
                                                                                                                           agaTTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCGACAGTGATTAt                       >  3:100140/1‑149 (MQ=255)
                                                                                                                                 tCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCGACAGTGATTATGCTaa                  <  5:365345/148‑1 (MQ=255)
                                                                                                                                  ccAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCGACAGTGATTATGCTAAg                 >  1:31163/1‑148 (MQ=255)
                                                                                                                                   ctgATACTCACGGGTTCCGTCGTTTTAACCGCCATAGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAg                                                                         <  8:75809/89‑1 (MQ=255)
                                                                                                                                   cAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAg                                                                         >  7:75809/1‑91 (MQ=255)
                                                                                                                                     gATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCGACAGTGATTATGCTAAGCaa              >  6:76717/1‑148 (MQ=255)
                                                                                                                                            aCGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCGACAGTGATTATGCTAAGCAACATGttc       <  4:9892/148‑1 (MQ=255)
                                                                                                                                                gTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCGACAGTGATTATGCTAAGCAACATGttcttatt  <  6:163260/149‑1 (MQ=255)
                                                                                                                                                gTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCGACAGTGATTATGCTAAGCAACATGttcttatt  <  3:293531/149‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                                   
AAAATATCGCCCGAATTGAGAAACAACGCGAACTTGCCTTGTGCCATCGCAATGCCTTTGTTCATGGCGTCGTAGATACCGTTATCTGGCTCGCTGACAAAGCGTAGGTTAAAGATACCATTGAGATTTTCCAGATACTCACGGGTGCCGTCGTTGGAACCGCCATCGACAACAATCCATTCGAAGCTGATATCTTCCACCTGCGCCAGATGCGCCAGCGAGGCATGTGTTTTGACTATCCCTTCGAGGTTACGAAACGCGACAGTGATTATGCTAAGCAACATGTTCTTATT  >  NC_000913/2129367‑2129659

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: