Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,520,487 T→A 5.8% T73S (ACT→TCT) ‡ hofQ ← DNA catabolic putative fimbrial transporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,520,4870TA 5.8% 245.1 / 3.7 80T73S (ACT→TCT) ‡hofQDNA catabolic putative fimbrial transporter
Reads supporting (aligned to +/- strand):  ref base T (35/40);  new base A (3/2);  total (38/42)
Fisher's exact test for biased strand distribution p-value = 6.64e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.59e-01

CTGCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCA  >  NC_000913/3520342‑3520627
                                                                                                                                                 |                                                                                                                                            
cTGCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTg                                                                                                                                           <  3:161811/149‑1 (MQ=255)
 tGCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTg                                                                                                                                           >  2:384466/1‑148 (MQ=255)
 tGCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTg                                                                                                                                           >  2:364058/1‑148 (MQ=255)
   cGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTa                                                                                                                                         <  5:323890/148‑1 (MQ=255)
   cGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTa                                                                                                                                         <  5:81204/148‑1 (MQ=255)
        ttCAAGCGGCAGATTTGCCTGCGCCCCCGCCTCCTCCCCCTCCTGGGGGGCGAAATTGATATTCTGCCAGGCAATCAAAACCACAGAGACAAAGTTGACTTCCTGCCGTGTTTTCAGTCCGGCGCTTTAAAATACAGAATGTAGAGCCt                                                                                                                                   >  8:58147/1‑149 (MQ=255)
         tCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTg                                                                                                                                  >  3:329716/1‑149 (MQ=255)
           cAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCTCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGAAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGc                                                                                                                                 <  7:1859/148‑1 (MQ=255)
           cAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCt                                                                                                                                >  6:245846/1‑149 (MQ=255)
             gCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGc                                                                                                                                     <  3:24635/142‑1 (MQ=255)
             gCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGc                                                                                                                                     >  4:24635/1‑142 (MQ=255)
             gCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTc                                                                                                                              <  4:329716/149‑1 (MQ=255)
               ggCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCa                                                                                                                            <  1:384466/149‑1 (MQ=255)
               ggCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCGCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATTCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCa                                                                                                                            <  2:309047/149‑1 (MQ=255)
                gCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAg                                                                                                                           >  6:321050/1‑149 (MQ=255)
                           tGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCa                                                                                                                            <  3:4442/137‑1 (MQ=255)
                           tGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGTAACATCTGt                                                                                                                >  6:141710/1‑149 (MQ=255)
                           tGCGCCCGCGCCTGCTCCGCCTCCTCGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCa                                                                                                                            >  4:4442/1‑137 (MQ=255)
                            gcgcCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGt                                                                                                                >  5:141408/1‑148 (MQ=255)
                             cgcCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGt                                                                                                                <  5:321050/147‑1 (MQ=255)
                             cgcCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTa                                                                                                              <  6:88851/149‑1 (MQ=255)
                                    gcCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTaa                                                                                                       <  1:364058/149‑1 (MQ=255)
                                         cTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGa                                                                                                    <  2:103069/147‑1 (MQ=255)
                                          tCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGc                                                                                                                                     <  3:118989/113‑1 (MQ=255)
                                          tCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGc                                                                                                                                     >  4:118989/1‑113 (MQ=255)
                                          tCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGacac                                                                                                 >  5:94731/1‑149 (MQ=255)
                                              cctcctGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTa                                                                                             <  1:377008/149‑1 (MQ=255)
                                                     gcgggcgATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGa                                                                                      <  6:220569/149‑1 (MQ=255)
                                                           gATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTACCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTg                                                                                >  3:84748/1‑149 (MQ=255)
                                                            atatTGTTATTCTGCCATGCAATTGAATGCACTGAGAGAAAGTTGCCTTCCTGCCGAGATATAACACCGGATCATTTCAATAAACTATGTAGAGCCAGCATCCAGGGAACATCTGTTACATGTAACGACACAGTACCGCAGACGACTgg                                                                               >  2:141085/1‑149 (MQ=255)
                                                            atatTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTgg                                                                               >  8:155430/1‑149 (MQ=255)
                                                             tatTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGAAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGc                                                                              >  8:364497/1‑149 (MQ=255)
                                                              atTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGACCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGc                                                                              >  7:226194/1‑148 (MQ=255)
                                                                tGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGac                                                                           >  3:263154/1‑149 (MQ=255)
                                                                  ttATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGacac                                                                         >  4:173741/1‑149 (MQ=255)
                                                                  ttATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGacac                                                                         >  6:125823/1‑149 (MQ=255)
                                                                     ttCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGAAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGAc                                                                      <  7:364497/149‑1 (MQ=255)
                                                                        tGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAAACTGTAGTgaaaa                                                                                                                                  >  1:58578/1‑82 (MQ=255)
                                                                        tGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAAACTGTAGTgaaaa                                                                                                                                  <  2:58578/86‑5 (MQ=255)
                                                                          ccAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAg                                                                   <  2:201311/147‑1 (MQ=255)
                                                                            aGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGtt                                                                >  5:215735/1‑148 (MQ=255)
                                                                            aGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTc                                                               <  5:245846/149‑1 (MQ=255)
                                                                              gCAATGGAAAGAACTGAGAGAATGTTTCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCaa                                                             <  1:313213/149‑1 (MQ=255)
                                                                               cAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTACCGACACCGTACCGCTGACGTCAGGCGCCACGACCAGGTTCAAc                                                            >  8:143718/1‑149 (MQ=255)
                                                                               cAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAAc                                                            <  5:141710/149‑1 (MQ=255)
                                                                                aaTGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGt                                                                 >  2:382593/1‑143 (MQ=255)
                                                                                aaTGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGt                                                                 <  1:382593/143‑1 (MQ=255)
                                                                                    gAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTtctc                                                       <  5:125823/149‑1 (MQ=255)
                                                                                          aCTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTg                                                                                <  8:295209/118‑1 (MQ=255)
                                                                                          aCTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTg                                                                                >  7:295209/1‑118 (MQ=255)
                                                                                          aCTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTc                                                 >  7:306415/1‑149 (MQ=255)
                                                                                          aCTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTc                                                 <  4:263154/149‑1 (MQ=255)
                                                                                          aCTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAAACTGTAGTgaa                                                                                                                                    >  2:353201/1‑64 (MQ=38)
                                                                                          aCTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAAACTGTAGTgaa                                                                                                                                    <  1:353201/66‑3 (MQ=38)
                                                                                           cTGAGTGATTGTTGCCTTCCTGCCGTGTTATCTGTCCTGGGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCATGGAACATCTTTTAGATGTAACGACACCGTAGCGCTGACGTGTGGGGACACGAGCAGGTTCTTGTTCTCCTGTTca                                                <  7:58147/149‑1 (MQ=255)
                                                                                            tGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCa                                                                                                                            <  4:47922/72‑1 (MQ=255)
                                                                                            tGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCa                                                                                                                            >  3:47922/1‑72 (MQ=255)
                                                                                              agagaATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTca                                                >  3:42399/1‑146 (MQ=255)
                                                                                              agagaATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTca                                                <  4:42399/146‑1 (MQ=255)
                                                                                               gagaATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTcagc                                              >  2:247649/1‑147 (MQ=255)
                                                                                                  aaTGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAgcg                                         >  2:34987/1‑149 (MQ=255)
                                                                                                        gCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGc                                    <  7:136121/148‑1 (MQ=255)
                                                                                                                gCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCt                             >  5:64858/1‑147 (MQ=255)
                                                                                                                 ccGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTg                            <  1:34987/147‑1 (MQ=255)
                                                                                                                 ccGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGa                           <  6:215735/148‑1 (MQ=255)
                                                                                                                 ccGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGa                           <  6:150948/148‑1 (MQ=255)
                                                                                                                  cgcgTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAAc                                                            <  1:330187/114‑1 (MQ=255)
                                                                                                                  cgcgTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAAc                                                            >  2:330187/1‑114 (MQ=255)
                                                                                                                     gTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTAc                      >  4:134688/1‑149 (MQ=255)
                                                                                                                          cAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGaa                 >  7:341170/1‑149 (MQ=255)
                                                                                                                           aGTCCGGCGCTTTTCTCTACATTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGCGCTACCGGaa                 <  8:226194/148‑1 (MQ=255)
                                                                                                                           aGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAAc                <  6:40317/149‑1 (MQ=255)
                                                                                                                                ggCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTTCCGGACCGTCAt           >  5:3134/1‑149 (MQ=255)
                                                                                                                                ggCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCAt           >  6:15985/1‑149 (MQ=255)
                                                                                                                                   gcTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGCCACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATcca        >  6:354769/1‑149 (MQ=255)
                                                                                                                                       ttCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATccacc      <  3:173741/147‑1 (MQ=255)
                                                                                                                                        tCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCAt    <  6:94731/148‑1 (MQ=255)
                                                                                                                                          aCTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCa  <  5:15985/148‑1 (MQ=255)
                                                                                                                                          aCTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCa  <  3:134688/148‑1 (MQ=255)
                                                                                                                                          aCTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCa  <  6:64858/148‑1 (MQ=255)
                                                                                                                                                 |                                                                                                                                            
CTGCGATTTTCCAGCGGCAGATTTGCCTGCGCCCGCGCCTGCTCCGCCTCCTGGCGGGCGATATTGTTATTCTGCCAGGCAATGGAATGCACTGAGAGAATGTTGCCTTCCTGCCGCGTTATCAGTCCGGCGCTTTTCACTACAGTTTGTAGTGCCTGCTTCCAGGGAACATCTGTTAGATGTAACGACACCGTACCGCTGACGTCTGGCGACACGACCAGGTTCAACTTCTCCTGTTCAGCCAGCGCCTGCAACACCTGAGCTACCGGAACGTCATCCACCATCA  >  NC_000913/3520342‑3520627

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