Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,008,409 A→C 5.8% intergenic (+32/+30) rhtC → / ← rhtB threonine efflux pump/homoserine, homoserine lactone and S‑methyl‑methionine efflux pump

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,008,4090AC 5.8% 257.4 / 7.0 87intergenic (+32/+30)rhtC/rhtBthreonine efflux pump/homoserine, homoserine lactone and S‑methyl‑methionine efflux pump
Reads supporting (aligned to +/- strand):  ref base A (48/34);  new base C (2/3);  total (50/37)
Fisher's exact test for biased strand distribution p-value = 6.47e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 3.88e-01

GTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAG  >  NC_000913/4008266‑4008554
                                                                                                                                               |                                                                                                                                                 
gTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGt                                                                                                                                              >  1:611739/1‑149 (MQ=255)
gTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGt                                                                                                                                              <  1:1461881/149‑1 (MQ=255)
 tttGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGtt                                                                                                                                             >  2:1051099/1‑149 (MQ=255)
     cccTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTAca                                                                                                                                         >  2:889991/1‑148 (MQ=255)
      ccTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTAcgc                                                                                                                                        <  2:219230/149‑1 (MQ=255)
            gCAAATGCGCCGTGGTTATGAACGTCGGGCGCAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCTGATGCGGCGTAAAcgccttatc                                                                                                                                  <  2:782041/149‑6 (MQ=255)
            gCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATGCGGCGTAAAcgccttatc                                                                                                                                  <  2:1547889/149‑6 (MQ=255)
                aTGCGCCGGGGTTAGCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATGCGGCGTAAACGCCTTATCCGAc                                                                                                                              <  1:279342/149‑1 (MQ=255)
                     ccGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATc                                                                                                                                  <  2:747804/140‑1 (MQ=255)
                     ccGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATc                                                                                                                                  >  1:747804/1‑140 (MQ=255)
                        gggTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTAtttt                                                                                                                      <  1:356000/148‑1 (MQ=255)
                          gTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTc                                                                                                                    >  2:681626/1‑149 (MQ=255)
                          gTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTc                                                                                                                    >  1:758709/1‑149 (MQ=255)
                           ttATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCa                                                                                                                   <  2:242049/149‑1 (MQ=255)
                             ataaaCGTATGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGAAGCTTATACAGAATAATTCGTATAAGGCTTTTAAGCCGAATCCTACATTATTTTTAAca                                                                                                                 >  2:1576459/4‑148 (MQ=255)
                               cAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCa                                                                                                               >  2:1254219/1‑149 (MQ=255)
                                   gTCTGGCGAAGGGGATTGATGGTTTTGCCGGGCCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGCTGCGAGCCGCGTCTTCATAGTAAGTCGGATAAGGCGTTTACGCCGCCTCCGACATTATTTTACACGCAAGca                                                                                                           >  1:259303/1‑148 (MQ=255)
                                          gAAGTGCATTGATTGTTTAGCAGGGGCGTAATTTTCCTGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGa                                                                                                    <  2:793098/149‑1 (MQ=255)
                                           aaGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGAt                                                                                                   >  1:820502/1‑149 (MQ=255)
                                            aGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATg                                                                                                  <  1:656363/149‑1 (MQ=255)
                                                gATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTaa                                                                                              >  1:1507981/1‑149 (MQ=255)
                                                                      ttATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAaac                                                                        >  1:292579/1‑149 (MQ=255)
                                                                        aTTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGTCTCGCCGATGCTAACAGCGCTCCCACCAGCATAaacaa                                                                      <  1:734805/149‑1 (MQ=255)
                                                                        aTTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAaacaa                                                                      >  1:856934/1‑149 (MQ=255)
                                                                           tGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGAAGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTAAGCCGCATCCGACATTATTTTTCTCGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGa                                                                   <  2:772665/149‑1 (MQ=255)
                                                                            gCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAg                                                                  >  1:582966/1‑149 (MQ=255)
                                                                             ccGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGc                                                                 <  2:1455616/149‑1 (MQ=255)
                                                                             ccGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGc                                                                 >  2:597576/1‑149 (MQ=255)
                                                                                gATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGa                                                              >  1:1066122/1‑149 (MQ=255)
                                                                                  tttGGCATTCATTTGAATATATCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTAAGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGaaa                                                            >  2:933312/1‑149 (MQ=255)
                                                                                    tGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAAt                                                          >  2:544758/1‑149 (MQ=255)
                                                                                     ggCATTCATTTGCTGATATCGAGGTGATGCCCGACGCTGCTGCCGCGGAAGTCGGAGAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATc                                                         <  2:811130/149‑1 (MQ=255)
                                                                                     ggCATTCATTTGATTATTTCGAGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGAATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGATAATc                                                         <  1:664182/149‑1 (MQ=255)
                                                                                      gCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCTACATTATTTTTCAAGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCt                                                        >  1:1058843/1‑149 (MQ=255)
                                                                                      gCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCt                                                        >  2:320760/1‑149 (MQ=255)
                                                                                      gCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCt                                                        >  2:680335/1‑149 (MQ=255)
                                                                                        attcattTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCTGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTa                                                      <  2:76140/149‑1 (MQ=255)
                                                                                        attcattTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTa                                                      >  2:823473/1‑149 (MQ=255)
                                                                                        attcattTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTa                                                      >  2:208904/1‑149 (MQ=255)
                                                                                        attcattTGATTATTTCGCGGGGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTa                                                      <  2:359430/149‑1 (MQ=255)
                                                                                         ttcattTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTAt                                                     >  2:375788/1‑149 (MQ=255)
                                                                                          tcattTGATTATTTCGCGGTGAGGCCAGACGCGTCTTCAGAGTAAGGCGGATAAGGCGTTTACGCCGCATCCGACATTATTATTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTAtt                                                    <  1:1504088/149‑1 (MQ=255)
                                                                                           aaatttcttattTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTc                                                   <  1:583195/142‑1 (MQ=255)
                                                                                              ttGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGc                                                                                                       <  2:296307/94‑1 (MQ=255)
                                                                                              ttGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGc                                                                                                       >  1:296307/1‑94 (MQ=255)
                                                                                                 attattTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGcc                                             >  2:1287412/1‑149 (MQ=255)
                                                                                                  ttattTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCt                                            >  2:1121828/1‑149 (MQ=255)
                                                                                                  ttattTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCt                                            >  1:1113967/1‑149 (MQ=255)
                                                                                                   tattTCGCGGTGATGCCAGACGCGGCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCtt                                           <  2:1073591/149‑1 (MQ=255)
                                                                                                       tCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGTCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATc                                       <  1:288609/149‑1 (MQ=255)
                                                                                                       tCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATc                                       >  2:470587/1‑149 (MQ=255)
                                                                                                       tCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATc                                       <  1:1290078/149‑1 (MQ=255)
                                                                                                         gcgGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTg                                     >  2:998015/1‑149 (MQ=255)
                                                                                                         gcgGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTg                                     >  1:1012627/1‑149 (MQ=255)
                                                                                                              gATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTg                                <  2:292579/149‑1 (MQ=255)
                                                                                                               aTGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGCGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTgg                               >  2:333733/1‑149 (MQ=255)
                                                                                                                tGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGt                              >  2:1549670/1‑149 (MQ=255)
                                                                                                                tGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGt                              <  1:1206151/149‑1 (MQ=255)
                                                                                                                 gCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGACTCGCCGATGATAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTc                             <  1:74249/149‑1 (MQ=255)
                                                                                                                  ccAGCCGCGTCTTCAGATGAAGTCGGATAAGGCGTTTACGCCGAATCCGACATTATTTTTCACGCAGTCCTCGGCGATGCTAACAGCGCTCCAACCAGGATAAACTACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTcc                            <  2:90024/149‑1 (MQ=255)
                                                                                                                   cAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCCCCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCt                           >  1:151195/1‑149 (MQ=255)
                                                                                                                    aGACGCTTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCTCATCCGACATTATTTTTCACGCATGACTCGCCTATGCTAACAGCGCTCCCACCAGCTTAAACAACGAGCAGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCtt                          <  1:1576459/149‑1 (MQ=255)
                                                                                                                    aGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCCCCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCtt                          <  1:877326/149‑1 (MQ=255)
                                                                                                                     gACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCttt                         >  2:742901/1‑149 (MQ=255)
                                                                                                                      aCGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTa                        <  2:1340839/149‑1 (MQ=255)
                                                                                                                      aCGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTa                        >  1:1072060/1‑149 (MQ=255)
                                                                                                                      aCGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTa                        <  2:758709/149‑1 (MQ=255)
                                                                                                                      aCGCGTCTACAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCAAGCCTCGCCGATGCAAACAGCGCTCCAACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGACCTTTa                        >  1:1046510/1‑149 (MQ=255)
                                                                                                                       cgcgTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTaa                       <  1:1448910/149‑1 (MQ=255)
                                                                                                                        gcgTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCTCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTTAt                      >  2:721013/1‑149 (MQ=255)
                                                                                                                        gcgTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAAt                      >  1:1567704/1‑149 (MQ=255)
                                                                                                                            cTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCAt                  >  1:1446313/1‑149 (MQ=255)
                                                                                                                             ttCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATa                 <  2:820502/149‑1 (MQ=255)
                                                                                                                               cAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATaga               <  2:703854/149‑1 (MQ=255)
                                                                                                                               cAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATaga               <  2:731865/149‑1 (MQ=255)
                                                                                                                               cAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATaga               >  1:703854/1‑149 (MQ=255)
                                                                                                                                agagTAAGTCGGATAAGGCGATTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGTGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATagag              >  2:1271371/1‑149 (MQ=255)
                                                                                                                                  agTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCAGCTGCTTTGGTCCTTTAATCCATAGAGCa            >  1:614680/1‑149 (MQ=255)
                                                                                                                                   taatgTCGGATGCGGCGTAAACGCCTTATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCaa           >  1:1246174/5‑149 (MQ=255)
                                                                                                                                    tAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAAt          <  2:611739/149‑1 (MQ=255)
                                                                                                                                      aGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAAAAACGAGCCGAAAATCTTATTCAGCGCCTACATCTGCTTTGGTCCTTTAATCCATAGAGCAATcc        >  1:207118/1‑149 (MQ=255)
                                                                                                                                       gTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCg       >  2:510849/1‑149 (MQ=255)
                                                                                                                                       gTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCg       >  1:837005/1‑149 (MQ=255)
                                                                                                                                        tCGGATGCGGCGTAAACGCCTTATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGt      >  1:1320111/1‑149 (MQ=255)
                                                                                                                                        tCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGt      >  1:1269857/1‑149 (MQ=255)
                                                                                                                                        tCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGt      >  2:438134/1‑149 (MQ=255)
                                                                                                                                            aTAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAg  <  2:1066122/149‑1 (MQ=255)
                                                                                                                                               |                                                                                                                                                 
GTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAG  >  NC_000913/4008266‑4008554

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: