Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,890,754 C→A 6.3% D172Y (GAT→TAT)  tsaB ← tRNA(ANN) t(6)A37 threonylcarbamoyladenosine modification protein; binding partner and protease for TsaD

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,890,7540CA 6.3% 338.4 / 11.1 99D172Y (GAT→TAT) tsaBtRNA(ANN) t(6)A37 threonylcarbamoyladenosine modification protein; binding partner and protease for TsaD
Reads supporting (aligned to +/- strand):  ref base C (43/49);  new base A (3/3);  total (46/53)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 5.15e-01

GTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCA  >  NC_000913/1890608‑1890898
                                                                                                                                                  |                                                                                                                                                
gttACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCg                                                                                                                                                >  1:68504/1‑149 (MQ=255)
 ttACGTAAATAAACCGGTTGGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAGACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCgg                                                                                                                                               >  4:553537/1‑149 (MQ=255)
  tACGTAAATAAACCGGTTCGGCAAGTTCCACCGCCACCGTTTTACCATCAGCAAACATCTGACACGCAATCGGCAGCATACCTTCAGCAGCAGGCACTAACACTTCGCCATCGCGCAAAACCAGCCCGCTTTCTTTACCGAGATCCGGc                                                                                                                                              <  1:237750/149‑1 (MQ=255)
     gTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCaa                                                                                                                                           >  1:134749/1‑149 (MQ=255)
     gTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCaa                                                                                                                                           >  3:696228/1‑149 (MQ=255)
      taaataaaCCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAg                                                                                                                                          >  4:686133/1‑149 (MQ=255)
           aaaCCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGc                                                                                                                                     >  1:494450/1‑149 (MQ=255)
            aaCCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGcc                                                                                                                                    >  1:186429/1‑149 (MQ=255)
                     ggCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTcc                                                                                                                           <  3:801759/149‑1 (MQ=255)
                     ggCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTAAGCAGCCGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGtgc                                                                                                                           >  2:128723/1‑147 (MQ=255)
                      gCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTccc                                                                                                                          >  4:390551/1‑149 (MQ=255)
                      gCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTccc                                                                                                                          >  4:412029/1‑149 (MQ=255)
                      gCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTccc                                                                                                                          >  3:662312/1‑149 (MQ=255)
                        aTGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTa                                                                                                                        >  3:194941/1‑149 (MQ=255)
                           ttCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCg                                                                                                                     >  2:265106/1‑149 (MQ=255)
                               accgccaccgTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCAc                                                                                                                 <  2:134749/149‑1 (MQ=255)
                                ccgccaccgTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCAcc                                                                                                                >  1:719830/1‑149 (MQ=255)
                                  gccaccgTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCa                                                                                                              >  2:713235/1‑149 (MQ=255)
                                    caccgTTTTACCCTCAGCAAACATCTGACGCGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATACGGCCAAGCt                                                                                                                                        >  1:646132/1‑121 (MQ=255)
                                    caccgTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTa                                                                                                                        <  4:482587/137‑1 (MQ=255)
                                    caccgTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTa                                                                                                                        >  3:482587/1‑137 (MQ=255)
                                    caccgTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCAtt                                                                                                            <  3:638658/149‑1 (MQ=255)
                                    caccgTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATACGGCCAAGCt                                                                                                                                        >  3:782048/1‑121 (MQ=255)
                                    caccgTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATACGGCCAAGCt                                                                                                                                        <  2:646132/121‑1 (MQ=255)
                                    caccgTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATACGGCCAAGCt                                                                                                                                        <  4:636078/121‑1 (MQ=255)
                                    caccgTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATACGGCCAAGCt                                                                                                                                        <  4:782048/121‑1 (MQ=255)
                                    caccgTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGAATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATACGGCCAAGCt                                                                                                                                        >  3:636078/1‑121 (MQ=255)
                                     accgTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTc                                                                                                           <  1:156089/149‑1 (MQ=255)
                                        gTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCAc                                                                                                                 <  4:264504/140‑1 (MQ=255)
                                        gTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCAc                                                                                                                 >  3:264504/1‑140 (MQ=255)
                                         ttttACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCg                                                                                                       <  3:686133/149‑1 (MQ=255)
                                         ttttACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCg                                                                                                       >  3:331653/1‑149 (MQ=255)
                                         ttttACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTGTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCg                                                                                                       >  4:225958/1‑149 (MQ=255)
                                              cccTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAg                                                                                                  >  2:619460/1‑149 (MQ=255)
                                               ccTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCACGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGc                                                                                                 >  4:165141/1‑149 (MQ=255)
                                                 tCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTACGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTg                                                                                               >  2:664151/1‑149 (MQ=255)
                                                   aGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGtt                                                                                             <  1:226803/149‑1 (MQ=255)
                                                   aGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGtt                                                                                             <  4:264015/149‑1 (MQ=255)
                                                   aGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGtt                                                                                             <  3:193691/149‑1 (MQ=255)
                                                   aGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGtt                                                                                             <  1:811853/149‑1 (MQ=255)
                                                      aaaCATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCa                                                                                          >  1:243090/1‑149 (MQ=255)
                                                      aaaCATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCa                                                                                          >  4:441743/1‑149 (MQ=255)
                                                      aaaCATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCa                                                                                          >  2:629070/1‑149 (MQ=255)
                                                      aaaCATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCCACCGTCACCCATTCGCCGGAAAGCTGTTGCa                                                                                          <  4:194941/149‑1 (MQ=255)
                                                       aaCATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCAt                                                                                         <  3:667732/149‑1 (MQ=255)
                                                       aaCATCTGACACGCAATCGGCAGCAAACCTTCAGCAGCAGGCCGTAACACTCCGCCCTCGCGCAAAACCAGCCCGCTCTCTTTTCCGAGATCCGGCCAAGCTTGCCAGCACGTCCCTACCGTCACCCATCCGCCGGAAAGCTGTTGCAt                                                                                         <  2:719830/149‑1 (MQ=255)
                                                        aCATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCAtt                                                                                        >  1:758579/1‑149 (MQ=255)
                                                             tGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTc                                                                                   >  2:687321/1‑149 (MQ=255)
                                                             tGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTc                                                                                   <  3:441743/149‑1 (MQ=255)
                                                              gACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCa                                                                                  <  1:619460/149‑1 (MQ=255)
                                                              gACACGCAATCGGCAGCATATCATCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCa                                                                                  >  1:672237/1‑149 (MQ=255)
                                                               acacGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGc                                                                                                                                     <  3:568195/97‑1 (MQ=255)
                                                               acacGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGc                                                                                                                                     >  4:568195/1‑97 (MQ=255)
                                                               acacGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCAt                                                                                 <  1:51386/149‑1 (MQ=255)
                                                               acacGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCAt                                                                                 >  3:420811/1‑149 (MQ=255)
                                                                cacGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGGGATCGGGCCAAGCTCGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATg                                                                                <  1:128723/149‑1 (MQ=255)
                                                                cacGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTCGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATg                                                                                <  3:18039/149‑1 (MQ=255)
                                                                   gCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCTCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGAc                                                                             <  3:553537/149‑1 (MQ=255)
                                                                         ggCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTc                                                                       <  4:257408/149‑1 (MQ=255)
                                                                          gcagcaTATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCg                                                                      >  3:205137/1‑149 (MQ=255)
                                                                           cagcaTATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCgg                                                                     <  4:331653/149‑1 (MQ=255)
                                                                            agcaTATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCggg                                                                    <  1:265106/149‑1 (MQ=255)
                                                                             gcaTATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGt                                                                   >  3:148067/1‑149 (MQ=255)
                                                                              caTATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGtt                                                                  <  3:165141/149‑1 (MQ=255)
                                                                              caTATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGtt                                                                  >  2:253261/1‑149 (MQ=255)
                                                                               atatCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGttt                                                                 <  4:274761/149‑1 (MQ=255)
                                                                               atatCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGttt                                                                 <  3:412029/149‑1 (MQ=255)
                                                                                tatCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTc                                                                                   <  4:559469/130‑1 (MQ=255)
                                                                                tatCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTc                                                                                   >  3:559469/1‑130 (MQ=255)
                                                                                    ttCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTa                                                            <  4:662312/149‑1 (MQ=255)
                                                                                     tCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTAc                                                           <  2:68504/149‑1 (MQ=255)
                                                                                          agcagGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCtt                                                      <  2:186429/149‑1 (MQ=255)
                                                                                           gcaggcagTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTc                                                     <  4:59262/149‑1 (MQ=255)
                                                                                                 agTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTttc                                               <  4:696228/149‑1 (MQ=255)
                                                                                                  gTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTttct                                              <  3:386584/149‑1 (MQ=255)
                                                                                                    aaCACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGGGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTttcttc                                            <  1:37631/149‑1 (MQ=255)
                                                                                                    aaCACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCgg                                                                     <  1:648468/124‑1 (MQ=255)
                                                                                                    aaCACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTATTGCATTCGTTCATGGACGATTTCgg                                                                     >  2:648468/1‑124 (MQ=255)
                                                                                                     acacTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCa                                           <  1:253261/149‑1 (MQ=255)
                                                                                                                  gcgcAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATAc                              >  2:144724/1‑149 (MQ=255)
                                                                                                                     cAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGa                                                               <  4:617223/113‑1 (MQ=255)
                                                                                                                     cAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGa                                                               >  3:617223/1‑113 (MQ=255)
                                                                                                                      aaaaCCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCTGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGtt                          <  4:577260/149‑1 (MQ=255)
                                                                                                                      aaaaCCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGtt                          <  1:713235/149‑1 (MQ=255)
                                                                                                                          ccAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATGCCGTTTtca                      <  1:629070/149‑1 (MQ=255)
                                                                                                                          ccAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTtca                      <  1:135568/149‑1 (MQ=255)
                                                                                                                           cAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGACCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCTGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTtcat                     <  1:335563/149‑1 (MQ=255)
                                                                                                                            aGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTtcatc                    >  1:787804/1‑149 (MQ=255)
                                                                                                                             gCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTTTTGCATTCGTTCATGGACGATTTCAGGTTTGGGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTtcatca                   >  1:25777/1‑149 (MQ=255)
                                                                                                                             gCCCGCGCGCTTTCCCGAGATTCGGCCAAGCTTGACAGACCGTCACTTCCGTCACCCATTCCGCGGAAAGCTGTTTCTTTCGTTCATGGACGATTTCTGGTTTGAGTACGGCTTCGGTTTCTGCTTCGTGCCAGATACCGATTtcatca                   <  2:494450/149‑1 (MQ=255)
                                                                                                                                cGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGc                <  4:260885/149‑1 (MQ=255)
                                                                                                                                   tctctTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGa                                                               >  2:674374/1‑99 (MQ=255)
                                                                                                                                   tctctTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGa                                                               <  1:674374/99‑1 (MQ=255)
                                                                                                                                   tctctTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGa             >  3:706207/1‑149 (MQ=255)
                                                                                                                                     tctTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGata           >  1:799918/1‑149 (MQ=255)
                                                                                                                                        ttACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTc        <  1:298448/149‑1 (MQ=255)
                                                                                                                                        ttACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTc        <  3:507028/149‑1 (MQ=255)
                                                                                                                                          aCCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCgg      >  4:476357/1‑149 (MQ=255)
                                                                                                                                              agaTCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCa  <  4:19738/149‑1 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GTTACGTAAATAAACCGGTTCGGCATGTTCCACCGCCACCGTTTTACCCTCAGCAAACATCTGACACGCAATCGGCAGCATATCTTCAGCAGCAGGCAGTAACACTTCGCCATCGCGCAAAACCAGCCCGCTCTCTTTACCGAGATCCGGCCAAGCTTGCCAGCCCGTCCCTACCGTCACCCATTCGCCGGAAAGCTGTTGCATTCGTTCATGGACGATTTCGGGTTTGAGTACGGCTTCGGTTTCTTCACCGTGCCAGATACCGTTTTCATCACGCTGATATTCGGCCCA  >  NC_000913/1890608‑1890898

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 16 ≤ ATCG/ATCG < 31 ≤ ATCG/ATCG < 38 ≤ ATCG/ATCG < 40 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: