Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,909,592 T→C 5.4% K170E (AAA→GAA)  kdgR ← KDG regulon transcriptional repressor

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,909,5920TC 5.4% 231.6 / 2.7 74K170E (AAA→GAA) kdgRKDG regulon transcriptional repressor
Reads supporting (aligned to +/- strand):  ref base T (43/27);  new base C (2/2);  total (45/29)
Fisher's exact test for biased strand distribution p-value = 6.42e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 8.30e-01

ACGCCAAAGCGATCAAATACCGGTACCGCAATGCATCGCAGCCCTTCTTCCTGCTCTTCATTATCTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGT  >  NC_000913/1909444‑1909735
                                                                                                                                                    |                                                                                                                                               
aCGCCAAAGCGATCAATTCCCGGTTCCGCAATGCATCGCAGCCCTTCTTCCTGCTCTTCATTATCTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGttt                                                                                                                                                 <  3:10729/149‑1 (MQ=255)
 cGCCAAAGCGATCAAATACCGGTACCGCAATGCATCGCAGCCCTTCTTCCTGCTCTTCATTATCTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTta                                                                                                                                                <  6:234106/149‑1 (MQ=255)
 cGCCAAAGCGATCAAATACCGGTACCGCAAGGCATCGCAGCCCTTCTTCCTGCTCTTCATTATCTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTta                                                                                                                                                <  3:154766/149‑1 (MQ=255)
  gCCAAAGCGATCAAATACCGGTACCGCAATGCATCGCAGCCCTTCTTCCTGCTCTTCATTATCTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTtat                                                                                                                                               <  7:237304/149‑1 (MQ=255)
  gCCAAAGCGATCAAATACCGGTACCGCAATGCATCGCAGCCCTTCTTCCTGCTCTTCATTATCTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTta                                                                                                                                                >  3:267067/1‑148 (MQ=255)
   ccAAAGCGATCAAATACCGGTACCGCAATGCATCGCAGCCCTTCTTCCTGCTCTTCATTATCTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTtata                                                                                                                                              <  4:151836/149‑1 (MQ=255)
    cAAAGCGATCAAATACCGGTACCGCAATGCATCGCAGCCCTTCTTCCTGCTCTTCATTATCTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATAc                                                                                                                                             <  3:158765/149‑1 (MQ=255)
                       tACCGCAATGCATCGCAGCCCTTCTTCCTGCTCTTCATTATCTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAAttt                                                                                                                          <  8:234170/149‑1 (MQ=255)
                         ccGCAATGCATCGCAGCCCTTCTTCCTGCTCTTCATTATCTACGCCATACCCCTGCTCGCGCCCCTGGTCCAGAACGGGTAATAACGCGTCTGTACTGGGGCTGGTCCGCTCGGTACTGCGGTTATACTCCCCGCCCTCAAGAATTTGc                                                                                                                        >  8:174027/1‑149 (MQ=255)
                           gCAATGCATCGCAGCCCTTCTTCCTGCTCTTCATTATCTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGctt                                                                                                                      <  1:414201/149‑1 (MQ=255)
                                gCATCGCAGCCCTTCTTCCTGCTCTTCATTATCTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGCGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGcttcactt                                                                                                                 >  2:415432/1‑149 (MQ=255)
                                  aTCGCAGCCCTTCTTCCTGCTCTTCATTATCTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTc                                                                                                                                  >  6:321896/1‑130 (MQ=255)
                                  aTCGCAGCCCTTCTTCCTGCTCTTCATTATCTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTc                                                                                                                                  <  5:321896/130‑1 (MQ=255)
                                        gCCCTTCTTCCTGCTCTTCATTATCTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCatcgcga                                                                                                         <  4:160274/149‑1 (MQ=255)
                                        gCCCTTCTTCCTGCTCTTCATTATCTTCGCCATACCCCTGCTCGCGCACCTAGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGGCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCatcgcga                                                                                                         <  4:105706/149‑1 (MQ=255)
                                         cccTTCTTCCTGCTCTTCATTATCTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGttt                                                                                                                                                 <  5:241823/108‑1 (MQ=255)
                                         cccTTCTTCCTGCTCTTCATTATCTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGttt                                                                                                                                                 >  6:241823/1‑108 (MQ=255)
                                          ccTTCTTCCTGCTCTTCATTATCTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCatcgcgatc                                                                                                       >  7:382095/1‑149 (MQ=255)
                                                   tGCTCTTCATTATCTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGc                                                                                              >  4:261044/1‑149 (MQ=255)
                                                    gCTCTTCATTATCTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGc                                                                                              >  2:216294/1‑148 (MQ=255)
                                                     ctctTCATTATCTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGcc                                                                                             <  4:267067/148‑1 (MQ=255)
                                                        ttCATTATCTTCGCCATACCCCTGATAGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGAGTTTATAATCCACGACCTCAAGAATTTGCTTCACGTCATCGCCAGCGCGCCATGCcag                                                                                           >  6:255718/1‑147 (MQ=255)
                                                           aTTATCTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTa                                                                                      >  7:256356/1‑149 (MQ=255)
                                                           aTTATCTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTAGGTACTGCGTTTATACTCAACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTa                                                                                      >  7:311704/1‑149 (MQ=255)
                                                                cTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCtt                                                                                  >  8:293457/1‑148 (MQ=255)
                                                                cTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCtt                                                                                  <  7:293457/148‑1 (MQ=255)
                                                                cTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCGTCACTTCATCGCGATCGCGCCACGCCAGCAGTACCTTa                                                                                 >  6:138018/1‑149 (MQ=255)
                                                                         aCCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATcgcg                                                                        >  6:310297/1‑149 (MQ=255)
                                                                         aCCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACGTCATCGCGATCGCGCCATGCCAGCAGTACCTTAcc                                                                               >  7:356930/1‑142 (MQ=255)
                                                                          ccccTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCgg                                                                       <  3:261044/149‑1 (MQ=255)
                                                                           cccTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCAt                                                                                                >  5:332066/1‑123 (MQ=255)
                                                                           cccTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGt                                                                      >  1:61262/1‑149 (MQ=255)
                                                                                cTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGt                                                                 <  8:382095/149‑1 (MQ=255)
                                                                                cTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGt                                                                 <  1:216294/149‑1 (MQ=255)
                                                                                 tCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTa                                                                >  4:213210/1‑149 (MQ=255)
                                                                                     gcACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGc                                                            >  8:155785/1‑149 (MQ=255)
                                                                                         cTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCt                                                                   >  1:432054/1‑138 (MQ=255)
                                                                                         cTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCt                                                                   <  2:432054/138‑1 (MQ=255)
                                                                                            gTCCAGAACGGGGAATAACGCTTCTGTACTGGTGAGGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTac                                                     <  7:174027/149‑1 (MQ=255)
                                                                                             tCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTac                                                     >  3:302682/1‑148 (MQ=255)
                                                                                              ccAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTacga                                                   >  5:420967/1‑149 (MQ=255)
                                                                                               cAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGa                                                         >  5:11972/1‑142 (MQ=255)
                                                                                               cAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGa                                                         <  6:11972/142‑1 (MQ=255)
                                                                                                 gAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGc                                                >  1:250783/1‑149 (MQ=255)
                                                                                                 gAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGc                                                >  6:293135/1‑149 (MQ=255)
                                                                                                      ggTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTCATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAAt                                           >  8:389971/1‑149 (MQ=255)
                                                                                                      ggTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTCATACTCCACGACCTCAAGAATTTGCTTCACTTCATCGAGATCGCGCCATTCCAGCAGTACCTTAACAATAGCGGTGCTGTACAGATGATTACGACGCCCaa                                            >  8:396064/1‑148 (MQ=255)
                                                                                                       gTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATc                                          >  7:374552/1‑149 (MQ=255)
                                                                                                        taataaCGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACAAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATcc                                         >  5:308169/1‑149 (MQ=255)
                                                                                                          ataaCGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCg                                        <  5:310297/148‑1 (MQ=255)
                                                                                                          ataaCGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGt                                       <  4:100479/149‑1 (MQ=255)
                                                                                                          ataaCGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGt                                       <  7:155785/149‑1 (MQ=255)
                                                                                                           taaCGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGt                                       >  8:14013/1‑148 (MQ=255)
                                                                                                           taaCGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGt                                       >  8:131790/1‑148 (MQ=255)
                                                                                                            aaCGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTAc                                                                                >  3:46788/1‑106 (MQ=255)
                                                                                                            aaCGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTAc                                                                                <  4:46788/106‑1 (MQ=255)
                                                                                                            aaCGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTg                                      >  7:363431/1‑148 (MQ=255)
                                                                                                            aaCGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTACACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTg                                      >  4:31614/1‑148 (MQ=255)
                                                                                                                 ttCTGTACTGGTGATGGTCCGCTCGTTACTGCGTTTAAACGCCACGCCCTCAAGAATTTGCTTCACTACAACGCGATCGCGCCAGGCCAGCAGTACATTACCAATCGCGGTGCTGTACAGCGGATTACGCCGCCCAAACCTTCAAAACa                                >  8:5189/1‑149 (MQ=255)
                                                                                                                 ttCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCCCTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACa                                >  2:116884/1‑149 (MQ=255)
                                                                                                                   cTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCCCAAGAATTTGCTTCACTTCATCGCGCTCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACa                                <  6:199372/147‑1 (MQ=255)
                                                                                                                       aCTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATgcgc                           >  3:220585/1‑148 (MQ=255)
                                                                                                                         tGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGAGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCaaa                        >  5:160517/1‑149 (MQ=255)
                                                                                                                           gTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTACACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAAtt                      <  3:31614/149‑1 (MQ=255)
                                                                                                                              aTGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTa                   <  2:61262/149‑1 (MQ=255)
                                                                                                                              aTGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTa                   <  6:160517/149‑1 (MQ=255)
                                                                                                                               tGGTCCGCTCGGTACTGCGTTCATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTTAATACATGCGCAAATTGTa                   <  7:389971/148‑1 (MQ=255)
                                                                                                                               tGGTCCGCTCGGTACTGCGTTCATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGAGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTa                   <  7:396064/148‑1 (MQ=255)
                                                                                                                                 gTCCGCTCTGTACTGCGTTTATACTCCACGCCCTCAAGAATTTTCTTCACTTCATCGCGATAGCGCAATGCCAGCAGTACCTTACCAATCTAGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACAt                <  5:255718/149‑1 (MQ=255)
                                                                                                                                 gTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACa                 >  2:134280/1‑148 (MQ=255)
                                                                                                                                  tCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATa               >  7:408962/1‑149 (MQ=255)
                                                                                                                                      cTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATagag            >  4:169848/1‑148 (MQ=255)
                                                                                                                                      cTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGt           >  4:14440/1‑149 (MQ=255)
                                                                                                                                         ggTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTg                     <  7:240124/136‑1 (MQ=255)
                                                                                                                                         ggTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTg                     >  8:240124/1‑136 (MQ=255)
                                                                                                                                              tGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAAtttt    >  8:23592/1‑148 (MQ=255)
                                                                                                                                              tGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTtg   >  1:113126/1‑149 (MQ=255)
                                                                                                                                              tGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTtg   >  5:20510/1‑149 (MQ=255)
                                                                                                                                                cGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAAtttt    <  7:23592/146‑1 (MQ=255)
                                                                                                                                                cGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTtgt  >  4:200151/1‑148 (MQ=255)
                                                                                                                                                cGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTtgt  <  8:363431/148‑1 (MQ=255)
                                                                                                                                                cGTTTATACTCCACGCCCTCAAGAATTTGCTTAACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACTCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTtgt  <  4:36581/148‑1 (MQ=255)
                                                                                                                                                    |                                                                                                                                               
ACGCCAAAGCGATCAAATACCGGTACCGCAATGCATCGCAGCCCTTCTTCCTGCTCTTCATTATCTTCGCCATACCCCTGCTCGCGCACCTGGTCCAGAACGGGTAATAACGCTTCTGTACTGGTGATGGTCCGCTCGGTACTGCGTTTATACTCCACGCCCTCAAGAATTTGCTTCACTTCATCGCGATCGCGCCATGCCAGCAGTACCTTACCAATCGCGGTGCTGTACAGCGGATTACGACGCCCAATCCGTGAATACATGCGCAAATTGTACATAGAGTCAATTTTGT  >  NC_000913/1909444‑1909735

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: ‑