Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,608,695 C→G 6.6% intergenic (+46/+56) zntA → / ← tusA zinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,608,6950CG 6.6% 224.5 / 4.0 68intergenic (+46/+56)zntA/tusAzinc, cobalt and lead efflux system/mnm(5)‑s(2)U34‑tRNA 2‑thiolation sulfurtransferase
Reads supporting (aligned to +/- strand):  ref base C (28/35);  new base G (2/3);  total (30/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

CACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAA  >  NC_000913/3608547‑3608842
                                                                                                                                                    |                                                                                                                                                   
cacGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCctc                                                                                                                                                     >  8:5986/1‑149 (MQ=255)
 acGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCctct                                                                                                                                                    >  8:387984/1‑149 (MQ=255)
  cGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTg                                                                                                                                                   <  6:272611/149‑1 (MQ=255)
       ttAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGga                                                                                                                                               <  3:212081/148‑1 (MQ=255)
         aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGgaga                                                                                                                                             <  2:348671/148‑1 (MQ=255)
         aGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGATAAGATTGTAGCGCAGGAGATAAGGCAAACCGATCGCTACCTTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGagag                                                                                                                                            >  3:11597/1‑149 (MQ=255)
              tcaCCGGGTTGTGGCTGGCAGTGCTGGCAGATCCGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTa                                                                                                                                       <  1:275341/147‑1 (MQ=255)
              tGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTa                                                                                                                                       <  3:90466/149‑1 (MQ=255)
              tGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTa                                                                                                                                       <  3:148187/149‑1 (MQ=255)
              tGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACAGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTa                                                                                                                                       <  6:190083/149‑1 (MQ=255)
              tGAACGGGTTGTGGCTGGCAGTGCTGGCAGATACTGGGGCGACGGTGCGTGTGACAGCGAATGCGTTAAGATTGTTGGGCAGGTGATTAGGCAAAGCGATCGCAACTTTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGTGTTa                                                                                                                                       <  1:230888/149‑1 (MQ=255)
               gACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTa                                                                                                                                       <  8:107121/148‑1 (MQ=255)
                acagggTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCTAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGCTCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGGTAgg                                                                                                                                     >  5:230818/4‑149 (MQ=255)
                aCCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                                                     <  7:78112/149‑1 (MQ=255)
                aCCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                                                     <  3:205512/149‑1 (MQ=255)
                  cGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGt                                                                                                                                   <  8:393619/149‑1 (MQ=255)
                   gggTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGt                                                                                                                                   >  7:340829/1‑148 (MQ=255)
                     gTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAg                                                                                                                                <  3:160627/149‑1 (MQ=255)
                      ttGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAgg                                                                                                                               >  1:240411/1‑149 (MQ=255)
                           gctggcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGAGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGAGAGGGGaaa                                                                                                                          <  1:347887/149‑1 (MQ=255)
                               gcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGTGTGAGGGGAAAAgg                                                                                                                       >  7:297403/1‑148 (MQ=255)
                               gcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggc                                                                                                                      <  8:52646/149‑1 (MQ=255)
                               gcAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAg                                                                                                                        <  8:297403/147‑1 (MQ=255)
                                cAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAggcg                                                                                                                     <  7:387984/149‑1 (MQ=255)
                                     cTGGCAGATACGGGGGCGACGGGGCTGGGGACAGCGAATGCGTTATGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTGTGGGGAGAGGGTTAGGGTGAGGGTAAAAGGCGGCATc                                                                                                                <  8:223434/149‑1 (MQ=39)
                                      tGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATc                                                                                                                >  2:48365/1‑148 (MQ=255)
                                        tcAGATACGGTGGCGACGGTGCTGGGGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCTCGATCGGTCCCCTGGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATc                                                                                                                <  4:199380/145‑1 (MQ=255)
                                           gATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGc                                                                                                           >  6:385070/1‑148 (MQ=255)
                                               cGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATc                                                                                                      >  2:234654/1‑149 (MQ=255)
                                                  gggCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGc                                                                                                   <  4:204067/149‑1 (MQ=255)
                                                    gCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGccc                                                                                                 <  2:240411/149‑1 (MQ=255)
                                                    gCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGccc                                                                                                 >  2:60466/1‑149 (MQ=255)
                                                    gCGACGGTGCTGGTCACAGCGAATGCGTTAAGAATGTAGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGGCCCCTCGCCCCTCTGGGGAGAGGGTTAGTGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGccc                                                                                                 <  7:304761/149‑1 (MQ=39)
                                                      gACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGcccc                                                                                                >  6:175496/1‑148 (MQ=255)
                                                          gTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGCAGATAAGGCAAACCGATCGCAACATTTAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATc                                                                                            >  8:158014/1‑148 (MQ=255)
                                                                 tGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCt                                                                                    >  1:154878/1‑149 (MQ=255)
                                                                 tGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCt                                                                                    >  7:105091/1‑149 (MQ=255)
                                                                   aCAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgggg                                                                                                                                                <  5:335611/87‑1 (MQ=255)
                                                                   aCAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTgggg                                                                                                                                                >  6:335611/1‑87 (MQ=255)
                                                                     aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                                                     >  1:6756/1‑96 (MQ=255)
                                                                     aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAgg                                                                                                                                     <  2:6756/96‑1 (MQ=255)
                                                                     aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTAc                                                                                <  5:385070/149‑1 (MQ=255)
                                                                     aGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTAc                                                                                <  8:340829/149‑1 (MQ=255)
                                                                         aaTGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCa                                                                                                                  <  1:305620/111‑1 (MQ=255)
                                                                         aaTGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCa                                                                                                                  >  2:305620/1‑111 (MQ=255)
                                                                              gTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCaaa                                                                        >  1:52675/1‑148 (MQ=255)
                                                                                  aGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCGGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAAc                                                                    >  5:163644/1‑148 (MQ=255)
                                                                                   gATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGa                                                                  >  7:30522/1‑149 (MQ=255)
                                                                                        ttgCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAgg                                                             <  1:90016/149‑1 (MQ=255)
                                                                                            gcAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGt                                                         >  1:340746/1‑149 (MQ=255)
                                                                                                 agaTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACTAATCAAATAACGAGAAGGCAGTCCATc                                                    >  4:344028/1‑149 (MQ=255)
                                                                                                   aTAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCGCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATcc                                                   <  8:355147/148‑1 (MQ=255)
                                                                                                     aaGGCACACCGATCGCAACTTTGAGCGCGTTCGGTCCCCTCGCCCCTCTGGGGTGAGGGTTAGGGTGTGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTa                                                                                 <  8:349016/116‑1 (MQ=38)
                                                                                                     aaGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCACAGGGGAGAGGCTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCACCCCCTATCAACCGCCTTTa                                                                                 >  7:349016/1‑116 (MQ=255)
                                                                                                      aGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCATCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGCAGCCAATCAGCCCATATCACCCGCCTTTACGAATCAAAAAACGATAAGGCAGTCCATCCGtct                                               >  7:237093/1‑149 (MQ=38)
                                                                                                         ctaaCCGTTCGCAACATTGAGGGCGATCGGTCCCCTCGCCCCTGTGGGGAGAGGGTTAGGGTGAGGGGAAAACGCGCCATCGAAGCCAATCAGCCCCTAACACCCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                            <  4:11597/147‑1 (MQ=38)
                                                                                                         cAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAGTCAAATAACGATAAGGCAGTCCATCCGTCTCtt                                            <  7:5986/149‑1 (MQ=255)
                                                                                                           aaCCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTa                                          <  2:55203/149‑1 (MQ=255)
                                                                                                                 tCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                     >  6:69713/1‑148 (MQ=255)
                                                                                                                         ttGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtg                            >  3:158749/1‑149 (MQ=255)
                                                                                                                         ttGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCgtg                            <  8:105091/149‑1 (MQ=255)
                                                                                                                                 gATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATa                     >  6:18360/1‑148 (MQ=255)
                                                                                                                                    cctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAgg                 >  6:389682/6‑149 (MQ=255)
                                                                                                                                     ctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtct                                               >  8:339285/5‑118 (MQ=255)
                                                                                                                                     ctaacccTCTCCCCAGAGGGGCGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGtct                                               <  7:339285/114‑1 (MQ=255)
                                                                                                                                         cccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                     <  7:154285/124‑1 (MQ=255)
                                                                                                                                         cccTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAAc                                     >  8:154285/1‑124 (MQ=255)
                                                                                                                                           cTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCACCCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACaaaa          <  2:340746/149‑1 (MQ=255)
                                                                                                                                               ccccTCTGTGGAGAGGGTTATGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCa      <  6:230818/149‑1 (MQ=255)
                                                                                                                                                 ccTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAg     >  1:294049/1‑148 (MQ=255)
                                                                                                                                                   tctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGaa  >  7:405602/1‑149 (MQ=255)
                                                                                                                                                   tctGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGaa  >  7:388912/1‑149 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
CACGCTGTTAGGGATGACCGGGTTGTGGCTGGCAGTGCTGGCAGATACGGGGGCGACGGTGCTGGTGACAGCGAATGCGTTAAGATTGTTGCGCAGGAGATAAGGCAAACCGATCGCAACATTGAGCGCGATCGGTCCCCTCGCCCCTCTGGGGAGAGGGTTAGGGTGAGGGGAAAAGGCGGCATCGAAGCCAATCAGCCCCTATCAACCGCCTTTACGAATCAAATAACGATAAGGCAGTCCATCCGTCTCTTTAGCAACCAGTTCGTGTTCCATAAAGGTACAAAACCCAGGAA  >  NC_000913/3608547‑3608842

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: