Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,008,423 C→T 5.8% intergenic (+46/+16) 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,4230CT 5.8% 271.0 / 4.6 88intergenic (+46/+16)rhtC/rhtBthreonine efflux pump/homoserine, homoserine lactone and S‑methyl‑methionine efflux pump
Reads supporting (aligned to +/- strand):  ref base C (46/34);  new base T (3/2);  total (50/38)
Fisher's exact test for biased strand distribution p-value = 1.00e+00
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 9.30e-01

ATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAGGGTGGCGTAAC  >  NC_000913/4008282‑4008569
                                                                                                                                             |                                                                                                                                                  
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:597576/1‑149 (MQ=255)
                                                             ccGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGc                                                                                <  2:1455616/149‑1 (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:208904/1‑149 (MQ=255)
                                                                        attcattTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTa                                                                     >  2:823473/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                                                    >  1:1012627/1‑149 (MQ=255)
                                                                                         gcgGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTg                                                    >  2:998015/1‑149 (MQ=255)
                                                                                              gATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTg                                               <  2:292579/149‑1 (MQ=255)
                                                                                               aTGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGCGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTgg                                              >  2:333733/1‑149 (MQ=255)
                                                                                                tGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGt                                             <  1:1206151/149‑1 (MQ=255)
                                                                                                tGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGt                                             >  2:1549670/1‑149 (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                                       >  1:1072060/1‑149 (MQ=255)
                                                                                                      aCGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTa                                       <  2:758709/149‑1 (MQ=255)
                                                                                                      aCGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTa                                       <  2:1340839/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:731865/149‑1 (MQ=255)
                                                                                                               cAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATaga                              <  2:703854/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                     >  2:438134/1‑149 (MQ=255)
                                                                                                                        tCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGt                     >  1:1269857/1‑149 (MQ=255)
                                                                                                                           gatgcggCGTAAACGCCTTATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGa                  <  2:1246174/144‑1 (MQ=255)
                                                                                                                            atgcggCGTAAACGCCTTATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAATCAACGAGACGAAAATCTTATTCAGCGCCTTCACCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAg                 >  1:232189/5‑149 (MQ=255)
                                                                                                                            aTAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAg                 <  2:1066122/149‑1 (MQ=255)
                                                                                                                                ggCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAg             <  1:1121828/149‑1 (MQ=255)
                                                                                                                                   gTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAGGGt          <  2:824756/149‑1 (MQ=255)
                                                                                                                                   gTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAGGGt          >  2:1555360/1‑149 (MQ=255)
                                                                                                                                   gTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAGGGt          <  2:1567704/149‑1 (MQ=255)
                                                                                                                                      tACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAGGGTGGc       >  1:1344127/1‑149 (MQ=255)
                                                                                                                                      tACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAGGGTGGc       <  2:1269857/149‑1 (MQ=255)
                                                                                                                                          ccgcATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAGGGTGGCGTaa   <  1:320760/149‑1 (MQ=255)
                                                                                                                                          ccgcATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAGGGTGGCGTaa   >  1:725313/1‑149 (MQ=255)
                                                                                                                                           cgcATCCGACATTATTTTTCACGCATGCCTCGCTGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAGGGTGGCGTAAc  <  1:1254219/149‑1 (MQ=255)
                                                                                                                                             |                                                                                                                                                  
ATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAGGGTGGCGTAAC  >  NC_000913/4008282‑4008569

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: