Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,079,615 C→T 6.0% intergenic (‑285/+29) cmtB ← / ← tktA putative mannitol‑specific enzyme IIA component of PTS/transketolase 1, thiamine triphosphate‑binding

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,079,6150CT 6.0% 241.2 / 6.5 102intergenic (‑285/+29)cmtB/tktAputative mannitol‑specific enzyme IIA component of PTS/transketolase 1, thiamine triphosphate‑binding
Reads supporting (aligned to +/- strand):  ref base C (41/52);  new base T (2/4);  total (45/57)
Fisher's exact test for biased strand distribution p-value = 6.95e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.60e-01

Only 100 of 104 total aligned reads displayed.

ATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAG  >  NC_000913/3079469‑3079763
                                                                                                                                                  |                                                                                                                                                    
ataGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGa                                                                                                                                                    >  1:846585/1‑149 (MQ=255)
  aGATATAACAACAAATGTCAATACGCATTTCGTGATGCATATTACATAATTGATATTTATTGCATTTAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAg                                                                                                                                                  <  1:662691/149‑1 (MQ=255)
    atatAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCg                                                                                                                                                >  1:345911/1‑149 (MQ=255)
       taACAACAAATGTCAATACGCATATCTTGATGCATATTACATTATTGATATTTATTGCATATAAATATTTGTGATCTAAAACACGCATTATCTATTGCTTGTCCGCAAACGGACATATCAAGGAAATAAAAAAGTTAGCCGAAGCGAcc                                                                                                                                             >  2:751313/1‑149 (MQ=255)
        aacaacAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGTTAATAAAAAAGTTCTCAGAAGAGACCt                                                                                                                                            >  2:594484/1‑149 (MQ=255)
         acaacaAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTAATGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAAAGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCtt                                                                                                                                           >  2:480919/1‑149 (MQ=255)
          caacaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCTTCGGCGACCttt                                                                                                                                          <  2:733753/149‑1 (MQ=255)
          caacaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCttt                                                                                                                                          >  2:935022/1‑149 (MQ=255)
           aacaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCTTCGGCGACCtttt                                                                                                                                         <  1:511729/149‑1 (MQ=255)
           aacaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCtttt                                                                                                                                         <  2:38295/149‑1 (MQ=255)
           aacaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCAAAGCGACCtttt                                                                                                                                         <  1:867239/149‑1 (MQ=255)
            acaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTCTTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCTTCGGCGACCttttt                                                                                                                                        <  1:788884/149‑1 (MQ=255)
            acaaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCttttt                                                                                                                                        <  2:846585/149‑1 (MQ=255)
             caaATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCTTCGGCGACCtttttt                                                                                                                                       <  1:1005140/149‑2 (MQ=255)
             caaATGTCAATACGCATATCGTGATCCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTa                                                                                                                                       >  2:175593/1‑149 (MQ=255)
              aaaTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCCAACGGACATATCAAGGTAATAAAAAAGGGCGCCGAAGCGACCTTTTTAc                                                                                                                                      <  1:496376/149‑1 (MQ=255)
              aaaTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCATTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCTCCGAAGCGACCTTTTTAc                                                                                                                                      >  1:954251/1‑149 (MQ=255)
                aTGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATCGCATATAAATATTTGTGATATACAACACGCCTTATCTATTGCTTGTCCGCAAACGGAAATATCAAGGTAATAAAAAAGGTCGCCGAAGTGACCTTTTTAcac                                                                                                                                    >  1:77190/1‑147 (MQ=255)
                    cAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGaaa                                                                                                                                <  1:158537/149‑1 (MQ=255)
                         cGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTa                                                                                                                           <  1:220532/149‑1 (MQ=255)
                            atatCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGCCGCCGAAGCGACCTTTTTACCCGAAATGCTAAtt                                                                                                                        <  1:476774/149‑1 (MQ=255)
                            aaatCGTGATGCATATTAAATAATTGATATTTATTGAATATAAATATTTGTGATCTACAACACGCATAATATATTCCTTGTCAGCAAACGGACATAAAAAGGTAATAAAAAAGGTCGCCGAAGCTACCTTATTACAAGAAATTCTAAta                                                                                                                        >  2:380708/3‑148 (MQ=255)
                              atagttaTGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACCACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTAc                                                                                                                      <  1:488181/143‑1 (MQ=255)
                              atCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACCCGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAACAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTAc                                                                                                                      <  2:960170/149‑1 (MQ=255)
                              atCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTAc                                                                                                                      <  1:1011010/149‑1 (MQ=255)
                               tCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTAca                                                                                                                     >  2:80734/1‑149 (MQ=255)
                                cGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTAcag                                                                                                                    <  2:345911/149‑1 (MQ=255)
                                  tGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTAcagca                                                                                                                  <  2:659172/149‑1 (MQ=255)
                                  tGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGAGACCTTTTTACCCGAAATGCTAATTAcagca                                                                                                                  >  1:189446/1‑149 (MQ=255)
                                   gATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTAcagcag                                                                                                                 <  1:642181/149‑1 (MQ=255)
                                      gCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTc                                                                                                              >  1:905553/1‑149 (MQ=255)
                                       cATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACACGAAATGCTAATTACAGCAGTTCt                                                                                                             >  2:644978/1‑149 (MQ=255)
                                        atatTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCCTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCtt                                                                                                            <  1:58369/149‑1 (MQ=255)
                                          atTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCtttt                                                                                                          <  1:512059/149‑1 (MQ=255)
                                          atTACATAATAGATATTTATTGCATATAAATATTTGTGATCTACAACACTCCTTATCTATTGCTTGTCCGCAAACGGAAATATCAAGGTAATAAAAAAGGGCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCtttt                                                                                                          <  2:175903/149‑1 (MQ=255)
                                             aCATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATAAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCt                                                                                                       >  1:115787/1‑149 (MQ=255)
                                                tAATTGATAATTATTGCATATAAATTTTTGTGATCTACAACACGCCTTATCTATTGCTTGGCCGCAAACGGACATATCAAGGTAATAAAAAAGGGCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTc                                                                                                    <  2:1015140/149‑1 (MQ=255)
                                                 aaTTGCTATTTATTGCATATAAATAATTGTTATCTACAAAACGCCTTATCTTTTGCTTGTCCGCATACGGACATATCAAGGTAATACAAAAGGTCTCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTct                                                                                                   <  1:575350/149‑2 (MQ=255)
                                                 aaTTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCg                                                                                                   <  1:575348/149‑1 (MQ=255)
                                                  aTTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGc                                                                                                  >  1:822832/1‑149 (MQ=255)
                                                    tGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCTCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGcaa                                                                                                <  1:17766/149‑1 (MQ=255)
                                                              tgaatatAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATc                                                                                      <  1:302139/146‑1 (MQ=255)
                                                              tGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATc                                                                                      <  1:889170/149‑1 (MQ=255)
                                                              tGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATc                                                                                      <  2:34863/149‑1 (MQ=255)
                                                               gCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATAACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCa                                                                                     >  2:518391/1‑149 (MQ=255)
                                                                  tataAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACa                                                                                  >  1:929626/1‑149 (MQ=255)
                                                                   ataAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATATAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAg                                                                                 <  2:817785/149‑1 (MQ=255)
                                                                    taAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGt                                                                                <  1:735570/149‑1 (MQ=255)
                                                                     aaaTATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTg                                                                               <  1:935022/149‑1 (MQ=255)
                                                                     aaaTATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTg                                                                               >  1:55908/1‑149 (MQ=255)
                                                                       atatTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGaa                                                                             >  2:20044/1‑149 (MQ=255)
                                                                             gtgATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACTTATCAAGGTAATAAAAAAGGTCGCCGCAGCGACCTGAGTACCCGAAATGCTAAGTACCGCCGTGCTTTTGCTTACGCAACAACGTTATCAACAGTGATGCCGac                                                                       >  1:575862/1‑148 (MQ=255)
                                                                              tgATCTACAACACACCTTATCTATTGCTTGTCCTCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCTACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAAc                                                                      <  2:444577/149‑1 (MQ=255)
                                                                                 tCTACAACTCGCTTTATGTTTTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAATGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTAGAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAActct                                                                   <  2:954251/149‑1 (MQ=255)
                                                                                 tCTACAACAAGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCAGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAActct                                                                   >  2:512190/1‑149 (MQ=255)
                                                                                  cTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCtt                                                                  <  2:951417/149‑1 (MQ=255)
                                                                                  cTACAACACGCCTTATCTAATGCTTGTCCGCAAACGGACATAACAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTAACCGAAATGCTAATTACAGCACTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCtt                                                                  >  1:706125/1‑149 (MQ=255)
                                                                                   ttcaacaCGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTc                                                                 <  2:464053/147‑1 (MQ=255)
                                                                                   tACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTc                                                                 <  2:11602/149‑1 (MQ=255)
                                                                                   tACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTc                                                                 <  2:799249/149‑1 (MQ=255)
                                                                                    acaacaCGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCa                                                                <  2:142640/149‑1 (MQ=255)
                                                                                     caacaCGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCaa                                                               >  1:154693/1‑149 (MQ=255)
                                                                                      aacaCGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCaaa                                                              >  1:810866/1‑149 (MQ=255)
                                                                                       acacGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAAc                                                             >  2:88778/1‑149 (MQ=255)
                                                                                         acGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAAcag                                                           <  1:69491/149‑1 (MQ=255)
                                                                                            ccTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAAcagcag                                                        <  2:651818/149‑1 (MQ=255)
                                                                                             ctgaTCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCAGAAGCGACCTTTTTACACGAATTGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGCACTCTTCAAACagcagc                                                       >  2:280788/4‑149 (MQ=255)
                                                                                                 tCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTg                                                   >  1:650207/1‑149 (MQ=255)
                                                                                                    aTTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCg                                                >  1:674680/1‑149 (MQ=255)
                                                                                                    aTTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCg                                                >  2:1005403/1‑149 (MQ=255)
                                                                                                        cTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCTGTTGTTTTGCTTTCGCAACTACGTTATCAACAGTGAAGGCGAACTCTTCAAACAGCAGCTCTGCCGGAGc                                            >  1:655059/1‑149 (MQ=255)
                                                                                                        cTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCAGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAgt                                            <  2:185022/149‑2 (MQ=255)
                                                                                                         ttGTCCGCATACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCa                                           <  1:589824/149‑1 (MQ=255)
                                                                                                            tCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGACCTCTTCAAACAGCAGCTCTGCCGGAGCAGAt                                        >  2:207401/1‑149 (MQ=255)
                                                                                                               gCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCa                                     >  2:331857/1‑149 (MQ=255)
                                                                                                               gCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCa                                     >  1:718086/1‑149 (MQ=255)
                                                                                                                cAAACGGACATATCAAGGTATTAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTATTTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCAc                                    <  2:164348/149‑1 (MQ=255)
                                                                                                                cAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCAc                                    <  1:175593/149‑1 (MQ=255)
                                                                                                                 aaaCGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCAcc                                   <  2:292363/149‑1 (MQ=255)
                                                                                                                 aaaCGGACATATCAAGGTAATAAAAAAGGACGCCGAAGCGACCGTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCACCATTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCCGATTCAcc                                   >  1:913167/1‑149 (MQ=255)
                                                                                                                       aCATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTCCCCGAAATGCTAATTACAGCAGTTGTTTTGCTTTCGCACCCACGTTATCAACAGTGAAGCCGAACTCTTGAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAggt                             >  1:676605/1‑149 (MQ=255)
                                                                                                                       aCATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAggt                             <  2:822832/149‑1 (MQ=255)
                                                                                                                        cATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAggtg                            <  2:674680/149‑1 (MQ=255)
                                                                                                                        aatatCAAGGTAATAAATTAGGTCGCCGAAGTGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAggtg                            <  2:77190/148‑1 (MQ=255)
                                                                                                                          tatCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAggtggt                          >  1:503538/1‑149 (MQ=255)
                                                                                                                          tatCAAGGTAATAAAAAAGGGCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAggtggt                          <  1:80734/149‑1 (MQ=255)
                                                                                                                           atCAAGTAAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCGAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCTGATTCACCGAAggtggta                         >  2:389913/1‑148 (MQ=255)
                                                                                                                           atCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTc                         >  2:512224/1‑149 (MQ=255)
                                                                                                                           atCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTAAAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTc                         <  2:87830/149‑1 (MQ=255)
                                                                                                                            tCAAGGTAATAAAAAAGGTCGCCGAAGCGCCCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCa                        >  2:887476/1‑149 (MQ=255)
                                                                                                                             cAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGct                                                      <  1:861854/118‑1 (MQ=255)
                                                                                                                             cAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGct                                                      >  2:861854/1‑118 (MQ=255)
                                                                                                                                    ggtaaaaaGGTCGCTTCGGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACcgacg                >  1:545789/4‑149 (MQ=255)
                                                                                                                                      taaaaaGGTCGCTTCGGCGACCGTTTTACCCGAAATGCTACTTACAGCAGTTCTTTTGCTTTCGCACCAACGGTATCAACAGTGAAGACGCACTCTTCAAACAGCAGCTCTGCCGGAGCAGAGTCCCCGCAGGTGGTCATACCGACGAt              >  1:472222/2‑149 (MQ=255)
                                                                                                                                      aaaaaaTGTCGACGAAGCGACCTTTTTACCCGAAATGCTAAATACAGCAGTTCTTTTGCTATCGCAACTACGTTATCAACAGTCAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGAt              <  1:751313/149‑1 (MQ=255)
                                                                                                                                      aaaaaaGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGAt              >  2:19350/1‑149 (MQ=255)
                                                                                                                                          aaGGTCGCCGAATCGACCTTTTTACCAGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGATGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCa          <  2:585696/149‑1 (MQ=255)
                                                                                                                                               cGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGtt     <  2:905553/149‑1 (MQ=255)
                                                                                                                                                 ccGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAATCAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCa   >  2:813699/1‑149 (MQ=255)
                                                                                                                                                  cGAAGCGACCTTTTTACCCGAATTGCTAATTACAGCAGTTCTTTTGCTTTCGCAACATCGTTATCAACAGTGAAGCCGAACTCTTCAATCAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAg  <  2:706125/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ATAGATATAACAACAAATGTCAATACGCATATCGTGATGCATATTACATAATTGATATTTATTGCATATAAATATTTGTGATCTACAACACGCCTTATCTATTGCTTGTCCGCAAACGGACATATCAAGGTAATAAAAAAGGTCGCCGAAGCGACCTTTTTACCCGAAATGCTAATTACAGCAGTTCTTTTGCTTTCGCAACAACGTTATCAACAGTGAAGCCGAACTCTTCAAACAGCAGCTCTGCCGGAGCAGATTCACCGAAGGTGGTCATACCGACGATAGCACCGTTCAG  >  NC_000913/3079469‑3079763

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