Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,628,195 A→G 49.5% Y294C (TAC→TGC)  nadR → trifunctional protein: nicotinamide mononucleotide adenylyltransferase, ribosylnicotinamide kinase, transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,628,1950AG49.5% ‑4.6 / 148.9 99Y294C (TAC→TGC) nadRtrifunctional protein: nicotinamide mononucleotide adenylyltransferase, ribosylnicotinamide kinase, transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base A (25/25);  new base G (30/19);  total (55/44)
Fisher's exact test for biased strand distribution p-value = 3.14e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.83e-01

Only 100 of 101 total aligned reads displayed.

AAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATC  >  NC_000913/4628058‑4628343
                                                                                                                                         |                                                                                                                                                    
aaacaaacTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTg                                                                                                                                           >  3:82610/1‑149 (MQ=255)
 aacatacTTGCGAATATCTTCTTCACGACCAGTGCGTTGGAATATGGCCTCGATTATGTCTTTTCACACCTCGGCGTTTATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATa                                                                                                                                                      <  2:375049/137‑1 (MQ=255)
 aacaaacTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGc                                                                                                                                          >  5:233927/1‑149 (MQ=255)
  acaaacTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCg                                                                                                                                         <  2:57692/149‑1 (MQ=255)
           ccAATATCTTCATCACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGTTGATGTGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCGCAATACATTGATTTTGCGGTGAAtta                                                                                                                                 <  6:22201/148‑1 (MQ=255)
           ccAATATCTTCAACACCACCAGTGCGTGGGAATGTGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAata                                                                                                                                 <  4:365916/148‑1 (MQ=255)
           ccAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAata                                                                                                                                 <  8:254676/148‑1 (MQ=255)
            cAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATg                                                                                                                               >  2:76998/1‑149 (MQ=255)
            cAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAGTATg                                                                                                                               >  1:97818/1‑149 (MQ=255)
               tatCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGcc                                                                                                                             >  6:30724/1‑148 (MQ=255)
                   ttCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATa                                                                                                                         >  1:119503/1‑148 (MQ=255)
                     cAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATa                                                                                                                         <  1:147962/146‑1 (MQ=255)
                      aaCACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGt                                                                                                                     >  1:190917/1‑149 (MQ=255)
                      aaCACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAg                                                                                                                      <  4:53096/148‑1 (MQ=255)
                      aaCACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGt                                                                                                                     >  2:276847/1‑149 (MQ=255)
                           caccaGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCAt                                                                                                                >  1:353908/1‑149 (MQ=255)
                           caccaGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCAt                                                                                                                >  6:114808/1‑149 (MQ=255)
                           caccaGTGCGTGGGAATATGGCCCCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCAt                                                                                                                >  4:243531/1‑149 (MQ=255)
                              caGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCTGCGGTGATGAGATCGCATTGCAGTATTCTTACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTa                                                                                                             <  6:335991/149‑1 (MQ=255)
                               aGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTAt                                                                                                            <  1:337600/149‑1 (MQ=255)
                                 tGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATc                                                                                                           >  8:384136/1‑148 (MQ=255)
                                     tGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGTGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATAc                                                                                                      >  4:314753/1‑149 (MQ=255)
                                     tGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAGTAAAGTGGCATTTATCGATAc                                                                                                      >  1:330356/1‑149 (MQ=255)
                                     tGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATa                                                                                                       >  6:336704/1‑148 (MQ=255)
                                     tGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATAc                                                                                                      >  1:71228/1‑149 (MQ=255)
                                        gAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGa                                                                                                   >  8:124303/1‑149 (MQ=255)
                                          atatGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGAtt                                                                                                 <  8:52787/149‑1 (MQ=255)
                                                ccGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTca                                                                                           >  5:241118/1‑149 (MQ=255)
                                                 cgcgATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGt                                                                                                                     <  5:104647/122‑1 (MQ=255)
                                                 cgcgATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGt                                                                                                                     >  6:104647/1‑122 (MQ=255)
                                                 cgcgATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTcac                                                                                          >  4:172177/1‑149 (MQ=255)
                                                   cgATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTcacca                                                                                        <  7:64141/149‑1 (MQ=255)
                                                       tATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTGTGACTACGATAAATTCGCGCTTGGCCACTCTCATTGCATTGATTTTGCGGTGAAATATGCGAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTc                                                                                     <  1:16856/148‑1 (MQ=255)
                                                        aTGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAg                                                                                   >  3:357551/1‑149 (MQ=255)
                                                         tGTCTTTTCACACCTCGGCGGTGATGAGCTCGCATTGCAGTATGCTGACTACGATAAAATCGCTCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAg                                                                                   <  1:135267/148‑1 (MQ=255)
                                                         tGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAATTATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAg                                                                                   <  2:97818/148‑1 (MQ=255)
                                                             ttttCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTcac                                                                                          <  6:292934/137‑1 (MQ=255)
                                                             ttttCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTcac                                                                                          >  5:292934/1‑137 (MQ=255)
                                                             ttttCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGtt                                                                              >  4:287609/1‑149 (MQ=255)
                                                             ttttCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGtt                                                                              >  5:122807/1‑149 (MQ=255)
                                                             ttttCACACCTCGGAGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGtt                                                                              >  8:17505/1‑149 (MQ=255)
                                                               ttCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTc                                                                             <  5:114808/148‑1 (MQ=255)
                                                               ttCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTc                                                                             <  3:243531/148‑1 (MQ=255)
                                                                  acacCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGc                                                                          >  2:60014/1‑148 (MQ=255)
                                                                  acacCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGc                                                                          >  5:231111/1‑148 (MQ=255)
                                                                  acacCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCa                                                                         <  7:124303/149‑1 (MQ=255)
                                                                    acCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCaa                                                                        >  4:142495/1‑148 (MQ=255)
                                                                       tCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAAGACATTGATTTTGCGGTGAAATATGCAAATAAAGTGGCAATTATCGAGACCGATTTTGTCACCACTCAGTCGTTCTGCAAAAAg                                                                    >  7:351522/1‑149 (MQ=255)
                                                                         ggcggTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCTGTGAAATATGCCAATAAAGTTGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGt                                                                   <  1:71240/148‑1 (MQ=255)
                                                                         ggcggTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTa                                                                  <  7:384136/149‑1 (MQ=255)
                                                                          gcggTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTAc                                                                 <  5:30724/149‑1 (MQ=255)
                                                                            ggTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGa                                                               <  6:233927/149‑1 (MQ=255)
                                                                             gTGATGAGATCGAATTGGAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGaa                                                              <  3:288301/149‑1 (MQ=255)
                                                                              tgatgaGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAg                                                             >  1:391120/1‑149 (MQ=255)
                                                                                atgaGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAggg                                                           <  2:190917/149‑1 (MQ=255)
                                                                                  gagaTCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCg                                                         >  7:366798/1‑149 (MQ=255)
                                                                                  gagaTCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCaaaaa                                                                     <  1:348627/137‑1 (MQ=255)
                                                                                  gagaTCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTACTGCaaaaa                                                                     >  2:348627/1‑137 (MQ=255)
                                                                                     aTCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCggg                                                       >  2:199294/1‑148 (MQ=255)
                                                                                     aTCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGa                                                      <  3:142495/149‑1 (MQ=255)
                                                                                         cATTGCTGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACAt                                                  <  7:194894/149‑1 (MQ=255)
                                                                                          aTTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATc                                                 >  8:6186/1‑149 (MQ=255)
                                                                                              cAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGtt                                             >  1:42234/1‑149 (MQ=255)
                                                                                               agtagtCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGACAATAAAGTGGGATTTATCGAGAACGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGTGCGGGAACATCCGTTc                                            >  8:302429/1‑149 (MQ=255)
                                                                                               aGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTc                                            <  2:353908/149‑1 (MQ=255)
                                                                                                 tATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGt                                          >  6:267211/1‑149 (MQ=255)
                                                                                                   ttCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAACGTACGAAGGGCGGGAACATCCGTTCGGGc                                        <  7:17505/149‑1 (MQ=255)
                                                                                                   ttCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCcgtgcgtg                                         >  6:169391/1‑148 (MQ=255)
                                                                                                     cTGACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCa                                       <  5:336704/148‑1 (MQ=255)
                                                                                                     cTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCa                                       <  6:122807/148‑1 (MQ=255)
                                                                                                       gACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTg                                         <  2:261132/144‑1 (MQ=255)
                                                                                                       gACTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTg                                         >  1:261132/1‑144 (MQ=255)
                                                                                                        aCTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGgcg                                   >  8:327637/1‑149 (MQ=255)
                                                                                                         cTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGgcgc                                  <  4:82610/149‑1 (MQ=255)
                                                                                                         cTACGATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGgcgc                                  <  1:76998/149‑1 (MQ=255)
                                                                                                             gATAAAATCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCtgat                              >  1:203473/1‑149 (MQ=255)
                                                                                                             gATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCtgat                              >  4:31515/1‑149 (MQ=255)
                                                                                                                aaaaTCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCtgattga                           >  4:122213/1‑149 (MQ=255)
                                                                                                                 aaatgccgcTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCtgattga                           <  6:241118/142‑1 (MQ=255)
                                                                                                                 aaaTCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAGTAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCtgattgat                          <  2:330356/149‑1 (MQ=255)
                                                                                                                 aaaTCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCtgattga                           <  1:60014/148‑1 (MQ=255)
                                                                                                                 aaaTCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCtgattgat                          <  3:31515/149‑1 (MQ=255)
                                                                                                                  aaTCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCtgattgat                          >  7:192111/1‑148 (MQ=255)
                                                                                                                  aaTCGCGCTGGGCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATtgatg                         <  5:139979/149‑1 (MQ=255)
                                                                                                                  aaTCGCGCTGGGCAACGCTCAATGAATTGATTATGCGGTGATATATGACAATAAAGTGGCATTTATCGATACCGATTTTGTCACAACTCAGTCGTACTGCAACAAGTACGAAGGGCTGGAACATCCGTTCGTGCTGTCTCttattgat                          >  2:325830/1‑148 (MQ=255)
                                                                                                                   aTCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATtgatg                         >  3:211798/1‑148 (MQ=255)
                                                                                                                    tCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGaa                       <  3:287609/149‑1 (MQ=255)
                                                                                                                         cTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCg                  >  6:153978/1‑149 (MQ=255)
                                                                                                                          tGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCTTTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGAAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGt                 >  8:326727/1‑149 (MQ=255)
                                                                                                                          tGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGt                 >  1:355093/1‑149 (MQ=255)
                                                                                                                             gCCACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTc              >  7:135169/1‑149 (MQ=255)
                                                                                                                               cACGCTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGa            <  2:391120/149‑1 (MQ=255)
                                                                                                                                  gCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCt         >  8:300189/1‑149 (MQ=255)
                                                                                                                                  gCTCAATACATTGATTTTGAGGTGAAATATCCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGAGGGAACATCCGTTCGTGCAGGCGCAGATTGATGAAAACCGATTCGATCt         >  2:334068/1‑149 (MQ=255)
                                                                                                                                   cTCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTg        <  2:203473/149‑1 (MQ=255)
                                                                                                                                   cTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTg        <  4:357551/149‑1 (MQ=255)
                                                                                                                                    tCAATGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTgg       >  1:198231/1‑149 (MQ=255)
                                                                                                                                      aaTGCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGt      >  7:299175/1‑148 (MQ=255)
                                                                                                                                      aaTACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGt      >  8:379652/1‑148 (MQ=255)
                                                                                                                                         aCATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGATCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATc  >  2:372823/1‑149 (MQ=255)
                                                                                                                                         |                                                                                                                                                    
AAACAAACTTGCCAATATCTTCAACACCACCAGTGCGTGGGAATATGGCCGCGATTATGTCTTTTCACACCTCGGCGGTGATGAGATCGCATTGCAGTATTCTGACTACGATAAAATCGCGCTGGGCCACGCTCAATACATTGATTTTGCGGTGAAATATGCCAATAAAGTGGCATTTATCGATACCGATTTTGTCACCACTCAGGCGTTCTGCAAAAAGTACGAAGGGCGGGAACATCCGTTCGTGCAGGCGCTGATTGATGAATACCGTTTCGATCTGGTGATC  >  NC_000913/4628058‑4628343

Alignment Legend
Aligned base mismatch/match (shaded by quality score): ATCG/ATCG < 3 ≤ ATCG/ATCG < 14 ≤ ATCG/ATCG < 21 ≤ ATCG/ATCG < 32 ≤ ATCG/ATCG < 36 ≤ ATCG/ATCG
Unaligned base: atcg    Masked matching base: atcg    Alignment gap:     Deleted base: