Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 3,581,737 T→A 5.0% intergenic (+111/‑126) yhhY → / → yhhZ aminoacyl nucleotide detoxifying acetyltransferase/putative Hcp1 family polymorphic toxin protein; putative colicin‑like DNase/tRNase activity

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009133,581,7370TA 5.0% 307.1 / 3.5 97intergenic (+111/‑126)yhhY/yhhZaminoacyl nucleotide detoxifying acetyltransferase/putative Hcp1 family polymorphic toxin protein; putative colicin‑like DNase/tRNase activity
Reads supporting (aligned to +/- strand):  ref base T (48/44);  new base A (2/3);  total (50/47)
Fisher's exact test for biased strand distribution p-value = 6.71e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

AATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAGGACTTTTCAGGGATGAGTAATATTGTTTACCTGA  >  NC_000913/3581589‑3581884
                                                                                                                                                    |                                                                                                                                                   
aaTATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACt                                                                                                                                                     <  2:16129/149‑1 (MQ=255)
  tatGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTc                                                                                                                                                   >  5:42639/1‑149 (MQ=255)
        gATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTaaa                                                                                                                                              <  2:241410/148‑1 (MQ=255)
         aTGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAgt                                                                                                                                            >  4:35730/1‑149 (MQ=255)
           gCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAgtg                                                                                                                                           >  6:193743/1‑148 (MQ=255)
           gCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAgtg                                                                                                                                           >  6:222569/1‑148 (MQ=255)
              tattatATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTaa                                                                                                                                        >  5:153925/1‑148 (MQ=255)
              tattatATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAat                                                                                                                                       <  2:70020/149‑1 (MQ=255)
                   tatGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAAtattt                                                                                                                                   >  3:39029/1‑148 (MQ=255)
                     tGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAAtatttta                                                                                                                                 <  7:60790/148‑1 (MQ=255)
                          cgGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTAtttttt                                                                                                                           <  2:45775/149‑1 (MQ=255)
                           gggTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTa                                                                                                                          <  5:117370/149‑1 (MQ=255)
                           gggTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTATTTaa                                                                                                                         <  6:96263/149‑1 (MQ=255)
                             gTGAAGTAAGATCGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAat                                                                                                                        <  2:95756/149‑1 (MQ=255)
                                 aGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAatata                                                                                                                     >  7:137410/1‑148 (MQ=255)
                                  ttAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACATCCTAAAGTGTAATATTTTATTTTTTAAtatata                                                                                                                   >  7:175996/2‑149 (MQ=255)
                                  gTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACATCCTAAAGTGTAATATTTTATTTTTTAAtatata                                                                                                                   >  7:175999/1‑149 (MQ=255)
                                   tAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAgt                                                                                                                                            <  7:190241/123‑1 (MQ=255)
                                     aGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGc                                                                                                                >  2:23569/1‑149 (MQ=255)
                                      gATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGcc                                                                                                               >  1:178734/1‑149 (MQ=255)
                                      gATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGcc                                                                                                               >  1:182152/1‑149 (MQ=255)
                                       aTAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGcc                                                                                                               >  2:33819/1‑148 (MQ=255)
                                        tAGTGCCCTTTTTCTTAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTa                                                                                                             <  5:122029/149‑1 (MQ=255)
                                                tttttCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTcc                                                                                                     >  1:119742/1‑149 (MQ=255)
                                                  tttCTGAGATGGAAAAAGGGTGTCATTCAAACTCGACATACCTTCCTTTAAGGTAGTTATTTGCCCAATACATATCTTGAGATTTAACTCATACTTACTTCCTAAAGTGTAAAATTTTATTTTTTAATATATACGCCTCCAATTTcc                                                                                                     <  4:222121/147‑1 (MQ=255)
                                                   ttCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTg                                                                                                   <  2:401/148‑1 (MQ=255)
                                                   ttCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTg                                                                                                   <  6:204708/148‑1 (MQ=255)
                                                   ttCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTg                                                                                                   <  4:39029/148‑1 (MQ=255)
                                                    tCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACTCCTACAATTTCCTg                                                                                                   >  6:235773/1‑147 (MQ=255)
                                                     cTGATATTGAAAAAGGGTGTCATTCAAAATCGACATACCTTACTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGa                                                                                                 <  5:222569/148‑1 (MQ=255)
                                                     cTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGa                                                                                                 <  7:100890/148‑1 (MQ=255)
                                                      tGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAg                                                                                                >  8:72465/1‑148 (MQ=255)
                                                        agaTGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTaaa                                                                                                                                              <  7:133464/100‑1 (MQ=255)
                                                        agaTGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGta                                                                                              >  8:55710/1‑148 (MQ=255)
                                                        agaTGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTAGAGta                                                                                              >  1:185440/1‑148 (MQ=255)
                                                        agaTGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTAATTTGGCCAATACATATATTGAGATTTAACTCAGAATTACTTCGTaaa                                                                                                                                              >  8:133464/1‑100 (MQ=255)
                                                          aTGGAAAAAGGGTGTCATACAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATAGAGATATAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTACTGGAGtaaa                                                                                            >  2:228302/1‑148 (MQ=255)
                                                             gAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGtaaataaa                                                                                        <  4:229069/149‑1 (MQ=255)
                                                             gAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTACCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGtaaataaa                                                                                        <  3:30476/149‑1 (MQ=255)
                                                                acaGGGTGTCATTCAAAATCGCCATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCAAAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAataaataa                                                                                     <  2:16062/147‑1 (MQ=255)
                                                                 aaGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAAc                                                                                    <  5:193743/149‑1 (MQ=255)
                                                                   gggTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACaa                                                                                  >  5:76634/1‑149 (MQ=255)
                                                                   gggTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACaa                                                                                  >  7:94565/1‑149 (MQ=255)
                                                                   gggTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACaa                                                                                  <  5:150020/149‑1 (MQ=255)
                                                                    ggTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAAt                                                                                 <  1:91350/149‑1 (MQ=255)
                                                                     gtgtCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAAt                                                                                 >  6:182515/1‑148 (MQ=255)
                                                                      tgtCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACATCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTa                                                                               <  8:175999/149‑1 (MQ=255)
                                                                      tgtCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACATCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTa                                                                               <  8:175996/149‑1 (MQ=255)
                                                                         cATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAgt                                                                                                                                            >  8:41729/1‑85 (MQ=255)
                                                                         cATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAgt                                                                                                                                            <  7:41729/85‑1 (MQ=255)
                                                                         cATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAAc                                                                             >  2:24606/1‑148 (MQ=255)
                                                                          aTTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACaa                                                                           <  8:32187/149‑1 (MQ=255)
                                                                          aTTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACaa                                                                           <  7:67447/149‑1 (MQ=255)
                                                                             cAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGc                                                                         <  5:182515/148‑1 (MQ=255)
                                                                                  tCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCataata                                                                   <  3:35730/149‑1 (MQ=255)
                                                                                     aCATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTg                                                                >  1:38055/1‑149 (MQ=255)
                                                                                       aTACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGcc                                                              <  4:104326/149‑1 (MQ=255)
                                                                                        tACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGACTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGcaa                                                             >  5:10281/1‑147 (MQ=255)
                                                                                         aCCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTAATTCCTAAAGTGTAATATTTTATTTTTTAATATATAAGCCTAAAATTTAATGGATTAAATAAATAAAAATTAAAAATAATAATAtttaaat                                                            >  1:248427/1‑144 (MQ=255)
                                                                                          ccttccttTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCAt                                                            <  1:23569/148‑1 (MQ=255)
                                                                                          ccttccttTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCAt                                                            <  4:126839/148‑1 (MQ=255)
                                                                                           cttccttTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCAtt                                                           >  4:40985/1‑148 (MQ=255)
                                                                                              ctatctaGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGAt                                                        <  3:2731/142‑1 (MQ=255)
                                                                                               cttTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATa                                                       >  6:179661/1‑148 (MQ=255)
                                                                                                 ggaaGCTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATaaaa                                                    <  8:157359/147‑1 (MQ=255)
                                                                                                  tAAGGTATTTATTTTCCCAATACATATATTGAGATTTAACTCATAATTACATCCTAAAGAGTAATAATATATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATaaaa                                                    <  3:56572/148‑1 (MQ=255)
                                                                                                     ggTATTTATTTTCCCAATACATATATTGAGATTTAACTCATAATTAGTTCCTAAAGTGTAATATTTTATTTTTTAATATATACCCCTACAATTTCCTGGAGTATATAAATAACAATTATCAAGCATAATATTGCCATTGATAAAATAGc                                                <  1:228302/149‑1 (MQ=255)
                                                                                                     ggTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGAt                                                        >  4:2731/1‑141 (MQ=255)
                                                                                                     ggTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGc                                                <  2:182152/149‑1 (MQ=255)
                                                                                                     ggTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGACTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGAATAATATTGCCATTGATAAAATAg                                                 >  4:142920/1‑148 (MQ=255)
                                                                                                     ggTATTTATTTGACCAAAACATATATTGATATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGc                                                <  2:178734/149‑1 (MQ=255)
                                                                                                      gTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGcc                                                              <  8:537/134‑1 (MQ=255)
                                                                                                      gTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGcc                                                              >  7:537/1‑134 (MQ=255)
                                                                                                      gTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCa                                               <  4:38666/149‑1 (MQ=255)
                                                                                                       tatttattTGCCAAATACATATATTCAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAAAAAGCATAATATTGCCATTGATAAAATAGCAt                                              >  8:157247/1‑149 (MQ=255)
                                                                                                        atttatttGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCGATGGATAAAATAGCATg                                             >  6:175155/1‑149 (MQ=255)
                                                                                                           tatttGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGcc                                           <  2:119742/148‑1 (MQ=255)
                                                                                                            atttGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCATACAATTTCCTGGAGTAAATAAATAAAAATTAACAAGAATAATATTGCCATTGATAAAATAGCATGCCa                                          >  3:168788/1‑148 (MQ=255)
                                                                                                             tttGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCAt                                         <  5:235773/148‑1 (MQ=255)
                                                                                                                   cAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAGGa                                  >  4:234317/1‑149 (MQ=255)
                                                                                                                    aaTACATATATTGAGATTTACCTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAGGAc                                 >  3:132256/1‑149 (MQ=255)
                                                                                                                      tACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAGGACtt                               <  4:132256/149‑1 (MQ=255)
                                                                                                                        cATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAGGACtttt                             <  6:76634/149‑1 (MQ=255)
                                                                                                                         atatatTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAGGACtttt                             >  5:195113/1‑148 (MQ=255)
                                                                                                                         atatatTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAGGACtttt                             >  5:109429/1‑148 (MQ=255)
                                                                                                                              ttGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAGGACTTTTCAGGGa                       >  8:242973/1‑149 (MQ=255)
                                                                                                                               tGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAGGACTTTTCAGGGAt                      >  8:157350/1‑149 (MQ=255)
                                                                                                                                gagaTTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAGGACTTTTCAGGGAt                      >  3:206404/1‑148 (MQ=255)
                                                                                                                                 agaTTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAGGACTTTTCAGGGAt                      >  5:198426/1‑147 (MQ=255)
                                                                                                                                   aTTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAGGACTTTTCAGGGATGAg                   >  5:47026/1‑148 (MQ=255)
                                                                                                                                          tCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAGGACTTTTCAGGGATGAGTAATATTg           >  2:153982/1‑149 (MQ=255)
                                                                                                                                            aTAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAGGACTTTTCAGGGATGAGTAATATTGtt         >  1:239854/1‑149 (MQ=255)
                                                                                                                                            aTAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAGGACTTTTCAGGGATGAGTAATATTGtt         >  8:234178/1‑149 (MQ=255)
                                                                                                                                               aTTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTTGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAGGACTTTTCAGGGATGAGTAATATTGTTTa       <  7:157247/148‑1 (MQ=255)
                                                                                                                                               aTTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAGGACTTTTCAGGGATGAGTAATATTGTTTAc      >  4:191579/1‑149 (MQ=255)
                                                                                                                                                taaCTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAGGACTTTTCAGGGATGAGTAATATTGTTTAcc     >  3:215931/3‑149 (MQ=255)
                                                                                                                                                  aCTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAGGACTTTTCAGGGATGAGTAATATTGTTTACCt    <  1:33819/148‑1 (MQ=255)
                                                                                                                                                   cTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAGGACTTTTCAGGGATGAGTAATATTGTTTACCTGa  <  5:44577/149‑1 (MQ=255)
                                                                                                                                                    ttCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAGGACTTTTCAGGGATGAGTAATATTGTTTACCTGa  >  6:56517/1‑148 (MQ=255)
                                                                                                                                                    |                                                                                                                                                   
AATATGTCGATGCATATTATATGGCGCGGGTGAAGTAAGATAGTGCCCTTTTTCTGAGATGGAAAAAGGGTGTCATTCAAAATCGACATACCTTCCTTTAAGGTATTTATTTGCCCAATACATATATTGAGATTTAACTCATAATTACTTCCTAAAGTGTAATATTTTATTTTTTAATATATACGCCTACAATTTCCTGGAGTAAATAAATAACAATTAACAAGCATAATATTGCCATTGATAAAATAGCATGCCATAAAAGGACTTTTCAGGGATGAGTAATATTGTTTACCTGA  >  NC_000913/3581589‑3581884

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: