Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,698,587 T→A 5.5% A198A (GCA→GCT) ‡ malI ← transcriptional repressor of Mal regulon

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,698,5870TA 5.5% 229.7 / 3.6 73A198A (GCA→GCT) ‡malItranscriptional repressor of Mal regulon
Reads supporting (aligned to +/- strand):  ref base T (34/35);  new base A (2/2);  total (36/37)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.32e-01

GGCACTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACT  >  NC_000913/1698443‑1698721
                                                                                                                                                |                                                                                                                                      
ggCACTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCAc                                                                                                                                    <  1:15496/149‑1 (MQ=255)
  cACTGATCGTCGGGTTATGACGTAACAGCCCCGTCATAGCTTCCGCGGCTTGCTTCTCGCTGGCAGTGCCCTCCAACACCCAATCGCTATGACACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACgg                                                                                                                                  <  4:300085/149‑1 (MQ=255)
   aCTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACggg                                                                                                                                 >  7:402287/1‑149 (MQ=255)
        tGGTCTGGTTATGACGTAAAAGCGCCGTGATCGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTaa                                                                                                                             <  3:254898/148‑1 (MQ=255)
         ggTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAAt                                                                                                                            >  7:345688/1‑148 (MQ=255)
          gTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGa                                                                                                                          <  3:239295/149‑1 (MQ=255)
          gTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGa                                                                                                                          <  4:107700/149‑1 (MQ=255)
          gTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGa                                                                                                                          <  5:387975/149‑1 (MQ=255)
            cGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCATGCTTCTGGCTGGAAGATCACTCCAACACCCAATCGATATGAAACGCCAGGCCAAATTTTAGTAGAGTAGCACAAAAGCCCCCAACCCGTTCTGACCGGGTTAATGAg                                                                                                                         >  5:161406/1‑148 (MQ=255)
               gTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAAc                                                                                                                     >  7:225124/1‑149 (MQ=255)
                   tGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTg                                                                                                                 <  3:190631/149‑1 (MQ=255)
                     aCGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGcc                                                                                                               >  2:130894/1‑149 (MQ=255)
                      cGTTAAAGCGGCGTGATAGCTTCCGCGTCTTGCTTCTGTCTGGAATGGTACTCCTACACCCAATCGTTATGAAAGTGAAGGCCAAATTTTAGTAGAGTTGCACAATTCCCCACAAGCCGTTCTGCACTGGTTAATGTGTAACTTTGccc                                                                                                              <  2:394614/149‑1 (MQ=255)
                      cGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGccc                                                                                                              >  3:391917/1‑149 (MQ=255)
                      cGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGcc                                                                                                               >  4:10854/1‑148 (MQ=255)
                        tAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACCCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTc                                                                                                            >  6:162448/1‑149 (MQ=255)
                         aaaaGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTcc                                                                                                           >  2:334489/1‑149 (MQ=255)
                                ccGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCccagccag                                                                                                    >  6:163372/1‑149 (MQ=255)
                                        gCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCg                                                                                            <  1:130894/149‑1 (MQ=255)
                                         cTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCCACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGc                                                                                           >  7:141977/1‑149 (MQ=255)
                                         cTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGc                                                                                           <  8:402287/149‑1 (MQ=255)
                                            ccGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTg                                                                                         <  3:125634/148‑1 (MQ=255)
                                            ccGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGa                                                                                        >  2:174715/1‑149 (MQ=255)
                                              gcgGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGATACGGCAGGCCAATTTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATg                                                                                      >  5:2288/1‑149 (MQ=255)
                                                  cttgcttCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCa                                                                                  <  1:319853/149‑1 (MQ=255)
                                                      cttCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGc                                                                              >  1:409149/1‑149 (MQ=255)
                                                            gctggAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATga                                                                        <  5:163372/149‑1 (MQ=255)
                                                              tggAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAAt                                                                          >  4:191369/1‑145 (MQ=255)
                                                                gAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATgaga                                                                      <  8:225124/147‑1 (MQ=255)
                                                                gAAGTGCACTCCAACAACCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGAt                                                                     <  5:346151/148‑1 (MQ=255)
                                                                   gTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTc                                                                 <  1:53289/149‑1 (MQ=255)
                                                                      cACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGt                                                              >  1:396074/1‑149 (MQ=255)
                                                                      cACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGt                                                              <  1:174715/149‑1 (MQ=255)
                                                                        cTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGt                                                              <  7:210764/147‑1 (MQ=255)
                                                                        cTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGt                                                              >  8:210764/1‑147 (MQ=255)
                                                                         tCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTGCAGAACGGGttg                                                                                                                              <  2:17700/82‑2 (MQ=255)
                                                                         tCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTGCAGAACGGGttg                                                                                                                              >  1:17700/1‑81 (MQ=255)
                                                                             acacCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCATTGCGAATGAGATGCTCCGTCAACAACt                                                      <  8:141977/149‑1 (MQ=255)
                                                                               acCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACtg                                                     >  1:357333/1‑149 (MQ=255)
                                                                                  cactCGCTATGAAACGGCAGGCCTAATTTTAGTAGAGTTTCACAATAGCCCCCAACCCGTTATGCACGGGTTAATGAGGAACTTTGCCGTCCCCGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCTACAACTGTGc                                                  <  6:161406/146‑1 (MQ=255)
                                                                                  cAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGc                                                  >  5:300770/1‑149 (MQ=255)
                                                                                  cAATAGCTATGAAACTGCAGGCCAAATTTTAGTAGAGTTGCACAAAAGACCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCTACCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACtgtc                                                   >  8:168673/1‑147 (MQ=255)
                                                                                   aaTCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACtgtg                                                   >  3:332549/1‑147 (MQ=255)
                                                                                     tCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGc                                               <  3:10854/149‑1 (MQ=255)
                                                                                       gCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGcc                                              <  3:191369/148‑1 (MQ=255)
                                                                                                  gCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTcc                                  <  5:162448/149‑1 (MQ=255)
                                                                                                   cAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCg                                 >  8:318380/1‑149 (MQ=255)
                                                                                                    aGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCgg                                <  1:334489/149‑1 (MQ=255)
                                                                                                       ccAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCACCAACTGTGCAGCCTGCATGTTGTCCGGGc                              >  6:273028/1‑148 (MQ=255)
                                                                                                         aaaTTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGaa                           <  4:391917/149‑1 (MQ=255)
                                                                                                          aaTTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACTGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGc                                                                                     <  6:165747/90‑1 (MQ=255)
                                                                                                                      gTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGGCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAAcatc              <  6:2288/149‑1 (MQ=255)
                                                                                                                          cacaATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACatcatc           <  2:357333/148‑1 (MQ=255)
                                                                                                                          cacaATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACatcatc           >  5:294388/1‑148 (MQ=255)
                                                                                                                          cacaATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACatcatc           >  4:375507/1‑148 (MQ=255)
                                                                                                                          cacaATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACatcatc           >  6:173181/1‑148 (MQ=255)
                                                                                                                           acaATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCAAGTTGTCCGGGCGAACCGTATCAACATCATCaa         <  2:396074/149‑1 (MQ=255)
                                                                                                                            caATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATgag                                                                       <  6:357599/86‑1 (MQ=255)
                                                                                                                            caATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATgag                                                                       >  5:357599/1‑86 (MQ=255)
                                                                                                                             aaTAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACatcatc           >  1:308526/1‑145 (MQ=255)
                                                                                                                                  cccccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAAc                 <  6:32234/134‑1 (MQ=255)
                                                                                                                                  cccccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAAc                 >  5:32234/1‑134 (MQ=255)
                                                                                                                                  cccccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACt  >  6:350221/1‑149 (MQ=255)
                                                                                                                                  cccccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACt  >  8:275821/1‑149 (MQ=255)
                                                                                                                                  cccccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACt  <  2:409149/149‑1 (MQ=255)
                                                                                                                                    cccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTAt                     <  5:238465/128‑1 (MQ=255)
                                                                                                                                    cccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTAt                     >  6:238465/1‑128 (MQ=255)
                                                                                                                                     ttaaCCCGTGCAGAACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCAt                                         <  3:331841/105‑1 (MQ=255)
                                                                                                                                     ttaaCCCGTGCAGAACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCAt                                         >  4:331841/3‑107 (MQ=255)
                                                                                                                                     ccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATc                    <  6:301977/128‑1 (MQ=255)
                                                                                                                                     ccAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATc                    >  5:301977/1‑128 (MQ=255)
                                                                                                                                        aCCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGc                                               >  7:305452/1‑98 (MQ=255)
                                                                                                                                        aCCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGc                                               <  8:305452/98‑1 (MQ=255)
                                                                                                                                                |                                                                                                                                      
GGCACTGATGGTCGGGTTATGACGTAAAAGCGCCGTGATAGCTTCCGCGGCTTGCTTCTGGCTGGAAGTGCACTCCAACACCCAATCGCTATGAAACGGCAGGCCAAATTTTAGTAGAGTTGCACAATAGCCCCCAACCCGTTCTGCACGGGTTAATGAGGAACTTTGCCCTCCCAGCCAGGCGATCCGCTGATGCCCATTGCGAATGAGATGCTCCGTCAACAACTGTGCAGCCTGCATGTTGTCCGGGCGAACCGTATCAACATCATCAAGATAACT  >  NC_000913/1698443‑1698721

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: