Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,909,732 T→A 6.4% K123I (AAA→ATA)  kdgR ← KDG regulon transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,909,7320TA 6.4% 219.8 / 3.7 72K123I (AAA→ATA) kdgRKDG regulon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (28/39);  new base A (3/2);  total (31/41)
Fisher's exact test for biased strand distribution p-value = 6.46e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 1.00e+00

CTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATCTGGATATCTGCGCTACGAATTAAATCGACGTTTTGTAACGCGCGAGCGCCCAGTTCAAACAATTTCAG  >  NC_000913/1909585‑1909880
                                                                                                                                                   |                                                                                                                                                    
cTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAAtttt                                                                                                                                                     <  1:54174/149‑1 (MQ=255)
    gTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCCTGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTtgtg                                                                                                                                                  >  6:381148/1‑148 (MQ=255)
         taCTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCTCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAATTTGTACATAGAGTCAATTTTGTGAATGt                                                                                                                                             <  1:180326/148‑1 (MQ=255)
            tCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTaaa                                                                                                                                          >  7:381319/1‑148 (MQ=255)
              cACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGt                                                                                                                                             <  3:129982/143‑1 (MQ=255)
              cACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGt                                                                                                                                             >  4:129982/1‑143 (MQ=255)
               aCGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAAt                                                                                                                                      <  7:373261/149‑1 (MQ=255)
                 gCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATAc                                                                                                                                    <  2:33181/149‑1 (MQ=255)
                   ccTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTg                                                                                                                                  <  5:343757/149‑1 (MQ=255)
                     tCAAGAATTTGCTTCACTTCATCGCTATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTTCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTggc                                                                                                                                <  1:301091/149‑3 (MQ=255)
                     tCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTc                                                                                                                                >  5:354723/1‑149 (MQ=255)
                     tCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTc                                                                                                                                >  3:315489/1‑149 (MQ=255)
                     tCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTATGTGAATGTAAACAATACTGTc                                                                                                                                >  5:227802/1‑149 (MQ=255)
                     tCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTATGTGAATGTAAACAATACTGTc                                                                                                                                >  7:387957/1‑149 (MQ=255)
                         gAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCg                                                                                                                            >  1:306062/1‑149 (MQ=255)
                           aTTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTc                                                                                                                          <  7:67006/149‑1 (MQ=255)
                            tttGCTTCACTTCATCGCGATCGAGCCATGCCAGCAGTACCTTTCGAATCGCTGTGCTGTAGAGCGGATTACGACGCCCTATCCGTGAATATATGCTCAAATTGTACATAGAGTCAATTTTGTGATTGTAAACAATAGTGTCTTCGTcc                                                                                                                         <  8:329496/149‑1 (MQ=255)
                               gCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGAGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGt                                                                                                                      >  4:116892/1‑149 (MQ=255)
                                cttcacttcaTCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTg                                                                                                                     >  3:263286/1‑149 (MQ=255)
                                 ttcacttcaTCTCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTgc                                                                                                                    >  1:200531/1‑149 (MQ=255)
                                 ttcacttcaTCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTgc                                                                                                                    >  7:200180/1‑149 (MQ=255)
                                   cacttcaTCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTgcgc                                                                                                                  >  2:269304/1‑149 (MQ=255)
                                   cacttcaTCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTgcgc                                                                                                                  >  2:385337/1‑149 (MQ=255)
                                        caTCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAg                                                                                                              <  8:42312/148‑1 (MQ=255)
                                           cgcgatcgcgCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTgg                                                                                                          >  4:353216/1‑149 (MQ=255)
                                            gcgttcgcgCCATTCAAGCAGTACCTTACCACTCGAGGTGCTGTACATCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGa                                                                                                         <  1:237519/149‑1 (MQ=255)
                                            gcgatcgcgCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGa                                                                                                         <  4:315489/149‑1 (MQ=255)
                                              gatcgcgCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATa                                                                                                       <  8:381319/149‑1 (MQ=255)
                                                tcgcgCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGt                                                                                                     >  1:32813/1‑149 (MQ=255)
                                                      cATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCttt                                                                                               <  6:160213/149‑1 (MQ=255)
                                                       aTGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTg                                                                                              <  7:365988/149‑1 (MQ=255)
                                                       aTGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAAGTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTg                                                                                              >  6:77126/1‑149 (MQ=255)
                                                                gTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCgg                                                                                     <  1:29389/149‑1 (MQ=255)
                                                                gTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCgg                                                                                     <  1:385337/149‑1 (MQ=255)
                                                                   ccTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCttt                                                                                               >  7:345727/1‑136 (MQ=255)
                                                                   ccTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCttt                                                                                               <  8:345727/136‑1 (MQ=255)
                                                                   ccTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGgag                                                                                  <  6:354723/149‑1 (MQ=255)
                                                                   ccTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGga                                                                                   <  4:263286/148‑1 (MQ=255)
                                                                        ccAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCt                                                                              <  6:110900/148‑1 (MQ=255)
                                                                          aaTCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCTCAc                                                                           >  6:13561/1‑149 (MQ=255)
                                                                              gcgGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGGCAGGCGGGAGAGCTCACGCAt                                                                       <  4:81505/149‑1 (MQ=255)
                                                                               cgGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATc                                                                      >  3:293505/1‑149 (MQ=255)
                                                                                ggTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATCt                                                                     <  3:297006/149‑1 (MQ=255)
                                                                                ggTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATCt                                                                     <  1:211021/149‑1 (MQ=255)
                                                                                    cggTACAGCGGAGTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACCAAGAGGCAATTTTGTGCATGTAAACAATCCTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATCTGGa                                                                  <  2:380558/146‑1 (MQ=255)
                                                                                    cTGTACAGCGGATTACGACGCCCAATCCGTGAATACCTGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATCTGGat                                                                 <  8:287059/149‑1 (MQ=255)
                                                                                    cTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATCTGGat                                                                 <  6:261923/149‑1 (MQ=255)
                                                                                    cTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATCTGGa                                                                  >  1:354866/1‑148 (MQ=255)
                                                                                    cTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATCTGGa                                                                  <  7:322615/148‑1 (MQ=255)
                                                                                             ggATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATCTGGATATCTGCGCt                                                        <  2:200531/149‑1 (MQ=255)
                                                                                                  acgacgCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTATGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATCTGGATATCTGCGCTACGa                                                    <  8:387957/148‑1 (MQ=255)
                                                                                                  acgacgCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTATGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATCTGGATATCTGCGCTACGa                                                    <  6:227802/148‑1 (MQ=255)
                                                                                                    gacgCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATCTGGATATCTGCGCTACGAAt                                                  >  1:137810/1‑148 (MQ=255)
                                                                                                        cccAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATCTGGATATCTGCGCTACGAATTAAAt                                             >  4:205840/1‑149 (MQ=255)
                                                                                                         ccAATCCGTGAATCCATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATCTGGATATCTGCGCTACGAATTAAATc                                            <  4:90582/149‑1 (MQ=255)
                                                                                                         ccAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATCTGGATATCTGCGCTACGAATTAAATc                                            <  2:103332/149‑1 (MQ=255)
                                                                                                         ccAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATCTGGATATCTGCGCTACGAATTAAATc                                            >  1:363983/1‑149 (MQ=255)
                                                                                                           aaTCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATCTGGATATCTGCGCTACGAATTAAATCGa                                          >  2:127240/1‑149 (MQ=255)
                                                                                                               cGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATCTGGATATCTGCGCTACGAATTAAATCGACGtt                                      >  3:359916/1‑149 (MQ=255)
                                                                                                                gTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATCTGGATATCTGCGCt                                                        <  7:331638/130‑1 (MQ=255)
                                                                                                                gTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATCTGGATATCTGCGCt                                                        >  8:331638/1‑130 (MQ=255)
                                                                                                                 tGAATACATGCGCAAAATGTACATAGAGTCAATTATGGGAATGTAAACAATACTGTCGTCGTCCAGTCCGCCGAGGTGGATAGTTTCGTTGGTCAGGCGGGAGAGCTCACGCATCGGGATATCTGGGCTACGATTTATATCGACGgttt                                    >  6:109918/1‑149 (MQ=255)
                                                                                                                  gAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATCTGGATATCTGCGCTACGAATTAAATCGACGTTTTg                                   >  1:203880/1‑149 (MQ=255)
                                                                                                                    aTACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTg                                                                                              >  2:5111/1‑88 (MQ=255)
                                                                                                                    aTACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTg                                                                                              <  1:5111/88‑1 (MQ=255)
                                                                                                                    aTACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCCAGGTGTATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATCTGGATATCTGCGCTACGAATTAAATCGACGTTTTGt                                  <  3:116892/148‑1 (MQ=255)
                                                                                                                        aTGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATCTGGATATCTGCGCTACGAATTAAATCGACGTTTTGTAAcgc                             <  1:127240/149‑1 (MQ=255)
                                                                                                                           cgcAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTTGATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATCTGGATATCTGCGCTACGAATTAAATCGACGTTTTGTAAcgcgc                           <  2:306062/148‑1 (MQ=255)
                                                                                                                            gcAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATCTGGATATCTGCGCTACGAATTAAATCGACGTTTTGTAACGCGCGa                         <  1:269304/149‑1 (MQ=255)
                                                                                                                            gcAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATCTGGATATCTGCGCTACGAATTAAATCGACGTTTTGTAACGCGCGa                         <  4:359916/149‑1 (MQ=255)
                                                                                                                                    tACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATCTGGATATCTGCGCTACGAATTAAATCGACGTTTTGTAACGCGCGAGCGCCCa                  <  2:32813/148‑1 (MQ=255)
                                                                                                                                          gagTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATCTGGATATCTGCGCTACGAATTAAATCGACGTTTTGTAACGCGCGAGCGCCCAGTTCaaa           <  3:353216/149‑1 (MQ=255)
                                                                                                                                               aaTTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATCTGGATATCTGCGCTACGAATTAAATCGACGTTTTGTAACGCGCGAGCGCCCAGTTCAAACAAtt      >  8:288054/1‑149 (MQ=255)
                                                                                                                                                   ttGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATCTGGATATCTGCGCTACGAATTAAATCGACGTTTTGTAACGCGCGAGCGCCCAGTTCAAACAATTtcag  >  1:335349/1‑149 (MQ=255)
                                                                                                                                                   |                                                                                                                                                    
CTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGTGAATGTAAACAATACTGTCTTCGTCCAGTGCGCCGAGGTGGATAGTTTCTTTGGTCAGGCGGGAGAGCTCACGCATCTGGATATCTGCGCTACGAATTAAATCGACGTTTTGTAACGCGCGAGCGCCCAGTTCAAACAATTTCAG  >  NC_000913/1909585‑1909880

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