Predicted mutation
evidence seq id position mutation freq annotation gene description
RA NC_000913 4,008,415 T→A 9.1% intergenic (+38/+24) rhtC → / ← rhtB L‑threonine exporter/L‑homoserine/L‑homoserine lactone/L‑threonine exporter

Read alignment evidence...
  seq id position ref new freq score (cons/poly) reads annotation genes product
*NC_0009134,008,4150TA 9.1% 352.6 / 20.1 121intergenic (+38/+24)rhtC/rhtBL‑threonine exporter/L‑homoserine/L‑homoserine lactone/L‑threonine exporter
Reads supporting (aligned to +/- strand):  ref base T (55/54);  new base A (3/8);  total (59/62)
Fisher's exact test for biased strand distribution p-value = 2.07e-01
Kolmogorov-Smirnov test that lower quality scores support variant p-value = 4.76e-01

Only 100 of 123 total aligned reads displayed.

GTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAGGGTG  >  NC_000913/4008266‑4008562
                                                                                                                                                     |                                                                                                                                                   
gTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGtt                                                                                                                                                     <  2:272001/150‑1 (MQ=255)
gTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGtt                                                                                                                                                     >  2:395624/1‑150 (MQ=255)
     cccTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTAcgc                                                                                                                                                >  2:1322630/1‑150 (MQ=255)
       cTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTAcgccg                                                                                                                                              >  2:108404/1‑150 (MQ=255)
         gCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCa                                                                                                                                            >  2:1350832/1‑150 (MQ=255)
            gCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                         <  1:252451/150‑1 (MQ=255)
            gCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATcc                                                                                                                                         >  1:1029344/1‑150 (MQ=255)
            gCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGAGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATGCGGCGTAAACGCCTTATcc                                                                                                                                         <  2:294137/150‑1 (MQ=255)
               aaTGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGAc                                                                                                                                      <  1:1134230/150‑1 (MQ=255)
                 tGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATGCGGCGTAAACGCCTTATCCGActt                                                                                                                                    <  1:127782/150‑3 (MQ=255)
                 tGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATGCGGCGTAAACGCCTTATCCGActt                                                                                                                                    <  1:1575060/150‑3 (MQ=255)
                 tGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATGCGGCGTAAACGCCTTATCCGActt                                                                                                                                    <  2:842032/150‑3 (MQ=255)
                 tGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACat                                                                                                                                    <  1:248611/150‑1 (MQ=255)
                  gcgcCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATGCGGCGTAAACGCCTTATCCGActta                                                                                                                                   <  1:855248/150‑4 (MQ=255)
                   cgcCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACatta                                                                                                                                  >  2:838127/1‑150 (MQ=255)
                     ccGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACattatt                                                                                                                                >  2:1136128/1‑150 (MQ=255)
                       gTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTAtttt                                                                                                                              <  1:1246371/150‑1 (MQ=255)
                         ggTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTAttttta                                                                                                                            >  2:1422786/1‑149 (MQ=255)
                          gTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCa                                                                                                                           <  1:71871/150‑1 (MQ=255)
                          gTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCa                                                                                                                           <  1:1458567/150‑1 (MQ=255)
                           ttATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGGTTACGCCGCATCCGACATTATTTTTCAc                                                                                                                          <  2:621427/150‑1 (MQ=255)
                             aTCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGc                                                                                                                        <  2:109272/150‑1 (MQ=255)
                                  cGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGcc                                                                                                                   <  1:371972/150‑1 (MQ=255)
                                      tGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGc                                                                                                               <  1:786128/150‑1 (MQ=255)
                                       ggCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGCCGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGcc                                                                                                              >  2:519633/1‑150 (MQ=255)
                                       ggCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGcc                                                                                                              >  2:159700/1‑150 (MQ=255)
                                       ggCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGcc                                                                                                              <  1:976370/150‑1 (MQ=255)
                                        gCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCg                                                                                                             >  2:839323/1‑150 (MQ=255)
                                          gAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGAt                                                                                                           <  1:1317239/150‑1 (MQ=255)
                                            aGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCAGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGc                                                                                                         >  2:750279/1‑150 (MQ=255)
                                             tttgattgattTTTTTTCCGGGGGGTTATTTGACGGATTTGTCATTAATTTTATTATTTCGCGGTGATGCCAGAGGCGTCTTCAGAGTAAGTCTGATAAGTAGTTTACGCCGGATCCTACATTATTTTTCAGGGATGCCTCGCCGATGCt                                                                                                        <  2:740946/147‑1 (MQ=255)
                                             gTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCGAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCt                                                                                                        <  2:735153/150‑1 (MQ=255)
                                             gTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCt                                                                                                        <  1:1320558/150‑1 (MQ=255)
                                               ggATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTaa                                                                                                      >  1:1338789/1‑150 (MQ=255)
                                                    gATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAgcg                                                                                                 >  2:1285814/1‑150 (MQ=255)
                                                    gATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAgcg                                                                                                 >  1:712587/1‑150 (MQ=255)
                                                    gATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAgcg                                                                                                 >  2:1472526/1‑150 (MQ=255)
                                                     aTGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAgcgc                                                                                                >  1:1389529/1‑150 (MQ=255)
                                                      tGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGCTGCTAACAGCGCt                                                                                               >  2:32310/1‑150 (MQ=255)
                                                       ggTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTc                                                                                              <  1:326156/150‑1 (MQ=255)
                                                        gTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTcc                                                                                             <  1:1025583/150‑1 (MQ=255)
                                                            tGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCccacc                                                                                         <  1:1409676/150‑1 (MQ=255)
                                                            tGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCccacc                                                                                         >  1:928457/1‑150 (MQ=255)
                                                            tGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCccacc                                                                                         <  1:1136128/150‑1 (MQ=255)
                                                              ccGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAg                                                                                       <  2:568360/150‑1 (MQ=255)
                                                                ggggCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCa                                                                                     <  1:1565430/150‑1 (MQ=255)
                                                                  ggCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATa                                                                                   >  1:909309/1‑150 (MQ=255)
                                                                  ggCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATa                                                                                   <  2:330235/150‑1 (MQ=255)
                                                                    cGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATaaa                                                                                 <  2:1564505/150‑1 (MQ=255)
                                                                         tttGCCTGTTTTTGCATTCATTTTATTATTTCGCGGTGATTCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAaacaaag                                                                            <  1:1399708/150‑3 (MQ=255)
                                                                          ttGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGa                                                                           >  1:1087597/1‑150 (MQ=255)
                                                                           tGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAg                                                                          >  1:1562933/1‑150 (MQ=255)
                                                                            gCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGc                                                                         <  2:909309/150‑1 (MQ=255)
                                                                             ccGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGcc                                                                        >  1:331839/1‑150 (MQ=255)
                                                                                    tGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATc                                                                 >  1:1258682/1‑150 (MQ=255)
                                                                                     ggCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCt                                                                <  1:814713/150‑1 (MQ=255)
                                                                                     ggCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCt                                                                >  2:848426/1‑150 (MQ=255)
                                                                                       cattcattTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTa                                                              >  2:31258/1‑150 (MQ=255)
                                                                                          tcattTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTc                                                           <  1:337441/150‑1 (MQ=255)
                                                                                           cattTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCa                                                          >  2:785977/1‑150 (MQ=255)
                                                                                            attTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAgcg                                                                                                 >  2:452904/1‑110 (MQ=255)
                                                                                            attTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAgcg                                                                                                 <  1:452904/110‑1 (MQ=255)
                                                                                              ttGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCGCAGCAGCATAAACAACGTGCCGAAAATCTTATTCAgcg                                                       <  2:118682/150‑1 (MQ=255)
                                                                                              ttGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAgcg                                                       <  1:724369/150‑1 (MQ=255)
                                                                                               tGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAgcgc                                                      >  2:775146/1‑150 (MQ=255)
                                                                                                gATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGcc                                                     >  1:1095758/1‑150 (MQ=255)
                                                                                                   tattTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTAGGATAAGGCGATTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGATAATCTTATTCAGCGCCTTc                                                  <  2:741884/150‑1 (MQ=255)
                                                                                                    attTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCa                                                 >  1:763917/1‑150 (MQ=255)
                                                                                                       tCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCt                                              <  2:1029344/150‑1 (MQ=255)
                                                                                                         gcgGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGc                                            >  1:840513/1‑150 (MQ=255)
                                                                                                            gTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCttt                                         <  2:1093369/150‑1 (MQ=255)
                                                                                                              gATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTgg                                       <  1:1350832/150‑1 (MQ=255)
                                                                                                                tGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTc                                     <  1:1050806/150‑1 (MQ=255)
                                                                                                                 gcgaGACGCTTCTTCAGAGTAAGTCTGATAAGGCGTTTACGCCGAATCCGACTTTATTTTTCACGCATGCCTCGCCGATGCTAACAGCTCTCCCACCAGCATAAACAACTAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTcc                                    <  2:1258682/147‑1 (MQ=255)
                                                                                                                 gCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCGTTGGTcc                                    >  2:1452212/1‑150 (MQ=255)
                                                                                                                   cAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCtt                                  <  2:712587/150‑1 (MQ=255)
                                                                                                                   cAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCtt                                  <  1:960625/150‑1 (MQ=255)
                                                                                                                    aGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCttt                                 <  1:218224/150‑1 (MQ=255)
                                                                                                                    aGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCttt                                 <  1:839323/150‑1 (MQ=255)
                                                                                                                     gACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACACCGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTa                                >  1:331362/1‑150 (MQ=255)
                                                                                                                      aCGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTaa                               >  2:1175077/1‑150 (MQ=255)
                                                                                                                        gcgTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATc                             <  1:1422786/150‑1 (MQ=255)
                                                                                                                         cgTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATcc                            >  2:1486830/1‑150 (MQ=255)
                                                                                                                              tCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATaga                       >  1:1601667/1‑150 (MQ=255)
                                                                                                                               cAGAGTAAGTCGGTTTAGCCGTGTACGCCGCATACGAGAGTATTTTTCACGCATGACACGCCTATGCTAACACCGCTCCCACCAGCATAAACCACGAGGCGAAAATCTTATTGAGCGCCTTAATCTGCTTTGGTCCTTTAATCCagagag                      >  2:811274/1‑150 (MQ=255)
                                                                                                                                    aatgTCGGATGCGGCGTAAACGCCTTATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATc                 >  1:1175078/4‑150 (MQ=255)
                                                                                                                                     atgTCGGATGCGGCGTAAACGCCTTATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCACCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCATTGGTCCTTTAATCCATAGAGCAATcc                >  2:1445513/3‑150 (MQ=255)
                                                                                                                                       gTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGt              >  1:179045/1‑150 (MQ=255)
                                                                                                                                       gTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGt              <  2:1389529/150‑1 (MQ=255)
                                                                                                                                         cGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTg            >  1:451790/1‑150 (MQ=255)
                                                                                                                                         cGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCAGGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTg            <  2:251893/150‑1 (MQ=255)
                                                                                                                                          ggATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGGCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCGGTTAATCCATAGAGCAAGCCGTTGa           >  2:1141712/1‑150 (MQ=255)
                                                                                                                                            aTAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGc         <  1:1472526/150‑1 (MQ=255)
                                                                                                                                             tAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCa        >  1:744234/1‑150 (MQ=255)
                                                                                                                                              aaGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCaa       >  1:315771/1‑150 (MQ=255)
                                                                                                                                               cggCGTAAACGCCTTATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAAt                  >  2:1526264/2‑138 (MQ=255)
                                                                                                                                                ggCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAgg     >  1:1022587/1‑150 (MQ=255)
                                                                                                                                                   gTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTg                                             >  1:41935/1‑107 (MQ=255)
                                                                                                                                                   gTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAGGGTg  >  1:78364/1‑150 (MQ=255)
                                                                                                                                                   gTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTAATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAGGGTg  >  1:292074/1‑150 (MQ=255)
                                                                                                                                                     |                                                                                                                                                   
GTTTGCCCTGCCGCAAATGCGCCGTGGTTATCAACGTCTGGCGAAGTGGATTGATGGTTTTGCCGGGGCGTTATTTGCCGGATTTGGCATTCATTTGATTATTTCGCGGTGATGCCAGACGCGTCTTCAGAGTAAGTCGGATAAGGCGTTTACGCCGCATCCGACATTATTTTTCACGCATGCCTCGCCGATGCTAACAGCGCTCCCACCAGCATAAACAACGAGCCGAAAATCTTATTCAGCGCCTTCATCTGCTTTGGTCCTTTAATCCATAGAGCAATCCGTTGAGCAAGGGTG  >  NC_000913/4008266‑4008562

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