Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 1,788,507 T→G 6.1% Y25D (TAT→GAT)  aroH → 3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tryptophan repressible

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009131,788,5070TG 6.1% 347.4 / 13.0 99Y25D (TAT→GAT) aroH3‑deoxy‑D‑arabino‑heptulosonate‑7‑phosphate synthase, tryptophan repressible
Reads supporting (aligned to +/- strand):  ref base T (44/49);  new base G (3/3);  total (47/52)
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.02e-01

GCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGTACGCCA  >  NC_000913/1788361‑1788651
                                                                                                                                                  |                                                                                                                                                
gcGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTAt                                                                                                                                                >  1:1063441/1‑149 (MQ=255)
  gATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATcc                                                                                                                                              <  4:1363409/149‑1 (MQ=255)
   aTTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATccc                                                                                                                                             <  4:303941/149‑1 (MQ=255)
        tttATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTa                                                                                                                                          <  1:1035539/147‑1 (MQ=255)
        tttATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTa                                                                                                                                          >  2:1035539/1‑147 (MQ=255)
        tttATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAAc                                                                                                                                        >  2:662725/1‑149 (MQ=255)
        tttATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAAc                                                                                                                                        >  3:332089/1‑149 (MQ=255)
           aTCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGTTACGGGATACCGTAgcgcg                                                                                                                                     >  4:591374/1‑148 (MQ=255)
           aTCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGcc                                                                                                                                     >  2:72024/1‑149 (MQ=255)
            tCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGTTACGGGATACCGTAgcgcga                                                                                                                                    <  3:591374/149‑3 (MQ=255)
            tCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCt                                                                                                                                    <  2:333070/149‑1 (MQ=255)
            tCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCt                                                                                                                                    <  2:1278160/149‑1 (MQ=255)
             cGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTg                                                                                                                                   <  1:1055695/149‑1 (MQ=255)
                 aTGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCt                                                                                                                                    <  4:501499/144‑1 (MQ=255)
                 aTGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCt                                                                                                                                    >  3:501499/1‑144 (MQ=255)
                       aTCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATc                                                                                                                                               <  1:990500/127‑1 (MQ=255)
                       aTCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATc                                                                                                                                               >  2:990500/1‑127 (MQ=255)
                           cTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCAt                                                                                                                     <  1:194079/149‑1 (MQ=255)
                            ttCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGTGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATg                                                                                                                    >  3:524491/1‑149 (MQ=255)
                            ttCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTaa                                                                                                                                         >  4:131770/1‑128 (MQ=255)
                            ttCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTaa                                                                                                                                         <  3:131770/128‑1 (MQ=255)
                            ttCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCa                                                                                                                          <  4:578553/143‑1 (MQ=255)
                            ttCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCa                                                                                                                          >  3:578553/1‑143 (MQ=255)
                            ttCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATg                                                                                                                    >  4:1274369/1‑149 (MQ=255)
                             tCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATcc                                                                                                                                              >  3:322459/1‑122 (MQ=255)
                             tCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATcc                                                                                                                                              <  4:322459/122‑1 (MQ=255)
                               ccGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCa                                                                                                                 <  1:662725/149‑1 (MQ=255)
                                cGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGTGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCAc                                                                                                                <  4:524491/149‑1 (MQ=255)
                                             gTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCcgccgc                                                                                                   >  2:355803/1‑149 (MQ=255)
                                               aGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAg                                                                                                 >  3:354165/1‑149 (MQ=255)
                                                  aGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAAt                                                                                              <  1:205144/149‑1 (MQ=255)
                                                  aGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAAt                                                                                              >  2:205144/1‑149 (MQ=255)
                                                     aacaaATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGa                                                                                           <  1:355803/149‑1 (MQ=255)
                                                      acaaATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGaa                                                                                          >  2:639199/1‑149 (MQ=255)
                                                      acaaATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGaa                                                                                          >  4:1167371/1‑149 (MQ=255)
                                                      acaaATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGaa                                                                                          >  4:985992/1‑149 (MQ=255)
                                                       caaATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGc                                                                                                                            <  4:1395512/114‑1 (MQ=255)
                                                       caaATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGc                                                                                                                            >  3:1395512/1‑114 (MQ=255)
                                                          aTTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGaaaaaa                                                                                      >  2:821477/1‑149 (MQ=255)
                                                             tCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATAc                                                                                   <  1:639199/149‑1 (MQ=255)
                                                                gagaCTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGa                                                                                >  3:656371/1‑149 (MQ=255)
                                                                 agccTTGTAATGAACAGAACTGACGACCGCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGaa                                                                               <  4:332089/146‑1 (MQ=255)
                                                                   aCTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATg                                                                             <  4:656371/149‑1 (MQ=255)
                                                                    cTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATgg                                                                            <  3:1274369/149‑1 (MQ=255)
                                                                     ttGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCGTGGc                                                                                                                                 >  1:801090/1‑95 (MQ=255)
                                                                     ttGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGc                                                                                                                                 <  2:801090/95‑1 (MQ=255)
                                                                        tAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGaa                                                                        <  3:1167371/149‑1 (MQ=255)
                                                                          aTGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATa                                                                                    >  2:761174/1‑135 (MQ=255)
                                                                          aTGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATa                                                                                    <  1:761174/135‑1 (MQ=255)
                                                                          aTGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGa                                                                      <  2:1063441/149‑1 (MQ=255)
                                                                          aTGAACACAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGCCTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGa                                                                      <  1:821477/149‑1 (MQ=255)
                                                                              aCAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATg                                                                             >  2:311344/1‑138 (MQ=255)
                                                                              aCAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATg                                                                             <  1:311344/138‑1 (MQ=255)
                                                                               cAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGc                                                                 <  1:72024/149‑1 (MQ=255)
                                                                                 gAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGa                                                               <  4:354165/149‑1 (MQ=255)
                                                                                  accTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAAt                                                                                              <  2:466566/115‑1 (MQ=255)
                                                                                  aaCTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAAt                                                                                              >  1:466566/1‑117 (MQ=255)
                                                                                      gACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGa                                                                                                             >  3:114552/1‑98 (MQ=255)
                                                                                      gACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGa                                                                                                             <  4:114552/98‑1 (MQ=255)
                                                                                       aCGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAg                                                                       <  3:14860/135‑1 (MQ=255)
                                                                                       aCGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAg                                                                       >  4:14860/1‑135 (MQ=255)
                                                                                       aCGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTg                                                         <  1:1152039/149‑1 (MQ=255)
                                                                                       aCGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACAGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTg                                                         >  1:251121/1‑149 (MQ=255)
                                                                                         gAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGt                                                       >  1:886294/1‑149 (MQ=255)
                                                                                           aCTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGaa                                                                        <  3:105089/130‑1 (MQ=255)
                                                                                           aCTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGaa                                                                        >  2:835223/1‑130 (MQ=255)
                                                                                           aCTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGaa                                                                        >  4:105089/1‑130 (MQ=255)
                                                                                           aCTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGaa                                                                        <  1:835223/130‑1 (MQ=255)
                                                                                            cTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGaa                                                                               <  1:1078218/122‑1 (MQ=255)
                                                                                            cTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGaa                                                                               >  2:1078218/1‑122 (MQ=255)
                                                                                             tCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCatt                                                   >  1:1208434/1‑149 (MQ=255)
                                                                                                gTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCattatt                                                <  2:1208434/149‑1 (MQ=255)
                                                                                                       gcgTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGccc                                           >  2:1253357/1‑147 (MQ=255)
                                                                                                       gcgTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCt                                         >  3:737762/1‑149 (MQ=255)
                                                                                                       gcgTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCAGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGccc                                           <  1:1253357/147‑1 (MQ=255)
                                                                                                          tATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCt                                      <  2:886294/149‑1 (MQ=255)
                                                                                                                  gCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCAc                              <  2:251121/149‑1 (MQ=255)
                                                                                                                  gCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCAc                              >  2:863773/1‑149 (MQ=255)
                                                                                                                     tGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTc                                                      <  3:823692/122‑1 (MQ=255)
                                                                                                                     tGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTc                                                      >  4:823692/1‑122 (MQ=255)
                                                                                                                       gTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGAtct                         <  3:985992/149‑1 (MQ=255)
                                                                                                                          aCGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGaa                                                                               >  3:83669/1‑92 (MQ=255)
                                                                                                                          aCGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGaa                                                                               <  4:83669/92‑1 (MQ=255)
                                                                                                                                  cgAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCaa              <  1:863773/149‑1 (MQ=255)
                                                                                                                                   gAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATAc                                                                                   >  4:328899/1‑79 (MQ=255)
                                                                                                                                   gAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATAc                                                                                   <  3:328899/79‑1 (MQ=255)
                                                                                                                                       tCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTTCTGGTCatta                                                  >  1:513181/1‑108 (MQ=255)
                                                                                                                                       tCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCatta                                                  <  2:513181/108‑1 (MQ=255)
                                                                                                                                        ggcgtTACGGGATACCGTAGCGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCAc                              <  2:979508/126‑1 (MQ=255)
                                                                                                                                        ggcgtTACGGGATACCGTAGCGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCAc                              >  1:979508/2‑127 (MQ=255)
                                                                                                                                         gcgtTACGGGATACCGTAGCGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCAc                              >  2:1358858/1‑126 (MQ=255)
                                                                                                                                         gcgtTACGGGATACCGTAGCGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCAc                              <  1:1358858/126‑1 (MQ=255)
                                                                                                                                         gcgcTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATa                                                                                    >  4:140588/1‑72 (MQ=255)
                                                                                                                                         gcgcTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATa                                                                                    <  3:140588/72‑1 (MQ=255)
                                                                                                                                         gcgcTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCatt                                                   <  1:579460/105‑1 (MQ=255)
                                                                                                                                         gcgcTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCatt                                                   >  2:579460/1‑105 (MQ=255)
                                                                                                                                         gcgcTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGAt                                  <  2:710884/122‑1 (MQ=255)
                                                                                                                                         gcgcTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGAt                                  >  1:710884/1‑122 (MQ=255)
                                                                                                                                              aCGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGTACGCCa  >  3:241197/1‑149 (MQ=255)
                                                                                                                                                  |                                                                                                                                                
GCGATTGGTTTATCGTGATGCGCATCACTTCCCGGCAGTCCTGCCGTAGAAGCAACAAATTTCTGAGACTTGTAATGAACAGAACTGACGAACTCCGTACTGCGCGTATTGAGAGCCTGGTAACGCCCGCCGAACTCGCGCTACGGTATCCCGTAACGCCTGGCGTCGCCACCCATGTCACCGACTCCCGCCGCAGAATTGAAAAAATACTGAATGGTGAAGATAAGCGACTGTTGGTCATTATTGGCCCCTGCTCGATCCACGATCTCACCGCTGCAATGGAGTACGCCA  >  NC_000913/1788361‑1788651

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