Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,568,919:1 +C 5.3% intergenic (‑96/+35) yddW ← / ← gadC liprotein, glycosyl hydrolase homolog/glutamate:gamma‑aminobutyric acid antiporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,568,9191.C 5.3% 603.0 / 28.1 133intergenic (‑96/+35)yddW/gadCliprotein, glycosyl hydrolase homolog/glutamate:gamma‑aminobutyric acid antiporter
Reads supporting (aligned to +/- strand):  ref base . (62/64);  new base C (4/3);  total (66/67)
Fisher's exact test for biased strand distribution p-value = 7.18e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 7.16e-01

Only 100 of 136 total aligned reads displayed.

ATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGT  >  NC_000913/1568774‑1569067
                                                                                                                                                  |                                                                                                                                                    
aTCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATc                                                                                                                                                   <  3:304761/149‑1 (MQ=255)
aTCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATc                                                                                                                                                   <  1:623715/149‑1 (MQ=255)
aTCGCTAGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATc                                                                                                                                                   <  3:538826/149‑1 (MQ=255)
 tCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGAGATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCTATGAATT‑ATCg                                                                                                                                                  <  1:936189/149‑1 (MQ=255)
       gTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCACTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCt                                                                                                                                            >  2:571483/1‑149 (MQ=255)
        tcttctAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCtt                                                                                                                                           >  2:583119/1‑149 (MQ=255)
           tctAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTc                                                                                                                                        <  4:177801/149‑1 (MQ=255)
                  gTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAAc                                                                                                                                 >  2:440345/1‑149 (MQ=255)
                   ttAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATc                                                                                                                                                   <  4:817829/130‑1 (MQ=255)
                   ttAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAAcc                                                                                                                                >  1:883411/1‑149 (MQ=255)
                   ttAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAAcc                                                                                                                                >  1:506030/1‑149 (MQ=255)
                       tttCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTCGTCTTATAACCATTc                                                                                                                            >  4:579111/1‑149 (MQ=255)
                       tttCTTGTTTCGGGAGCAGATATCCATAAGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGAAGAATT‑ATCGCTCCCTTGTCTTATAACCATTc                                                                                                                            >  1:51132/1‑149 (MQ=255)
                        ttCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCa                                                                                                                           <  4:411366/149‑1 (MQ=255)
                           ttGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTc                                                                                                                                               <  4:1084130/126‑1 (MQ=255)
                           ttGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGAc                                                                                                                        >  1:924973/1‑149 (MQ=255)
                               ttCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATgg                                                                                                                    <  1:447130/149‑1 (MQ=255)
                                tCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGt                                                                                                                   <  4:167735/149‑1 (MQ=255)
                                 cGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGtt                                                                                                                  <  3:1035508/149‑1 (MQ=255)
                                 cGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGtt                                                                                                                  <  3:1111464/149‑1 (MQ=255)
                                 cGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGtt                                                                                                                  <  3:659392/149‑1 (MQ=255)
                                  gggAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTa                                                                                                                 <  2:883411/149‑1 (MQ=255)
                                  gggAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTa                                                                                                                 >  1:895357/1‑149 (MQ=255)
                                     aGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAgtg                                                                                                              <  2:895357/149‑1 (MQ=255)
                                         gATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTc                                                                                                          >  3:902525/1‑149 (MQ=255)
                                         gATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCGCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTc                                                                                                          <  2:213549/149‑1 (MQ=255)
                                           tatCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCtt                                                                                                        >  2:878747/1‑149 (MQ=255)
                                             tCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGt                                                                                                      <  1:990664/149‑1 (MQ=255)
                                              ccATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGtc                                                                                                     <  1:583119/149‑1 (MQ=255)
                                                 tatGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGtcatt                                                                                                  >  1:935262/1‑149 (MQ=255)
                                                   tgtgCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGtcattca                                                                                                <  2:660312/149‑1 (MQ=255)
                                                    gtgCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGtcattcat                                                                                               <  3:208391/149‑1 (MQ=255)
                                                      gCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTccc                                                                                                                                             >  1:236668/1‑101 (MQ=255)
                                                      gCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTccc                                                                                                                                             <  2:236668/101‑1 (MQ=255)
                                                      gCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATtcatca                                                                                             <  1:244378/149‑1 (MQ=255)
                                                        tGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATcaca                                                                                           >  4:319416/1‑149 (MQ=255)
                                                        tGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGCCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAgtg                                                                                                              >  1:1012449/1‑130 (MQ=255)
                                                         ggTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACaa                                                                                          <  4:902525/149‑1 (MQ=255)
                                                          gTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAat                                                                                         <  1:440345/149‑1 (MQ=255)
                                                               cGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATAATTC‑ATCGCTcc                                                                                                                                              <  1:148619/91‑1 (MQ=39)
                                                               cGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGAT‑AATTCATCGCTcc                                                                                                                                              >  2:148619/1‑91 (MQ=39)
                                                               cGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGAT‑AATTCATCGCTc                                                                                                                                               <  1:1011832/90‑1 (MQ=255)
                                                               cGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGAT‑AATTCATCGCTc                                                                                                                                               >  2:1011832/1‑90 (MQ=255)
                                                                ggTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAAtata                                                                                      >  4:1119871/1‑146 (MQ=255)
                                                                ggTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAAtata                                                                                      <  3:1119871/146‑1 (MQ=255)
                                                                ggTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAgtg                                                                                   >  2:201998/1‑149 (MQ=255)
                                                                 gTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAgtgt                                                                                  <  1:878747/149‑1 (MQ=255)
                                                                 gTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACAACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAgtgt                                                                                  >  3:617417/1‑149 (MQ=255)
                                                                  tAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGtgtg                                                                                 <  3:319416/149‑1 (MQ=255)
                                                                    aaCAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTgtggt                                                                               >  3:864818/1‑149 (MQ=255)
                                                                          tgtgCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGt                                                                         <  3:216857/149‑1 (MQ=255)
                                                                          tgtgCACTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGt                                                                         <  1:571483/149‑1 (MQ=255)
                                                                            tgCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGc                                                                       <  3:322708/149‑1 (MQ=255)
                                                                                cTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCAcgc                                                                   >  3:741036/1‑149 (MQ=255)
                                                                                 tCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTgtgg                                                                                >  1:161713/1‑135 (MQ=255)
                                                                                 tCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTgtgg                                                                                <  2:161713/135‑1 (MQ=255)
                                                                                 tCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCAcgcg                                                                  >  4:539687/1‑149 (MQ=255)
                                                                                   aTTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCggg                                                                >  2:756356/1‑149 (MQ=255)
                                                                                     tctcGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGtcattc                                                                                                 >  4:208480/1‑114 (MQ=255)
                                                                                       tcGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTta                                                                                                                                     <  4:362169/76‑1 (MQ=255)
                                                                                       tcGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTta                                                                                                                                     >  3:362169/1‑76 (MQ=255)
                                                                                       tcGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCa                                                            >  3:1089176/1‑149 (MQ=255)
                                                                                       tcGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCa                                                            >  3:196332/1‑149 (MQ=255)
                                                                                       tcGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCa                                                            >  2:1102814/1‑149 (MQ=255)
                                                                                        cGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAg                                                           <  2:969909/149‑1 (MQ=255)
                                                                                                tCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGt                                                   >  1:896722/1‑149 (MQ=255)
                                                                                                     gTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCt                                              >  1:128714/1‑149 (MQ=255)
                                                                                                        aaGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCtttt                                           >  4:417771/1‑149 (MQ=255)
                                                                                                          gTCAATACGACAGCAAGCACGAAAAAGTGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTgg                                         >  1:274810/1‑149 (MQ=255)
                                                                                                          gTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTgg                                         >  2:663166/1‑149 (MQ=255)
                                                                                                           tCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGt                                        <  4:315098/149‑1 (MQ=255)
                                                                                                           tCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGt                                        >  2:556361/1‑149 (MQ=255)
                                                                                                             aaTACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTgtggt                                                                               <  2:994782/108‑1 (MQ=255)
                                                                                                             aaTACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTgtggt                                                                               >  1:994782/1‑108 (MQ=255)
                                                                                                              aTACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCg                                     >  4:175137/1‑149 (MQ=255)
                                                                                                                aCGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCgg                                                                 >  3:1059967/1‑119 (MQ=255)
                                                                                                                aCGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCgg                                                                 <  4:1059967/119‑1 (MQ=255)
                                                                                                                    cAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGa                               <  4:196332/149‑1 (MQ=255)
                                                                                                                     agcaagcaCGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGAc                              >  4:811281/1‑149 (MQ=255)
                                                                                                                        aagcaCGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGt                           <  4:617417/149‑1 (MQ=255)
                                                                                                                        aagcaCGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGt                           <  1:663166/149‑1 (MQ=255)
                                                                                                                          gcaCGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCgg                                                                 >  3:512334/1‑109 (MQ=255)
                                                                                                                              gAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGAt                       >  3:15229/1‑147 (MQ=255)
                                                                                                                              gAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGAt                       <  4:15229/147‑1 (MQ=255)
                                                                                                                              gAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTg                     >  2:195584/1‑149 (MQ=255)
                                                                                                                                   aGGGAGCGAT‑AATTCATCGCTCTCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTcc                >  1:1098132/1‑149 (MQ=255)
                                                                                                                                    gggAGCGATAATTC‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGt                                        <  2:628007/124‑1 (MQ=255)
                                                                                                                                     ggAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTgtgg                                                                                <  4:919831/83‑1 (MQ=255)
                                                                                                                                     ggAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTgtgg                                                                                >  3:919831/1‑83 (MQ=255)
                                                                                                                                     ggAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCag              >  3:501653/1‑149 (MQ=255)
                                                                                                                                     ggAGCGATAATTC‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCag              >  1:217538/1‑149 (MQ=255)
                                                                                                                                     ggAGCGAT‑AATTCATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGt                                                               >  4:781396/1‑100 (MQ=255)
                                                                                                                                     ggAGCGAT‑AATTCATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTATTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGt                                                               <  3:781396/100‑1 (MQ=255)
                                                                                                                                      gAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCaga             >  1:524680/1‑149 (MQ=255)
                                                                                                                                      gAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCaga             <  1:201998/149‑1 (MQ=255)
                                                                                                                                      gAGCGATAATTC‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCaga             >  1:37966/1‑149 (MQ=255)
                                                                                                                                       aGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCagag            <  2:935262/149‑1 (MQ=255)
                                                                                                                                         cGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTg          <  2:896722/149‑1 (MQ=255)
                                                                                                                                          gATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGa         >  3:1023044/1‑149 (MQ=255)
                                                                                                                                                tt‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCgg   <  3:811281/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                    
ATCGCTGGTCTTCTAATCGTTAATTTCTTGTTTCGGGAGCAGATATCCATATGTGCTGGTTTCCGGTAAACAGATGTGCGCTCATTCTCGTACTTATCCCCGTTAAGTCAATACGACAGCAAGCACGAAAAAGGGAGCGATGAATT‑ATCGCTCCCTTGTCTTATAACCATTCAGACATGGTTAGTGTTTCTTGTCATTCATCACAATATAGTGTGGTGAACGTGCACGCGGGTGCAGGAAGAAGTGACCTTTTGGTGCGTTCTGACTGTTGATTGGCTCCAGAGTGACGCCGGT  >  NC_000913/1568774‑1569067

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 27 ≤ ATCG/ATCG < 34 ≤ ATCG/ATCG < 39 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap: ‑    Deleted base: ‑