Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 2,124,684 T→C 5.5% I355V (ATT→GTT)  wcaI ← putative glycosyl transferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009132,124,6840TC 5.5% 229.2 / 2.4 73I355V (ATT→GTT) wcaIputative glycosyl transferase
Reads supporting (aligned to +/- strand):  ref base T (33/36);  new base C (2/2);  total (35/38)
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

TCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCG  >  NC_000913/2124536‑2124832
                                                                                                                                                    |                                                                                                                                                    
tCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGTCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAAt                                                                                                                                                      >  4:35519/1‑149 (MQ=255)
tCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAAt                                                                                                                                                      >  6:231516/1‑149 (MQ=255)
  aTTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGc                                                                                                                                                    >  4:240497/1‑149 (MQ=255)
     tataAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGccc                                                                                                                                                  <  8:372840/148‑1 (MQ=255)
       taAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCgg                                                                                                                                                >  6:199495/1‑148 (MQ=255)
         aaTTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGaaa                                                                                                                                             <  5:231516/149‑1 (MQ=255)
          aTTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAg                                                                                                                                            >  2:155374/1‑149 (MQ=255)
             gACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGGGTTGTGTTTGGTCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGt                                                                                                                                          <  3:35519/148‑1 (MQ=255)
               cGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTc                                                                                                                                       <  1:174612/149‑1 (MQ=255)
               cGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTc                                                                                                                                       <  3:240497/149‑1 (MQ=255)
                gTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCg                                                                                                                                      >  3:47148/1‑149 (MQ=255)
                 tAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCg                                                                                                                                      >  4:134267/1‑148 (MQ=255)
                  aaCACGTTCTCTTTATCGCGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTTGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGc                                                                                                                                     >  7:89276/1‑148 (MQ=255)
                     acGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAg                                                                                                                                 >  5:106329/1‑149 (MQ=255)
                             tttATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGTTTCAACGCAAACCGCAACGCCCGGAAAGGTTTCGCAAAGCTGccc                                                                                                                           >  5:347485/1‑147 (MQ=255)
                             tttATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGccc                                                                                                                           >  7:314738/1‑147 (MQ=255)
                             tttATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAACGCCCGGAAAGGTTTCGCAAAGCTGccc                                                                                                                           >  5:353013/1‑147 (MQ=255)
                                 tCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTc                                                                                                                     >  1:186873/1‑149 (MQ=255)
                                   gAGCGTGCGTTCGGCATATTCAGGTGCCACCGTGGTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCAGCAGCGCCTCGACCGATTCCGGTTCAACGCTAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCt                                                                                                                    <  5:29486/148‑1 (MQ=255)
                                    aGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCtgt                                                                                                                  <  2:26681/149‑1 (MQ=255)
                                       gtgcgtTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTa                                                                                                               <  8:314738/149‑1 (MQ=255)
                                                gCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTcagc                                                                                                      <  6:106329/149‑1 (MQ=255)
                                                 cATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTcagca                                                                                                     <  8:90519/149‑1 (MQ=255)
                                                  atatTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTcagcag                                                                                                    <  5:199495/149‑1 (MQ=255)
                                                     ttCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTa                                                                                                  <  3:134267/148‑1 (MQ=255)
                                                            gCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACcgc                                                                                          <  1:147642/149‑1 (MQ=255)
                                                               aCCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGt                                                                                        >  5:44689/1‑148 (MQ=255)
                                                               aCCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGt                                                                                        >  8:19644/1‑148 (MQ=255)
                                                                 cGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAACGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTAc                                                                                     <  6:353013/149‑1 (MQ=255)
                                                                 cGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAACGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTAc                                                                                     <  6:347485/149‑1 (MQ=255)
                                                                    gttgtgttTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGc                                                                                  >  1:156251/1‑149 (MQ=255)
                                                                        tgttTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTAc                                                                              >  2:382222/1‑149 (MQ=255)
                                                                        tgttTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGACCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGACTAc                                                                              >  2:249531/1‑149 (MQ=255)
                                                                                    agGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGTTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGattatt                                                                  <  2:293331/149‑1 (MQ=255)
                                                                                    agGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATAtt                                                                  <  7:19644/149‑1 (MQ=255)
                                                                                     ggAGCGCCTGACGGATCCCCTCCAACAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTTGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTg                                                                 <  1:249531/149‑1 (MQ=255)
                                                                                        gcgcCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATAtt                                                                  <  1:13947/145‑1 (MQ=255)
                                                                                        gcgcCTGAAGGATCCCCGCCACCAGCTCCTCGACCGATTCCGGATCAACGCAAACCGCAATGCCCGGAAAGGATTCGCAAAGAAGCCCCAGTTCAGTGTAGGCTACAGCAGTAATAACCGCGATACCGCCTACTGCCAGAATATTgag                                                               >  2:13947/1‑146 (MQ=255)
                                                                                              gACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCTGTTCTGTGTAGGCTTCGGCAGTAATCACCGCGTTAACGACTAGTGCCTGGATATTGGTCAGTTTc                                                        >  6:256218/1‑149 (MQ=255)
                                                                                              gACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATc                                                                                               >  5:281435/1‑110 (MQ=255)
                                                                                              gACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATc                                                                                               <  6:281435/110‑1 (MQ=255)
                                                                                               aCGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTc                                                        <  2:156251/148‑1 (MQ=255)
                                                                                                cGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGa                                                      <  4:47148/149‑1 (MQ=255)
                                                                                                   aTCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGTCAGAATATTGGTCAGTTTc                                                        <  1:112871/144‑1 (MQ=255)
                                                                                                   aTCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGTCAGAATATTGGTCAGTTTc                                                        >  2:112871/1‑144 (MQ=255)
                                                                                                    tCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACgg                                                   >  7:123843/1‑148 (MQ=255)
                                                                                                      cccGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCaa                                                >  4:316663/1‑149 (MQ=255)
                                                                                                      cccGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCaa                                                >  3:35667/1‑149 (MQ=255)
                                                                                                      cccGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCaa                                                >  7:110362/1‑149 (MQ=255)
                                                                                                             ccaGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCa                                         <  3:316663/149‑1 (MQ=255)
                                                                                                              caGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCAt                                        >  3:246227/1‑149 (MQ=255)
                                                                                                              caGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCAt                                        <  8:123843/149‑1 (MQ=255)
                                                                                                               aGCGCCTCGACCGATTACGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGATTATTGGTCAGTTTCGACGGCAATACGGCATc                                       <  8:89276/149‑1 (MQ=255)
                                                                                                                    cTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTgccg                                  >  2:332984/1‑149 (MQ=255)
                                                                                                                     tCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACg                                                    <  3:302117/130‑1 (MQ=255)
                                                                                                                     tCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACg                                                    >  4:302117/1‑130 (MQ=255)
                                                                                                                             ttCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGtt                         <  4:35667/149‑1 (MQ=255)
                                                                                                                              tCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGttt                        >  7:385219/1‑149 (MQ=255)
                                                                                                                                  gTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGCCGGCAATACGGCATCTGCCGCGCCGCGTTTTTGc                    <  1:332984/149‑1 (MQ=255)
                                                                                                                                  gTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGc                    <  6:44689/149‑1 (MQ=255)
                                                                                                                                  gTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGc                    <  8:110362/149‑1 (MQ=255)
                                                                                                                                   ttCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTACTCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGca                   >  4:295370/1‑149 (MQ=255)
                                                                                                                                      aaCGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTCTGcacca                >  7:53003/1‑149 (MQ=255)
                                                                                                                                      aaCGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTCTGcacca                >  7:58978/1‑149 (MQ=255)
                                                                                                                                       aCGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGcaccac               >  5:376521/1‑149 (MQ=255)
                                                                                                                                       aCGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGcaccac               >  3:192463/1‑149 (MQ=255)
                                                                                                                                         gCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTgtgt                                                                                                                >  7:205316/1‑50 (MQ=255)
                                                                                                                                         gCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTgtgt                                                                                                                <  8:205316/50‑1 (MQ=255)
                                                                                                                                             aCCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATg         <  6:376521/149‑1 (MQ=255)
                                                                                                                                             aCCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATg         <  2:186873/149‑1 (MQ=255)
                                                                                                                                             aCCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATg         <  1:155374/149‑1 (MQ=255)
                                                                                                                                              ccGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCTCGCCGCGTTTCTGCACCACCAGATgg        <  8:58978/149‑1 (MQ=255)
                                                                                                                                              ccGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTCTGCACCACCAGATgg        <  8:53003/149‑1 (MQ=255)
                                                                                                                                               cGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGc                    <  2:14557/136‑1 (MQ=255)
                                                                                                                                               cGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGc                    >  1:14557/1‑136 (MQ=255)
                                                                                                                                                    tGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGCCGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCg  >  2:244503/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                    
TCATTTATAAATTGACGTAACACGTTCTCTTTATCGAGCGTGCGTTCGGCATATTCACGTGCCACCGTGTTGTGTTTGGGCAGCAGGAGCGCCTGACGGATCCCCGCCACCAGCGCCTCGACCGATTCCGGTTCAACGCAAACCGCAATGCCCGGAAAGGTTTCGCAAAGCTGCCCCAGTTCTGTGTAGGCTTCAGCAGTAATCACCGCGTTACCGCCTACTGCCAGAATATTGGTCAGTTTCGACGGCAATACGGCATCTGCCGCGCCGCGTTTTTGCACCACCAGATGGCAATCG  >  NC_000913/2124536‑2124832

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